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Nitride Fuels: Properties and Potentials 


By Arthur A. Bauer* 


Abstract: Interest in nitride fuels centers principally on the use 
of (U,Pu)N in fast breeder reactors, A review of nitride fuel 
properties is presented including properties and irradiation 
behavior of UN, whose major potential application is as a space 
reactor fuel. Information on UN as presented here serves 
principally as background to, and to supplement, data on 
mixed-nitride, (U,Pu)N, fuels. 

Although the mixed nitrides have received less con- 
centrated study, a sufficient technology base appears to exist 
that will permit an evaluation of mixed nitrides as liquid- 
metal-cooled fast breeder reactor fuels. The properties of 
(U,Pu)N related to its irradiation performance, safety char- 
acteristics, and fuel-cycle costs are reviewed, Fuel-pin design is 
discussed, and fuel-performance potential and limitations are 
outlined, Emphasis is placed on potential advantages resulting 
from the use of (U,Pu)N to fuel an advanced liquid-metal- 
cooled fast breeder reactor. 


The search for high-performance uranium and plu- 
tonium fuels as alternates to the oxides has con- 
centrated on the nitrides and carbides. Oxides, because 
of their compatibility with the coolant, are uniquely 
suited for use in water reactors. It is only with other 
coolants that nitrides and carbides come into serious 
consideration. 

There is a general tendency to refer to nitrides and 
carbides interchangeably because of their similarities in 
properties and behavior. However, there are significant 
differences between the two that eventually will 
determine which is chosen for exploitation. It is not 
the intention in this article to compare these fuels, 
although comparisons at times will be unavoidable, but 
to review the available data concerning nitride fuels 
which permit an evaluation of their commercial po- 
tential. Since this potential is most evident in a 
liquid-metal-cooled fast breeder reactor (LMFBR) 
application, the review will focus on mixed-nitride, 
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(U,Pu)N, fuels. However, interest in uranium mononi- 
tride, UN, for high-surface-temperature operation in 
space power reactors remains active. Consequently UN 
data also will be reviewed, although principally for the 
purpose of providing background to the evaluation of 
the mixed-nitride fuels. 


PHYSICAL PROPERTIES 





Constitution and Phase Stability 


Compounds associated with UN, U,N3, and UN, 
compositions have been identified in the uranium— 
nitrogen system.'~* UN, the most stable and the only 
one of interest as a fuel material, has a face-centered 
cubic (fcc) NaCl structure*® and a theoretical density 
of 14.32 g/cm*. U,N3 has been found to exist in two 
forms, the body-centered cubic (bcc) alpha and the 
hexagonal beta structure.2*7'° The alpha structure, 
the lower temperature form, exists at compositions 
slightly hyperstoichiometric with respect to U,N3 and 
decomposes peritectoidally, with loss of nitrogen, to 
the beta form, which is slightly deficient in nitrogen 
with respect to U,N3. Under 1 atm nitrogen the 
alpha-to-beta decomposition on heating occurs at 
about 1130°C, whereas the beta decomposes to form 
UN with the evolution of nitrogen at about 1350°C. 
Alpha U,;N3 forms a complete series of solid solutions 
with UN, which has an fcc UO, fluorite structure; 
this structure is observed at substoichiometric UN, 
compositions, 

The nitrogen decomposition pressure of UN has 
been studied by a number of investigators with 
generally good agreement.?*!!*!? These data have 
been shown!® to be consistent with a standard heat of 
formation for UN of —70 kcal/mole, in agreement with 
calorimetric determinations.'4~! © UN is found to melt 
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congruently at 2850°C under 2.5 atm nitrogen. At 
temperatures below this, and down to the melting 
point of uranium, molten uranium is the equilibrium 
product of decomposition. Under dynamic conditions, 
where diffusion of nitrogen in the UN becomes rate 
controlling, apparent volatilization of UN can 
occur.2’!7 Two studies of phase stability in the 
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and in PuN (Ref. 26), and the conclusion of complete 
interchangeability of carbon and nitrogen in the (U,Pu) 
(C,N) system appears warranted. Oxygen, on the other 
hand, exhibits only limited solubility in UN (Ref. 27). 
However, complete solution between UO, and UN, 
and of UO, with U,N3 at temperatures above 1200 to 
1300°C is reported.27-?9 


Table 1 Vapor Pressures** over UN and Up Pug ..N 





Pressure of indicated 
species over UN, atm 


Pressure of indicated species 
over U, , Pu, ..N,atm 











Temperature, 

°K Nitrogen Uranium Nitrogen Plutonium Uranium 
1400 19°'° 10°** | | ie 2x07? LO 
1600 Sx 107° Sxi0°** 25 %10-- 15x10°° 3x 10°"? 
1800 13 x0 3:0:x10” 8.8 x 10° eo 6:3 x10*° 
2000 61:x10"" 75x10 1a x20 6.7 x 10°° 4.0 x 10° 
2200 1:3:x10-° 1:0: x10° 2:0.x.10°° 5:8)x,10"° 1:0:%10-" 
2400 23:%10 * Li x40" 5.4x 10° 3.1 x 10 8.1 x 1076 





U,N3—UN, region have been reported.'*"!? Agree- 
ment between the studies is generally poor, apparently 
related to the differences in experimental techniques 
employed. 

The UN phase apparently exists at essentially the 
stoichiometric composition at room temperature. At 
temperatures of 1500°C and above, solubility ranges, 
including N : U ratios, from 0.96 and 0.92 up to 1.04 
have been reported.?*!7 


PuN is soluble at all compositions in UN (Ref. 20). 
The congruent melting point of PuN is estimated?’ as 
2830 + 50°C at a nitrogen pressure of 50+ 20 atm. 
Thus, based on an essential ideality of solution 
between UN and PuN, the congruent melting point of 
(U,Pu)N alloys is expected to exhibit little variation 
with U : Pu ratio. The mononitride is the only phase in 
the plutonium—nitrogen system, and negligible solu- 
tion of plutonium in U,N3 is indicated.?? 

In contrast to UN, (U,Pu)N does not decompose at 
elevated temperatures, but instead volatilizes such that 
the metal-to-nitrogen ratio in the solid material remains 
constant.23*?4 Volatilization is not congruent in a 
Uo.gPuo.2N alloy, since the plutonium vaporizes at a 
more rapid rate than the uranium. The result is that, 
under dynamic conditions, the rate of volatilization is 
controlled by the diffusion of plutonium in the mixed 
nitride. The vapor pressures of individual species over 
UN and Ug 3 Pug ..N are summarized in Table 1. 

The major contaminants of nitride fuels are oxygen 
and carbon. UC is completely soluble in UN (Ref. 25) 
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Heat Capacity and Thermal Conductivity 


Heat-capacity measurements of UN have been 
made by a number of investigators, and several 
reviews! *!7°3° of the resulting data have been made. 
The most recent review°°® provides the data listed in 
Table 2. Values measured”? for Ug.gPug..N also are 
given in Table 2. 

Thermal-conductivity data for UN have been re- 
viewed,?°’3! with the data given in Table 2 represent- 
ing preferred values.?® As a matter of interest, the 
conductivity of UN is less than that of UC below 
1000°K but greater at higher temperatures.?° Also 
included are data?* for Up.gPu9.2.N corrected to 
theoretical density by the equation?! 


AT.D. ” Ameas (; “) 





where 6 is 0.5 and P is the fuel porosity. The 
conductivity of UN is decreased by the addition of 
plutonium. 


Thermal Expansivity 


Linear thermal-expansion coefficients for UN 
(Ref, 32) and Up.gPuo..N (Ref. 22) as a function of 
temperature are presented in Table 3. Slight increases 
in expansivity accompany the plutonium addition. The 
expansivity of nitride fuels is approximately one-half 
of that of stainless steel, the preferred cladding 
material for LMFBR application. 
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Table 2 Heat Capacity and Thermal Conductivity of UN and Up ,Puy .N 
UN (Ref. 30) U, ., Pu, .. N (Ref. 23) 
Temperature, Cp, A, Cp, Ay 
a cal/(mole)(° K) W/(cm)(° K) cal/(mole)(° K) W/(cm)( K) 
10 1.33 0.03 
50 3.76 0.073 
100 6.17 0.095 
200 9.80 0.123 
300 11.4 0.139 11.61 
400 42.3 0.149 11.88 0.135 
600 13.1 0.160 12.40 0.149 
800 13.7 0.189 12.92 0.163 
1000 14.2 0.206 13.44 0.177 
1200 14.7 0.226 13.96 0.186 
1400 ES:2 0.245 14.48 0.195 
1600 15.8 0.250 15.00 0.204 
1800 15.52 
Diffusivity = l T l 
The coefficient of self-diffusion for nitrogen in UN O, Ref. 33 
(Refs. 33 and 34) as a function of reciprocal tempera- a x, Ref. 44 a 
ture is plotted in Fig. 1. The data are presented for a 4 
nitrogen pressure of 100 torrs, the diffusivity exhibit- 8 sol _ = 
ing a dependence on nitrogen pressure to the }, power. oo “ 
The apparent activation energy indicated for constant he on 
nitrogen pressure is 62 kcal/mole. Other data have been | ae ve = 
reported giving diffusion coefficients about an order of 2 
magnitude greater? > than those shown here. 
-12 — 4 
Table 3 Mean Linear Thermal-Expansion Coefficients | l | l l 
-13 
for UN and Uo .gPuo.2N 40 44 48 5.2 5.6 6.0 64 





Thermal expansion, 10°°/°C 








Temperature UN UlerusN 
range, °C (Ref. 32) (Ref. 23) 
20-100 7.4 43 
20—200 7.6 8.1 
20-300 8.0 8.2 
20—400 8.2 8.4 
20—500 8.4 8.6 
20-600 8.6 8.8 
20-—700 8.8 9.0 
20—800 9.0 9:2 
20—900 F2 9.4 
20—1000 9:3 a 
20—1100 9.4 2 U 
20-1200 35 93 
20-1300 9.6 10.1 
20-1400 yl 10.2 
20—1500 9.8 10.3 





RECIPROCAL TEMPERATURE, 10°/°x 


Fig. 1 Self-diffusion of nitrogen in UN (Ref. 33). Data ad- 
justed to a pressure of 100 torrs assuming Dap’. 


The self-diffusivity of nitrogen in Ug.gPug..N has 
been found to be the same as in UN, within the limits 
of experimental error.?3 


The self-diffusivity of uranium in UN has been 
determined for low and high nitrogen activities over 
the temperature ranges 1100 to 1600°C and 1420 to 
1830°C, respectively.2® For low nitrogen activity the 
diffusivity was described by D=3.24x 107 exp 
(—60,000/RT), whereas for high nitrogen activity the 
diffusivity could be expressed as D = 7.54 x 10? exp 
(—105,000/R7). 
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MECHANICAL PROPERTIES 


Elastic Properties 


The modulus of elasticity for UN has been mea- 
sured by anumber of investigators,??°’?7*! with values 
ranging from 32 x 10° to 38.4 x 10° psi being deter- 
mined. The modulus decreases*°**! essentially linearly 
with temperature to about 0.9 of its room-temperature 
value at 1200°C. 

Values of Poisson’s ratio? 7°3°*4°'4? for UN have 
been given in the range 0.26 to 0.28. Shear and bulk 
modulus values ranging from 15.1 x 10® to 16 x 10° 
psi and from 26.7 x 10® to 29.9 x 10° psi have been 
reported,?7*39»40 

The elastic modulus for Ug.g5Pug.35N has been 
measured as 29.9 x 10° psi, with Poisson’s ratio being 
0.275 for the same material.*? 


Hot Hardness 


Hot-hardness data for UN are summarized in Fig. 2. 
Impurities are believed responsible for the higher 
hardness of some specimens at elevated tempera- 
tures.*4 Compared with UC, UN is softer at tempera- 
tures below 1000°C but retains hardness to a greater 
degree at higher temperatures. 

No data are available for mixed nitrides. 


Creep 


46:47 and tensile*® creep measure- 


Compressive 
ments have been made with the results shown*® in 
Fig. 3. The compressive creep data, in particular, show 
good agreement when effects of oxygen content and 
test atmosphere are taken into consideration. 

A pronounced dependence of creep resistance on 
oxygen content of the UN is indicated by the data 
obtained*? at about 1100°C under a pressure of 
6000 psi. The increase in oxygen content from 400 to 
2600 ppM is accompanied by an increase in creep rate 
of over one order of magnitude. 

Comparison of the activation energies for creep 
obtained from Fig. 3 with those for diffusion in UN 
suggests that the creep process involves uranium 
volume diffusion. Creep data obtained in the 1100 to 
1350°C range showed activation energies of 62.0 
kcal/mole at a stress of 6000 psi and 60.8 kcal/mole at 
a stress of 8000 psi. These tests were performed in a 
vacuum of 10° torr. The higher temperature creep 
tests, performed under a nitrogen pressure of 200 torrs, 
provided activation energies averaging 75 + 10 kcal/ 
mole. The increase in activation energy with nitrogen 
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pressure is consistent witn uranium self-diffusion data 
where an increase from 60 to 105 kcal/mole was deter- 
mined as the nitrogen pressure increased. 

In-pile creep tests have been performed from which 
the degree of irradiation-induced creep has been 
calculated.5® The results are shown in Fig. 4, which 
includes for comparison similar data obtained for UO,. 
It is seen that irradiation-induced creep in UN is an 
order of magnitude below that in UO,. The reduced 
irradiation effect in UN is believed related to the more 
rapid quenching of thermal spikes or rods in the higher 
conductivity UN as compared with UQ,. 


FUEL FABRICATION 


Fuel Synthesis 


The most widely used method for laboratory 
preparation of both UN (Refs. 20, 51, and 52) and 
(U,Pu)N (Refs. 53 and 54) powders is the hydride— 
nitride—vacuum degas cycle. The function of the 
hydride step is to provide a high-surface-area powder 
product that will react more rapidly with nitrogen than 
the starting material, The vacuum degas serves to 
reduce the higher nitride, U,N3, formed in the nitride 
step to the mononitride. The hydriding is accomplished 
at temperatures in the vicinity of 230 to 250°C, 
whereas nitriding is performed at temperatures be- 
tween 400 and 850°C. Reduction to the mononitride is 
achieved under vacuum at temperatures in excess of 
800°C and below 1500°C. 

UN powder has also been produced by the direct 
nitriding of uranium metal followed by reduction to 
the mononitride.3?°5* The nitriding is apparently best 
achieved at 850°C, either at reduced pressure or under 
1 atm nitrogen. Reduction of the mononitride has been 
accomplished at 1150 to 1400°C. 

Both processes are capable of producing nitride 
with relatively low levels (200 to 500 ppM) of oxygen 
contamination. However, depending on the subsequent 
method of pellet fabrication, high density can be 
achieved with the relatively coarse, directly nitrided 
powder only by hot pressing, whereas the finer powder 
from the hydride—nitride cycle is suitable for the cold 
pressing and sintering of pellets. Cast alloys of uranium 
containing 20 wt.% or more of plutonium also are not 
amenable to the direct nitride cycle but require the 
intermediate hydride step. 

Massive UN forms of the highest purity obtainable, 
containing less than 100 ppM oxygen, can be routinely 
fabricated by arc melting.*°*5® The UN can be cast as 
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either buctons or ingots by tungsten-electrode or 
consumable-electrode techniques. Producing essentially 
stoichiometric UN requires a nitrogen overpressure of 
20 atm or greater, while melting is performed at a 
current of 130 to 150 A and an arc potential of 100 V. 
Attempts to produce mixed nitride by arc melting were 
discontinued because of the high pressure required and 
difficulty in producing stoichiometric mononitride.>7 
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The process considered to exhibit the greatest 
commercial potential at the present time is the 
carbothermic conversion of oxide by reaction with 
carbon in a flowing nitrogen atmosphere. The reaction 
is carried out in a fluidized-bed reactor at temperatures 
of about 1500 to 1700°C. The process has been 
applied to both UN (Refs. 58 to 60) and (U,Pu)N 
(Refs. 54 and 61). The relative levels of oxygen and 
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Fig. 2. Hot hardness of various samples of UN corrected to theoretical density.** 
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TEMPERATURE, °C 
1800 1700 1600 1500 1400 1300 1200 1100 

















T T = 
ey tt) YT oe, &. Io 
— _Psi_ kcal/mole Ref. _ 
*o, 3000 80.6 
z= *a , 4000 70.3 46-4 
a *o,5000 742) _| 
-_— +8, 6000 62.0 + 
- +4, 8000 60.8 47-4 
= *@, 6000 62.0 = 
x, As shown —— 48 
twas * 400 ppM oxygen 
210° +, 2600 ppM oxygen 4 
¢ im me 
= ae 
@ 
x 
a =f 1700 psi 
10-3} — 
oO — 
10°“ — 
x 2500 psi 
Y ee ee oo 
0.45 0.50 0.55 0.60 0.65 0.70 0.75 
1000/T, °K 


Fig. 3 Steady-state creep rate as a function of temperature for 
uranium mononitride tested in compression.*° 


carbon retained in the resultant product are functions 
of the partial pressures of nitrogen and of the carbon 
monoxide generated in the reaction. It appears that, 
typically, the lowest residual oxygen and carbon 
contents achievable are®®*®! in the range 1000 to 
2000 ppM. 

The use of chemical vapor deposition techniques, 
as a means of producing UN whereby uranium hexa- 
fluoride is reduced by lithium in an atmosphere 
containing a partial pressure of nitrogen, has also been 
examined for feasibility.6? However, work was insuf- 
ficient to evaluate the potential of this process. 


Perhaps one of the more interesting synthesis 
procedures is applicable to the reprocessing of nuclear 
fuels and yields nitride as the direct product in a 
pyrometallurgical separation process.°* Although 
applied to date only to UN, the process appears to be 
applicable also to mixed (U,Pu)N. The process involves 
the dissolution of the fuel in liquid tin followed by 
pressurization with nitrogen such that pure UN pre- 
cipitates from the uranium—tin solution. Because of its 
higher density, the UN sinks to the bottom where it is 
subsequently recovered by mechanical separation from 
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the melt. The UN initially produced can be processed 
through subsequent stages by redissolution in tin under 
vacuum to decompose the UN followed by reprecipita- 
tion to achieve any required decontamination factor. It 
was estimated that a decontamination factor of 10° or 
higher would be achieved in six stages. 


Consolidation 


As-synthesized nitride powders are quite resistant 
to sintering densification and must be comminuted to 
an average size of <2y to obtain sinter densities of 
95% of theoretical. Uniaxial and isostatic hot pressing 
can be used to obtain similar densities with coarser 
starting powders. As noted earlier, massive nitride 
forms can be produced directly by arc melting. 

Uniaxial and isostatic hot-pressing techniques have 
been applied to the fabrication of both UN (Refs. 32, 
51, 55, and 64) and (U,Pu)N (Refs. 57 and 65) fuel 
bodies. Uniaxial hot pressing is accomplished in a 
graphite die that can be lined with a tantalum foil to 
act as a barrier to carbon diffusion into the specimen. 
Uniaxial hot pressing has been performed under both 
vacuum®* and nitrogen atmospheres.°° Degassing of 
the nitride powder to reduce U,N3 before uniaxial hot 
pressing is not necessary, since the reduction will occur 
during the hot-pressing cycle. 

In the case of isostatic hot pressing?” the U.N; 
must be reduced prior to pressing. The nitride is 
normally wrapped in tantalum foil and is then inserted 
in graphite-lined niobium or tantalum cans for the 
isostatic, hot-gas-pressure compaction. 

Cold-pressing and sintering techniques have been 
studied and applied to the fabrication of UN (Refs. 20, 
53, 58, 66, and 67) and (U,Pu)N (Refs. 53, 54, 68, and 
69) pellets at a number of sites. The sintering of UN is 
extremely sensitive to particle size and sintering atmo- 
sphere as well as temperature, as illustrated by Figs. 5 
and 6. Although vacuum sintering is more effective 
than either argon or nitrogen sintering, decomposition 
of UN limits®? the maximum temperature to below 
1600°C. 

Attempts to develop sintering aids for nitride fuels 
have proved unsuccessful.?°*7 Reports concerning the 
effect of oxygen have been inconsistent,°*’7° and it 
may be that only when UO, is present as a second 
phase will sinterability increase, whereas small amounts 
of oxygen decrease sinterability. 

Powder prepared from  prealloyed uranium— 
plutonium metal synthesized via the hydride—nitride 
cycle has proved to be much more amenable to 
sintering than mechanically mixed powders of UN and 
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Fig. 4 Ranges of measured in-pile creep rates of UU, and UN samples standardized to a fuel density 
of 96% of theoretical, a fission rate of 2.5 x 10'* fissions/(cm*)(sec), and a stress of 2 kgf/mm’. 
Ranges of estimated rates based on theory are also shown, as are out-of-pile thermal creep rates.°* ° 
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Fig. 5 UN densification as a function of particle size after sintering 3 hr in nitrogen.° ’ 
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Fig.6 Effect®’ of atmosphere on the sintering of UN 
powders ball milled for 32 hr. 


PuN (Refs. 53, 54, and 68); in the first case, densities 
in the range 96 to 97% of theoretical can be obtained 
by sintering at 1900°C, whereas the density of 
mechanically mixed powder appears to be limited to a 
maximum of 90% of theoretical. This limitation 
appears to result from the Kirkendall effect associated 
with unequal rates of interdiffusion of the metal atoms 
in UN and PuN such that voids form at the interface 
between particles and act as a barrier to further 
diffusion.”! Sintering of prealloyed powder is** re- 
portedly enhanced by reducing nitrogen pressure at 
temperatures up to 2150°C. 


COMPATIBILITY 


Water 


The hydrolysis of UN proceeds to completion by 
the reaction: 


UN + 2H,0 > UO, + NH; + ¥,H, 


However, this reaction is seldom complete, and U,N3 
is generally observed along with the oxide.7? Further- 
more, the U,N3 forms as an intermediate layer 
between UO, and UN. It has been observed7? in 
300°C water that as U.N; hydrolyzes, the structure 
changes from the bcc Mn, 03-type structure to the fcc 
UN, fluorite structure, although nitrogen is lost in the 
process so that substitution of oxygen for nitrogen is 
suggested. This indicates that the hydrolysis proceeds 
through a continuous series of solid solutions from 
U,N3 to UN, to UO,. 
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The rate of corrosion of UN has been found to be 
extremely dependent on stoichiometry,’? the presence 
of either free uranium or U,N3 being detrimental to 
corrosion resistance. However, regardless of material 
quality, corrosion penetration rates exceeding 0.01 
mm/hr are reached7* at temperatures in excess of 
250°C. In 100°C water the film that forms on UN 
affords protection that is lost at higher temperatures as 
its thickness increases and cracking occurs. The rela- 
tively good corrosion resistance of UN as compared 
with UC is ascribed??’7* to the coherency of a duplex 
film of bec U,N3—U(N,O), superimposed on UO). 
The fcc array of atoms common to UN, U,N3, and 
UO, leads to epitaxial orientation and coherency. 


Oxidizing Atmospheres 


The oxidation of UN proceeds as in the hydrolysis 
reaction, with the formation of U,N3 as an inter- 
mediate product in an adherent, epitaxially related 
oxidation film. However, unlike the hydrolysis reaction 
where the final oxidation product is UO,, U30x, 
and/or UO; (Refs. 73 and 75 to 77) may appear as the 
final product, although decreasing oxygen pressures 
favor the formation’* of UO}. 

The ignition temperature of UN in dry oxygen and 
air has been found to vary from room temperature for 
finely divided powders to about 300°C for sintered 
products and is only slightly sensitive to oxygen 
pressure.” ° 

The oxidation of UN single crystals has been 
studied in dry oxygen, CO,, and wet air.7% The rate of 
oxidation in wet air was more rapid than in the 
oxygen, indicating that the rate of reaction in air is 
controlled by its water content. 

Of significance in the storage of these fuels, 
although UC undergoes severe corrosion due to mois- 
ture in air at room temperature, massive UN has 
shown’* no significant attack after being exposed to 
laboratory air for times up to 3 years. Although 
(U,Pu)N alloys appear less corrosion resistant, the 
resistance apparently decreasing with plutonium con- 
tent, sintered specimens of Ug.gPug..N have been 
stored in typical tank-nitrogen glove-box atmospheres 
for several weeks with only darkening of the specimen 
surfaces.” ® 


Liquid Metals 


No indication of reaction between either UN or 
(U,Pu)N and either NaK or sodium has been evidenced 
in tests conducted at temperatures ranging from 450 to 
as high as possibly 1600°C (Refs. 32, 73, and 79 to 
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81). This is true of observations from irradiation 
experiments®7’*? as well as from laboratory tests. 
However, oxide in the sodium, or oxygen in capsule 
test components which is transported through the 
sodium, has been observed’? to react with the fuel to 
form an oxide phase plus U,N3. 

Sodium wetting, as indicated by contact-angle 
measurements in sessile-drop experiments, is found to 
increase markedly at temperatures above 300°C for 
UN, Up.gPuo.2N, and 304 stainless steel.2® The 
contact angle in all cases approaches 20° or less at a 
temperature of 550°C, ensuring satisfactory wetting at 
that temperature or higher in sodium bonds involving 
these components. 


Cladding Materials 


Only limited data for nitride fuel—refractory metal 
systems exist, but these indicate excellent compati- 
bility with UN at and above 1000°C for extended 
periods of time for all but titanium and zir- 
conium.®**®5 Tungsten exhibits®® excellent com- 
patibility with UN for times to 1000 hr at 1700°C and 
for short periods of time to temperatures®” as high as 
2800°C. Molybdenum shows®**®? no reaction with 
UN after at least 50 hr at 1200°C, whereas tantalum is 
used as a diffusion barrier at 1500°C in the hot 
pressing of UN (Ref.32) and (U,Pu)N (Ref. 65). 
Niobium does not react significantly with UN at 
1000°C for extended periods of time.** 

Reaction occurs®*’8* rapidly between UN and 
nickel at temperatures above 800°C and between UN 
and aluminum above 500°C. With iron and chromium 
the only reactions with UN observed at temperatures 
up to 1000°C have been ascribed to the presence either 
of U,N3 in the fuel or of oxide in the cladding 
material. 5 

The compatibility of 304 and 316 stainless steels 
with UN (Refs. 79, 84, 85, and 90) and with (U,Pu)N 
(Refs. 79, 91, and 92) has been studied extensively, 
with excellent compatibility indicated for stoichio- 
metric nitride at temperatures up to 1000°C in tests 
involving both helium- and sodium-bonded fuel. Car- 
bon®* and oxygen”® affect this compatibility, as will 
stoichiometry of the fuel. Free uranium in the fuel 
undergoes a eutectic melting reaction with the 
stainless-steel components, whereas chromium nitride 
appears as a grain-boundary precipitate in stainless steel 
that is in contact with fuel containing U,N3. 

The excellent out-of-pile compatibility is con- 
firmed in irradiation tests involving both 304 and 316 
stainless steels. Irradiations of UN in 316 stainless- 
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steel?? cans proceeded to burnups up to 6.3 at.% at 
cladding temperatures in the range 410 to 760°C with 
no evidence of reaction, except for several specimens 
with free uranium in the fuel which reacted with the 
cladding to form a molten eutectic structure. A series 
of irradiations involving sodium- and helium-bonded 
Up .gPug.2N clad with 304 stainless steel and proceed- 
ing to burnups of 150,000 MWd/metric ton (~18 at.%) 
at cladding temperatures in the vicinity of 700°C has 
proceeded with no evidence of reaction between fuel 
and cladding.®? The only reaction observed between 
fuel and cladding occurred when the cladding tempera- 
ture exceeded 1000°C. These observations are in 
agreement with results from other irradiation experi- 
ments in which cladding temperatures of about 560°C 
for 316 stainless-steel cladding?* containing 
Uo.ssPug.15N fuel and of 1000°C for 304 stainless- 
steel cladding?! containing Uy.gPug..N were 
estimated. 

The excellent compatibility at high fuel burnups is 
observed despite predictions, based on thermodynamic 
considerations, that U,N3 will be produced in the fuel 
as burnup proceeds.*7 As in out-of-pile compatibility 
tests, evidence of fuel—cladding reaction has been 
obtained in-pile when Up.gPug..N fuel contained 
U,N; prior to irradiation.2? However, no indication of 
U,N3 formation has been obtained in (U,Pu)N fuel at 
burnups ranging up to 150,000 MWd/metric ton. 


IRRADIATION PERFORMANCE 


Nuclear Properties 


The (n,p) and (n,a) reactions are of significance 
only with '*N. The thermal-neutron cross section for 
the (n,p) reaction is about 1.75 barns, and the cross 
section decreases as 1/y with neutron energies up to 
around 0.1 MeV. Above that energy, resonance effects 
occur.57 The (n,a) reaction, on the other hand, is®* 
characteristically resonant with a fission-spectrum- 
averaged cross section of 0.041 barn. 

An estimate of the relative amounts of helium and 
hydrogen production, as compared with the amount of 
fission gas generated in the core of a mixed-nitride- 
fueled fast reactor, has been made using group-averaged 
cross-section data for the (n,p) and (n,a) reactions.?® 
A 16% increase in gas generation was calculated, of 
which about one-third was hydrogen. On the basis of 
this calculation, a summed nitrogen (n,p) and (n,qa) 
cross section about },,, of the average plutonium 
fission cross section can be estimated for a nitride- 
fueled fast reactor core. 


REACTOR TECHNOLOGY, Vol. 15, No. 2, Summer 1972 





96 NITRIDE FUELS: PROPERTIES AND POTENTIAL 


Irradiation Behavior 


Relatively large density changes have been observed 
in low-temperature (i.c., below 250°C) irradia- 
tions.27°°® Quenched-in radiation-induced defects and 
damage associated with the nitrogen (n,p) reaction are 
suspected as being responsible for the growth observed. 

Irradiations of stainless-steel-clad UN specimens, 
where fuel centerline temperatures in the range 840 to 
1320°C were achieved in conjunction with burnups in 
the range 16 x 107° to 22 x 10?° fissions/cm*, were 
accompanied by fuel swelling rates that ranged®* from 
0.3 to 0.9% per 10?° fissions/cm?. UN irradiated at 
linear power ratings near 22 kW/ft operated at center- 
line temperatures between 840 and 1100°C and ex- 
hibited swelling rates of 0.3 to 0.7% per 107° 
fissions/em?. A power rating near 30 kW/ft, where 
centerline temperatures ranging from 1100 to 1320°C 
were calculated, was accompanied by swelling rates of 
0.7 to 0.9% per 10?° fissions/cm? . Fission-gas bubbles 
appear to contribute to the swelling at the higher 
temperatures. Almost total release of hydrogen pro- 
duced in the nitrogen (n,p) reaction occurred, and 
fission-gas release rates ranged up to 26%. 

Irradiations of UN, clad with either tungsten or 
W—26 wt.% Re, at linear heat ratings in the range 3 to 
10 kW/ft in which cladding temperatures have ex- 
ceeded 1100°C have shown a strong dependence of 
swelling and of fission-gas release on tempera- 
ture.2 9193 At 2100 to 2200°C a gas release from the 
fuel of 63.5% has been measured,?? and, at tempera- 
tures below 1400°C, releases were below 10%. Single 
crystals and large-grained cast samples of UN exhibit 
significantly reduced swelling and fission-gas-bubble 
coalescence and growth as compared with normal, 
fine-grained, pressed and sintered UN (Refs. 99 and 
101 to 103). It has been observed that fission-gas- 
induced fuel swelling can be accommodated in the 
porosity of low-density (76 to 78% of theoretical) 
fuel.1°1+1°2 The swelling rate of UN, expressed in 
terms of external dimensional change of the cladding, 
is shown'°? in Fig.7. The figure includes swelling 
rates for UC and UO, determined under similar 
conditions and indicates considerably lower rates for 
UN as compared with either UO, or UC for equivalent 
surface cladding temperatures. Taking into account 
differences in thermal conductivity, for equivalent fuel 
temperatures, UN and UQ, swelling rates will be 
essentially the same, whereas the spread between UN 
and UC remains unchanged from that shown in Fig. 7. 


Swelling rates for Ugg Pug ..N fuel, from a series of 
irradiations in the Materials Testing Reactor (MTR) 
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Fig. 7 External swelling rate of UN compared'°* with those 
for UO, and UC. All specimens were clad in W—26 wt.% Re. 


and the Engineering Test Reactor (ETR) where the fuel 
was either helium- or sodium-bonded to 304 stainless- 
steel cladding, are shown®? in Fig. 8. Except for two 
groups of data points, the indicated swelling rate for 
the mixed-nitride fuel is 0.48% per 10?° fissions/cm? 
or less in irradiations performed at linear heat ratings 
up to 40kW/ft and to burnups of 150,000 MWd/ 
metric ton. The exceptions occur with fuel that 
operated with centerline temperatures above 1260°C 
and with fuel containing U,N3 as a second phase in 
conjunction with relatively high levels of oxygen 
impurity (3000 to 4000 ppM) as compared with fuels 
containing only about 500 ppM oxygen. 

The effect of fuel centerline temperature on the 
rate of fuel swelling is illustrated in Fig.9. At 
temperatures below about 1300°C, fuel swelling results 
only from the introduction of individual fission- 
product atoms into the fuel crystal structure. However, 
in the range 1300 to 1400°C, fission gases begin to 
agglomerate, and fuel swelling associated with the 
growth of fission-gas bubbles is observed. In an analysis 
performed on a single Ug g5Pug.,5N fuel pin irra- 
diated under conditions where a fuel centerline tem- 
perature of 1900°C and a fuel burnup of 12,500 
MWd/metric ton were achieved, the onset of fission- 
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Fig. 9 Temperature dependence of (U,Pu)N swelling.®? 


the presence of high levels of oxygen impurity. 

A single fuel pin containing 81.4% dense 
Up.gPuo.2N fuel, helium-bonded to 304 stainless-steel 
cladding, has been irradiated®* in Experimental 
Breeder Reactor II (EBR-II) to a burnup of 32,100 
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MWd/ metric ton (3.6 at.%) at a peak linear heat rating 
of 24.8 kW/ft. The original 4.6-mil radial gap between 
fuel and cladding was consumed by the swelling fuel, 
and an additional 4,- to 1-mil cladding deformation 
associated with fuel—cladding mechanical interaction 
was measured, In view of the initial high fuel centerline 
temperature (1760°C) estimated before fuel swelling 
closed the fuel—cladding gap, the pin swelling is judged 
to be consistent with fuel swelling rates observed in 
thermal irradiations. Despite the initial high fuel 
temperatures, fission-gas release in the pin also was 
consistent with gas-release data obtained in thermal 
irradiations (see Fig. 10). 

Poor capsule performance in irradiation experi- 
ments involving UN—20 wt.% PuN fuel clad in 316 
stainless steel and irradiated to burnups in the range 
2500 to 20,000 MWd/metric ton prevented an accurate 
assessment of fuel performance.?' However, it was 
concluded that thermal decomposition of mixed 
nitrides was not a problem at fuel centerline tempera- 
tures up to 2000°C. This conclusion concerning sta- 
bility of the mixed nitrides is confirmed by short-term 
irradiations (~8 hr) of Up.g5Puo.45N at linear heat 
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Fig. 11 Fission-gas release as a function of burnup for 
high-density (>93.8%.of theoretical) (U,Pu)N fuel.®? 


ratings of 100 kW/ft (3300 W/cm). Maximum fuel 
centerline temperatures of 2500°C were estimated,’ °* 
but no evidence of fuel decomposition was observed. 


DISCUSSION 


Interest in UN, as noted at the beginning of this 
article, is restricted to its possible use at high surface 
temperatures in space reactor applications. The 
diametral swelling data presented in Fig. 7 provide 
excellent support for this interest. Mixed nitrides, on 
the other hand, are of potential interest as advanced 
LMFBR fuels. It is this interest which will serve as the 
focus of discussion in this section. 


Fuel-Pin Design 


If the high thermal conductivity of nitride fuels is 
to be used to advantage, two possible approaches exist. 
In one approach the fuel is helium bonded to the 
cladding, as with oxide fuel, and its melting point is 
accepted as the parameter limiting the fuel-pin power 
rating. In the other approach, sodium bonding is 
employed between fuel and cladding to reduce the 
cladding-to-fuel temperature increase, and the tempera- 
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ture at which swelling rates become excessive, nomi- 
nally the temperature at which swelling associated with 
fission-gas-bubble agglomeration occurs, is accepted as 
limiting pin performance. In either case, significantly 
higher fuel-pin ratings can be achieved than with oxide 
fuels. 

The first choice has certain inherent limitations. 
Because of the high fuel temperatures, fuel swelling 
rates will be high, and, once fuel—cladding contact is 
established, reliance must be placed on cladding re- 
sistance to restrain further swelling. Fission-induced 
creep, which lowers the fuel creep resistance and 
permits the cladding to force the fuel to deform into 
available fuel porosity, is apparently a prime con- 
tributor to the successful performance of clad oxide- 
fuel pins. However, even though fission-induced creep 
is observed in nitride fuel, the present indication, based 
on data®® for UN, is that the creep resistance is an 
order of magnitude greater than in oxide fuels. 
Consequently fuel swelling will be more difficult to 
restrain, Furthermore, relatively low fuel densities may 
be required if excessive cladding deformation is to be 
avoided, resulting in a loss in the fuel-loading advantage 
afforded by high-density nitride fuel. 

A sodium bond also results in lower fuel tempera- 
tures that minimize the possibility of cladding de- 
formation as a consequence of differences in thermal 
expansion between fuel and cladding. An increase in 
average fuel temperature about twice that of the 
cladding can be tolerated without the occurrence of 
mechanical interaction. Thus fuel-management concern 
that power changes accompanying fuel relocation will 
result in fuel—cladding interaction should be elimi- 
nated. 

The principal disadvantage associated with the 
sodium bond arises from the additional cost involved in 
the loading of a sodium-bonded fuel pin. However, this 
cost should be more than offset by the potential use of 
low-cost cladding—possible because the cladding acts 
only as a container and not as a mechanical structure 
restraining the swelling fuel. The swelling fuel simply 
displaces the sodium in the gap as it expands. A 20-mil 
diametral gap has been estimated®? as being sufficient 
with 0.275-in.-diameter fuel pellets to accommodate 
the fuel swelling associated with a minimum burnup of 
150,000 MWd/metric ton. The lack of fuel—cladding 
interaction also should permit the use of thinner 
cladding, which will have a beneficial effect on neutron 
economy. 

The use of a sodium gap may introduce an addi- 
tional problem that has already been encountered with 
carbide fuels,195*1°° that is, oval cladding deforma- 
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tion produced as the result of fuel-pellet cracking due 
to thermal stresses. Cracked pellet pieces lodge in the 
sodium-filled gap and deform the cladding nonuni- 
formly as the fuel grows with increasing burnup. This 
may not be as serious a problem with nitride fuels, 
which appear to crack less drastically than carbide 
fuels. Also, the low swelling rate of nitride fuels 
permits the use of a relatively small sodium gap, and a 
smaller gap size has been found to decrease the severity 
of fuel cracking in nitride fuels.2? Nevertheless, if fuel 
cracking does prove to be a problem, a fuel shroud, 
within the cladding to contain the fuel pieces, can be 
inserted in the fuel pin. 

The performance-boundary limits applicable to 
mixed-carbide and mixed-nitride fuels are shown!” in 
Fig. 12. The limits are based on a minimum pin inside 
diameter of 0.200 in. (considered the minimum practi- 
cal pin size fabricable), a minimum specific fuel power 
of 100 W/g (doubling times become excessive at lower 
ratings), a maximum heat flux of 2 x 10° Btu/(hr)(ft?) 
(limited by cladding thermal stress and sodium heat- 
removal pumping requirements), and a maximum linear 
pin power rating of 45 kW/ft. The linear power limit is 
based on an estimate that excessive fuel swelling can be 
avoided if fuel centerline temperatures are kept below 
1300 to 1400°C. A somewhat more conservative and 
easily defined limit is the minimum temperature at 
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Fig. 12 Performance boundaries of stainless-steel-clad sodium- 
bonded (U,Pu)C and (U,Pu)N fuel pins in a fast reactor.’ °’ 
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which fission-gas-bubble swelling begins; for 
Up .sPug.2N, this is®? about 1350°C, and, for mixed 
carbide, it appears'®5*!°® to be below 1300°C. This 
would tend to decrease the linear power boundary to 
40 kW/ft or slightly less, still well above the limit of 
19 kW/ft for oxide fuels. 

There is a serious need for irradiation data on 
(U,Pu)N fuel in a fast-flux environment. With the 
exception of a single helium-bonded pin, all data have 
been obtained in thermal reactors.2? However, it 
should be noted that the performance data reported 
for the Ug.gPup..N fuel are all based on the latter 
one-half or one-third of the irradiation at which point, 
because of fuel burnout in the outer region of the fuel, 
temperature gradients in the fuel have become very 
similar to those obtained in a fast reactor. Any 
influence of helium generation in a fast flux is 
expected to be minor. Also, as long as fuel per- 
formance is limited to prohibit fission-gas-bubble swell- 
ing in the fuel, and fuel—cladding mechanical interac- 
tions are avoided by the use of a sodium bond, it is 
difficult to see any major effect of differences between 
thermal- and fast-flux irradiations on the cladding. 
Thus the thermal irradiations appear a valid means of 
projecting mixed-nitride fuel behavior. 


Safety 


A major concern with nitride fuels has been their 
thermal stability. This concern is based on the behavior 
of UN, which decomposes in a vacuum with loss of 
nitrogen at temperatures above about 1600°C. How- 
ever, as noted earlier, the addition of PuN produces a 
considerable change in vaporization behavior such that 
PuN is lost preferentially from (U,Pu)N materials and 
decomposition is prevented, presumably until such 
time as all the PuN is vaporized. 

The essential stability of the mixed nitrides is 
confirmed by irradiation experience. Irradiations of 
Up .gPug.2N (Ref.91) and Up.g5Pug.;5N (Ref. 104) 
have proceeded at fuel centerline temperatures of 2000 
and 2500°C, respectively, with no evidence of fuel 
decomposition. 

High-energy neutron capture by '*N has been 
determined as providing a more negative sodium-void 
coefficient for a nitride-fueled reactor, compared with 
a carbide-fueled one.?® The relative effect in the two 
fuels is illustrated in Table 4. The importance of this 
effect is difficult to evaluate, but, if sodium-void 
characteristics become an important feature in core 
design, an optimized nitride core could have a bene- 
ficial effect on reactor safety and fuel-cycle costs. 
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Table 4 Sodium-Void Coefficients 
of Nitride- and Carbide-Fueled Reactors? ° 





Sodium-void coefficients, % Ak 





Sodium loss from Nitride core Carbide core 





Total reactor —1.4 —0.8 
Core only +1.0 go) pes. 





Fuel Synthesis and Fuel-Cycle Costs 


There appears little likelihood that the hydride— 
nitride cycle will provide an economic method for the 
production of mixed-nitride fuels unless new methods 
for reduction to metallic uranium and plutonium are 
developed. Carbothermic reduction coupled with 
nitridation has received sufficient study to indicate its 
potential as an economic process, although consider- 
able additional development is required. However, a 
major question concerning this process is the residual 
carbon and oxygen impurity levels and their influence 
on irradiation performance. High levels of oxygen in 
conjunction with U,N3 have been shown to have a 
detrimental effect on fuel performance that U,N3 
alone does not have.®? The effects either of high levels 
of oxygen alone or of carbon impurity have not been 
defined. 

Other techniques of fuel synthesis certainly deserve 
examination. Pyrometallurgical reprocessing, _ per- 
mitting direct conversion to nitride, appears extremely 
promising. One of the attractive features of this 
process is that it should permit the economic utiliza- 
tion of '5N in mixed-nitride fuel. Such utilization 
requires '5N recovery and reuse, which are readily 
accomplished in pyrometallurgical reprocessing. 

The interest in '5N stems from the parasitic 
neutron captures in 1*N that have a detrimental effect 
on breeding characteristics and fuel-cycle costs. Com- 
parative breeding characteristics and fuel-cycle costs 
for carbide and nitride cores, determined for the same 
reactor model?® and assuming equivalent fuel- 
performance capabilities for both fuels, are shown in 
Table 5, Fuel-cycle costs for the nitride and carbide 
cores are essentially equivalent, except for credit for 
the plutonium produced in the fuel cycle. Slightly 
increased fuel-inventory requirements are also asso- 
ciated with the nitride core. 

These relative differences could be essentially 
eliminated by the use of '5N. However, costs for its 
production (90% '*N) are estimated’®® to be about 
$1.00/g. Although this represents a negligible cost 
increment to the first core cost, the cost becomes 
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Table 5 Comparative Breeding Characteristics and 
Fuel-Cycle Costs for Carbide- and 
Nitride-Fueled Reactors? ° 








Carbide Nitride 
core core 
Breeding characteristics 
Core conversion radio 0.81 0.75 
Total breeding radio 1.39 1.33 
Doubling time, years 8.4 it 
Fuel-cycle cost, mills/kWh 
Gross cost 1.1242 1.1170 
Plutonium credit 0.3743 0.2992 
Net cost 0.7499 0.8178 





significant relative to the fuel-cycle costs unless the 
15N is recycled along with the other fuel components. 


CONCLUSIONS 


Available fuel-property and -performance data in- 
dicate that the nitrides are fuel materials of great 
potential interest, UN for space reactor and (U,Pu)N 
for LMFBR applications. The data for (U,Pu)N show 
that it is thermally stable at least to 2500°C, that 
swelling rates are extremely low and predictable as long 
as fuel temperatures that permit the agglomeration of 
fission gases are avoided, that fission-gas release per- 
centages are low (less than 10% at burnups to 150,000 
MWd/metric ton), and that burnups of 150,000 
MWd/metric ton can be achieved without degradation 
in fuel behavior. The data suggest that (U,Pu)N fuels 
can be best utilized in sodium-bonded fuel pins. 

When (U,Pu)N is compared with its principal 
competitor, (U,Pu)C, as an advanced LMFBR fuel, 
certain advantages to the use of (U,Pu)N fuel appear to 
be evident—no evidence of fuel—cladding incom- 
patibility, lower swelling rate, and potentially higher 
linear power ratings without the initiation of swelling 
resulting from fission-gas agglomeration. A detrimental 
effect on breeding characteristics and fuel-cycle costs 
results from parasitic neutron captures in '*N but may 
be overcome by the use of '*N. 
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ANS GUIDELINES FOR PREPARING CRITICAL REVIEWS 
FOR REACTOR TECHNOLOGY 


One interpretation of the words “Critical Review’ emphasizes the word critica/. That type of 
article would have as its purpose the discussion of a single subject, which has been uncertain and 
perhaps controversial, in considerable depth. An example might be a critical review of the values of 
alpha for plutonium, assessing all the work done and arriving at a best current estimate. This paper 
is primarily addressed to specialists on the particular subject, and serves as an authoritative source 
for the information they use. 

In the second interpretation of a Critical Review, the emphasis is on the word review, in the 
sense of survey. Such a paper would be broader and probably more descriptive in scope. A paper 
on “Solubility of Metallic Elements in Liquid Sodium’’ is an example of this category. Such a 
paper would be addressed to a much broader group of readers and written in a fashion that would 
be of interest to the majority of reactor technologists. It would provide enough of an introduction 
that the specialist from another field could immediately appreciate why the subject is important. 

The criteria for these two types of articles may not be very different: 

1. The paper should be based on a thorough coverage of relevant work on the subject from 
many sources. A review based on the work of only one individual or group is better suited for 
publication in one of the regular society journals. The review should be critically selective, 
reporting only the most valid results and indicating why some prior results have a questionable 
status, The review should call attention to significant gaps where more work is required in the 
subject of the article. 

2. The paper should be timely, on a subject of active current interest. 

3. The scope of articles acceptable as Critical Reviews includes all the subject areas identified 
for Reactor Technology (Economics, Physics, Mechanics, Construction, Fuel Elements, Fuel 
Cycles, Fluid and Thermal Technology, Fuel Processing, Components, Operating Performance), as 
well as Materials (including Source and Special Nuclear), Environmental Effects, and Effluent 
Management. Not desired are (1) reports of original research proposed for first publication and 
(2) review articles directed toward the specialist in the field of nuclear safety (which is covered by 
the AEC’s bimonthly review Nuclear Safety). 

4. The paper should be organized so that it is of immediate practical use to the readers. Such 
Organization requires: attention-to consistent use of units, presentation of important data in 
summary curves or tables, and a fully adequate bibliography to the original literature. Details on 
size, honorarium, style, etc., can be obtained from ANS. 
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Discrete-Ordinates Methods for the 
Numerical Solution 
of the Transport Equation 


By K. D. Lathrop* 


Abstract: The theory, advantages, disadvantages, and capabili- 
ties of the discrete-ordinates method for numerical solution of 
the Boltzmann transport equation are reviewed. To illustrate 
the simplicity of the method, the author derives the theory 
directly for the specific case of one-dimensional spherical 
geometry and thereby provides a clearer picture of the method 
than is available in more general expositions. On this founda- 
tion the essential generalizations of the techniques are made. 
The major limitations of the method, including orthogonal 
geometries, negative fluxes, ray effects, convergence require- 
ments, and data storage are discussed in detail, and current 
work to remedy these limitations is reviewed. 

The capabilities of the discrete-ordinates method are 
surveyed at length, and numerical programs applying the 
discrete-ordinates techniques are tabulated and compared. An 
ancestral chart of current discrete-ordinates programs is pro- 
vided. 

A general guide to the usage of discrete-ordinates programs 
is given, and the application of the method is illustrated in 
examples including a shielding bench-mark problem, a reactor 
core-design problem, a control-rod-worth and radiation heating 
study, a pulsed-neutron problem, and a shielding parameter 
survey. 


You have a need to know. As a reactor physicist, as a 
shielding designer, as a weapons physicist, as an 
astrophysicist, or as a manager in these or related areas, 
you need to know the number of neutrons or photons 
within some assembly or in some volume. The form in 
which you need this information may vary from the 
power distribution in a reactor to the number of 
particles emerging from a shield to a stellar albedo, but, 


*Los Alamos Scientific Laboratory, University of Cali- 
fornia, Los Alamos, N. Mex. 87544. 


to satisfy your need, you must solve the Boltzmann 
transport equation. Perforce your solution will be 
numerical. The method of discrete ordinates, often 
synonymously but incorrectly termed the S, method, 
is the most powerful and general tool available for the 
deterministic solution of the complicated (seven in- 
dependent variables) and refractory transport equation. 

In many situations the diffusion approximation to 
the transport equation provides adequate solutions; in 
many others the geometric versatility of Monte Carlo 
methods provides the only means of solution. But, in a 
broad class of transport physics calculations, the 
method of discrete ordinates provides the only feasible 
means to answer accurately the myriad questions posed 
in the design of assemblies that produce or contain 
neutrons or photons. Programs based on the method 
are available to treat as many as six of the seven 
independent variables of the transport equation. Pro- 
grams are available to provide direct or adjoint solu- 
tions in each of three common orthogonal geometries 
in both one and two space dimensions. Both static and 
time-dependent calculations are possible, and solutions 
are possible for a variety of isotropic and anisotropic 
inhomogeneous or scattering sources, q variety of 
boundary conditions including the specification of 
incoming particle distributions, and a variety of eigen- 
value options. 

Although the large number of options available 
complicates the description of many programs that 
apply the discrete-ordinates method, the basic method 
is simple. The method is simple because it is designed 
to minimize the number of computer arithmetic 
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operations. Although several expositions of the theory 
of the method are available, they are written mathe- 
matically and technically for the specialist. The basic 
ideas of the discrete-ordinates method in simple 
examples, working from the specific to the general, are 
explained in the next section. A clear understanding of 
the basic principles of the method will enable the 
manager or user to see how easily the method can be 
adapted to his specific problem. 

To alert potential users to possible pitfalls and to 
help them understand and interpret their results, the 
limitations and defects of the discrete-ordinates 
method are described in the following section, and in a 
subsequent section the capabilities of the method as 
applied in currently available computer programs are 
reviewed. Both the virtues and defects of these 
programs are assessed so an intelligent choice of 
available tools may be made. 


A general guide to the usage of discrete-ordinates 
programs is provided in a later section. Partly philoso- 
phy and partly suggestions for reducing cost, this 
section contains much already known to users and 
specialists but should be of interest to those newly 
desiring to use discrete-ordinates programs. Examples 
of a broad range of applications are also given later, 
and the review is summarized in the final section. 


SIMPLIFIED THEORY OF DISCRETE- 
ORDINATES METHODS 


Presently used discrete-ordinates programs are the 
evolutionary descendants of S, programs’ *? in which 
the angular variable was represented by n connected 
straight-line segments, giving rise to the appellation 
S(egment),,. The label S, is often used as a synonym 
for discrete ordinates, although the angular representa- 
tion is different. In this section the theory of presently 
used methods is explained without reference to pre- 
vious expositions. 


The analytic transport equation expresses the rate 
of change of the number of particles in a system as the 
difference between the production rate and the loss 
rate for a point in space at a particular velocity. Thus, 
in addition to time, the equation involves three space 
variables and three velocity variables. Instead of veloc- 
ity variables, the energy and two angular variables are 
frequently used. These three variables plus the three 
space variables define a six-dimensional phase space in 
an infinitesimal volume of which, in a nuclear reactor, 
neutrons may be produced by fission or scattering 
from other elements of phase space. From the same 
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volume, neutrons may be lost by streaming, absorp- 
tion, or scattering to other volume elements. For most 
reactor physics problems, it is sufficient to deal with 
the linear transport equation, a careful derivation of 
which is given in Ref. 3, along with an explanation of 
the assumptions made in the derivation. 

The transport equation is complicated, and almost 
all attempts at its solution begin by simplifying the 
equation. The greatest simplification is to remove 
independent variables. The time variable is eliminated 
by considering steady-state situations when production 
balances losses. Spatial variables and some angular 
direction variables are eliminated by considering simple 
geometries with convenient symmetry properties. The 
treatment of energy variables is often simplified by 
making use of the multigroup approximation, in which 
the transport equation is averaged over energy ranges 
and the detailed variation of cross sections is replaced 
by suitable averages. Let us begin our description of 
discrete-ordinates methods by considering a particu- 
larly simple form of the analytic transport equation. 

In a sphere in which material properties vary only 
as a function of distance from the origin, the steady- 
state transport equation involves only two independent 
variables: the distance, r, from the origin of the sphere, 
and the cosine, u, of the angle between the direction of 
particle motion and the direction of the vector from 
the origin to the point r. In this geometry the angular 
flux W(ry.£), that is, the number of particles per unit 
volume, angle, and energy times the particle speed, is 
given by the solution to the following equation’ : 





wary), 1 a[U —u’ )W] 
ror “— Ou + ow 
Os 





= 20) fyi’ du! + Sir) (1) 


In this equation, 0; and o, represent the macroscopic 
total and scattering cross sections, respectively. The 
lower case o is used by conveniion to avoid confusion 
with the capital sigma that occurs repeatedly in sums in 
numerical work. The scattering is assumed to be 
isotropic (the probability of a scattering event produc- 
ing a particle in a direction uw after a collision when 
traveling in the direction y’ is uniformly ',), but only 
to simplify the discussion; the source of particles, S, is 
allowed to be anisotropic, i.e., to depend on yw. The 
energy variable has been suppressed. All discrete- 
ordinates programs make use of the multigroup 
approximation, although a pointwise energy repre- 
sentation is also feasible. In either case, the contribu- 
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tion of scattering events or other sources that produce 
particles in the energy range of interest is included in 
the source, S. This is one of the major simplifications 
of present discrete-ordinates methods. Each energy 
range or group can be dealt with separately. In effect, 
the energy variable is eliminated. Although a solution 
for each group may be obtained separately, all the 
fluxes needed for calculating the source to the group 
are not usually known, and an iterative cycle through 
the groups is usually performed. 

The discrete-ordinates approximation to Eq | 
evaluates the transport equation at a set of discrete 
directions Um, m=1,2,...MM. The value of the 
angular flux, or ordinate, in these directions is denoted 
by Wm(r) = v(r,um). Because Eq. 1 contains both an 
angular derivative and an angular integral, the range of 
is divided into intervals Au, = Um+% — Um —%, where 
—1 =yy<uy...<uyy+y=tl. This enables repre- 
sentation of the derivative as a difference and repre- 
sentation of the integral as a sum. Thus the discrete- 
ordinates approximation to Eq. | is written 


Km a(r? Vm) , ( —Bin+%)Vm+ — (1 — Bm —%)Wm—% 
r or rAum 





MM 


FE ve Aue +Sm(r) (2) 


+ oWim = 
This is neither the form of the discrete-ordinates 
equation as used by Chandrasekhar® nor that used in 
present discrete-ordinates codes; both are more general. 
For the moment, a description of this generality is 
postponed. 

Equation 2 does preserve one important property 
of both the currently used methods and the analytic 
equation. If Eq. 2 is multiplied by Au,, and summed 
over all directions, the term representing the angular 
derivative vanishes because ui, = ufgm+y = 1. The same 
is true of the u integral of Eq. 1. The angular derivative 
accounts for the change of direction of particles as 
they stream in straight lines between collisions. The net 
number of particles per unit volume and energy is not 
affected by this redistribution. If the numerical 
approximation to Eq. 1 did not preserve this property, 
local conservation of particles would be lost. 

Equation 2 remains a continuous function of r. Let 
us now discretize this variable and evaluate the 
equation at the points r,, r2,...;, in the intervals 
Ar; = ri+% — Vi-y with 0 al ry< ore <Try+%: The 
value of the discrete-ordinate angular flux at these 
points 7; is denoted by W;,,, and Eq. 2 is approximated 
by 
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Um (ri+yVity,m ry Vi-y,m) 
r? Ar; 





F a —Min+%)Vi,m+%—( —in ~% )Vi,m—% 
rj Aum 
MM 
Evie Aue + Sim (3) 





a 
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Q 


+ O4;Vim = 


The domain of the variables u and r has now been 
divided into MM x IM boxes of size Au,, x Ar;, and 
Eq. 3 involves y evaluated inside and on the edges of 
these mesh boxes. Figure | shows how these fluxes are 
oriented for a typical mesh box. The location of r; 
within Ar; and the location of u,, within Au,, does 
not need to be specified unless the accuracy with 
which Eq.3 approximates Eq.1 is going to be ana- 
lyzed. 

For each box the five angular fluxes shown in 
Fig. 1 are coupled by Eq. 3. Two of these five fluxes 
are known, either from boundary conditions or cal- 
culations in adjoining boxes. For example, for u,, <0 
(incoming directions), the commonly used “vacuum” 
boundary condition means that Yyy+y m =0. Two 
other equations are required for use witn Eq.3 to 
allow calculation of the three unknown fluxes. The 
simple, symmetric, and accurate arithmetic mean rela- 
tions 


Vim = (Wi+y,m + Vi-y,m)/2 
(4) 
Vim = (Vim+y%t Vim—y)/2 


are frequently used to provide the two extra equations. 
If it is assumed that Yj+% m and Wj; »—y are known, 
then the use of Eq. 4 in Eq. 3 permits the center flux 
to be calculated in terms of these two known fluxes: 


Vim = GimVity,m + bimVi,m—% + Cim Qin (5) 


where the form of ajm, Dim, and Cjm is easily 
determined by performing the substitution and where 
Qim includes both Sj; and the scattering source. If 
Vi-y,m is known instead of Wj+y, m, a similar re- 
cursion can be found. For this one-dimensional spheri- 
cal geometry, it is always possible to proceed from 
small values of m to large values. For the first value of 
m, Eq.5 is used either with values of ¥; 4 obtained 
from a special calculation for 4 =—1 (see Ref. 6) or 
with a difference scheme such that djg = 0. 

Simple equations like Eq. 5 form the inmost part 
of almost all numerical programs that solve the 
discrete-ordinates equations. The complicated partial 
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Fig. 1 Typical discrete-ordinates mesh in one-dimensional spherical geometry. Location of cell edge 


and center fluxes is shown. 


differential equation of Eq. 1, or its counterpart in 
other geometries, is reduced to an algebraic system that 
can be solved recursively. Every flux on the mesh can 
be computed by stepping from box to box. Besides 
being simple, equations like Eq. 5 are very general. The 
cross sections may be continuously variable from one 
space interval to the next. Although not indicated in 
Eq. 3, because most discrete-ordinates codes do not 
allow it as an option, the cross sections could also 
depend on angle. A large variety of useful boundary 
conditions are easily applied. Similar recursions can be 
developed for the adjoint transport equation so that 
most discrete-ordinate numerical programs can be used 
to solve the adjoint transport equation as well. The 
versatility of such programs will be illustrated later. 

Most discrete-ordinates programs for numerical 
solution of the one-dimensional transport equation use 
a generalization of Eq. 3, which can be written 
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V; 
Ait — Ai-%)Om+4¥i,m+ — Um —4Vi,m—#) 
Vidum 
0. MM 
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+ O1iVim => x Wig Aug + Sim (6) 


In this form, Aj+ is the area of the surface at rj+%, and 
V; is the volume of the cell of width Ar; in whichever 
geometry is being considered. For the sphere, Aj+%= 
4arz,y and V;= 4n(r3.y4— r?_y), which leads to a 
slightly different representation of r; than that used in 
Eq. 3. The @ coefficients are the analog of (1 — u?,+%) 
and are defined so that the two approximate deriva- 
tives in Eq. 6 cancel when y is a constant, as would the 
analytic derivatives in Eq. 1. This requires that 
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Omn+% — An—% = —Mm Am (7) 


Like 1 — yi, the first a, ay is required to be zero so that 
Eq. 7 serves as a recursion to determine the rest of the 
coefficients. With Eq. 7, the use of the arithmetic mean 
difference relations in Eq. 6 produces a recursion like 
Eq.5 which has the same a, b, c coefficients for 
stepping from left to right or from right to left. This is 
not true of Eq.3 unless Um» = (Um+% + Um—%)/2, a 
requirement not satisfied by many useful quadrature 
sets (Um, Aum). 

With equations like Eqs. 4, 6, and 7, one can make 
a numerical program that treats plane, cylindrical, or 
spherical geometry simply by using the appropriate 
definition of A and V, an elegant embedding. Further- 
more, when Eq. 6 is multiplied by Au,,V;, each term 
in the equation is identifiable as a number of particles 
contributing to the gain or loss in the phase space box. 
For example, the term AUmUmAj+%Vi+%,m is the 
number of particles that flow out of the box through 
the surface at rj+%, Auy,Vj04;Vim is the total number 
removed from the box by collisions, and so on. The 
approximation to the analytic equation is boxwise 
conservative, a very useful feature. 

The two-dimensional orthogonal geometries (x,y) 
planar, (0,0) cylindrical planar, and (,z) finite cylindri- 
cal can all be treated by using a generalization of Eq. 6: 





Um (Ai+%,jVit4,im—Ai_%,/Vi_-%, jm) 








Vij 
+ Aire, jAi-% )Gm+ 2Viim+4—%m—*Vijm—) 
VigWm 
Nm Bij(Vi,j+%,m—Vi,j—-%,m) 
+ 
Vi; 
Osij x 
+ OrijVij = 57 Py We Wijo + Sim (8) 


Here the subscript j indicates the second spatial 
direction, with n,», being the cosine of the angle 
between the direction of particle motion and the 
direction of the axis along which j distances are 
measured. In these geometries the angular integral is 
actually a double integral, but it can be accomplished 
as a single sum with quadrature weight W. The extra 
space derivative in Eq. 8 requires the use of another 
arithmetic mean relation, that is, 2Wijm = Wi,j+%m + 
Vi,j-%,m> but otherwise the solution of Eq. 8 is much 
like the solution of Eq. 6. In both cases, one starts by 
assuming that some system exterior-edge values of the 
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flux are known and that the tight-hand, or source, sides 
of the equations are given. Then one uses Eqs. 6 or 8 to 
calculate the box center flux and equations like Eq. 4 
to calculate the other edge fluxes. Having these fluxes, 
one can move to the next cell to the interior and repeat 
the process until all fluxes in the system are known. If 
Os is zero, this completes the calculation. If not, the 
scattering portion of the source is calculated using the 
latest fluxes. Rather than store the angular flux, one 
usually accumulates, during the sweep through the 
mesh, those integrals of the angular flux which are 
necessary to calculate the scattering source. With the 
new estimate of the scattering source, the process is 
iterated. We illustrate this procedure with a very 
simple, but two-dimensional, example. 

For the simple four-space box system shown in 
Fig. 2, we assume that 0; = 1, o, = 0, and that S = Q in 


Angular quadrature 
m WM em 1m 
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Fig. 2. Geometry for simple two-dimensional test problem. 


each box in each direction. That is, the source is 
constant in space and isotropic in direction. For this 
geometry, Vj;= AxAy, Aj+y; = Ay, and By = Ax, so 
that Eq. 8 becomes 


Um (Wi+y,j, m — Vi- %,j,m) 





Ax 


as Vij-em) y=Q (9) 
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If the arithmetic mean relations like Eq. 4 are used, 
Eq. 5 becomes 


4 2|MmlAy Vi+y,j,m + 2lNm Ax Wi j+%,m + AxAyQ 
= 2M Ay + 2In»,IAx + AxAy 


= aVi+y jm + bWi,j+y,m +cQ (10) 





where for the problem at hand, a, b, c are the same for 
each box for each direction and a+ b+ c=1. The 
quadrature uses four symmetrically oriented directions 
with cosines uy = +d and n = +X. If X = 14% , the quadra- 
ture is the commonly used S,,set, which has the fewest 
discrete directions of any S, quadrature. 

To simplify our nomenclature, for this problem, we 
consider only direction 1 and replace the i subscripts 
with a box subscript and a letter, using C, L, R, T, and 
B to denote center, left, right, top, and bottom fluxes 
in the box. Thus Ws5,;2.2,; becomes 4 Rr, or the 
box 4, right-edge flux in direction 1. In box 4, for this 
direction, Y4,7 and W4 R are needed to evaluate 
Eq. 10 for 4,c. But, with vacuum boundary condi- 
tions, as indicated in Fig. 2, both of these fluxes are 
zero, and W4 ¢ is simply cQ. From Eq. 4, W4,7 = 2cQ 
and W4 p= 2cQ. We can now move to box 3, where 
the boundary condition implies 3 7 =0 and where 
W3,.R =Wa4,~, Which was just calculated to be 2cQ. 
Thus 


W3,c = a(2cQ) + cO = (2a + 1)cO 


Having this value, we can calculate Y3 7 and Y3 p. At 
this point, we have generated all the top-row fluxes for 
direction 1. We now can start direction 2 from box 3, 
where both boundary fluxes are known, and step to the 
right. For this simple problem, these fluxes are the 
same, from left to right, as those we just calculated 
from right to left. Having traversed from left to right 
and right to left along the top row of boxes, we have 
generated all the top-edge fluxes for boxes 1 and 2. 
These two boxes can then be treated just as if they 
were at the top of the system, so we can traverse from 
right to left for direction 1 and from left to right for 
direction 2 to compute all of these fluxes. We have 
now reached the bottom of the system, where we can 
begin solution for the fluxes in directions 3 (right to 
left traverse) and 4 (left to right traverse). Then we 
move up to boxes 3 and 4 and complete the calcula- 
tion. In Fig.3, we show the magnitudes of all the 
fluxes for direction 1, in units of cQ. Because of the 
symmetry of the system, this diagram gives us every 
angular flux. In fact, the angle-integrated or scalar flux 
at each cell center is given by 
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4 
dij = z We Wig = (1 tatbta® +b*\cQ 

It is this flux that would be used to form an isotropic 
scattering source, if o, were not zero. Then the traverse 
through the mesh would have to be repeated, adding 
this source to Q. The whole process would have to be 
iterated until the amplitude of the scattering source 
converged. 
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Fig.3 Angular fluxes for the edges and centers of the 
in—down direction (units of cQ). 


With more directions the traverses are made: (1) in- 
ward, beginning at the top-right corner for all in—down 
directions, and (2) outward, beginning at the top-left 
corner for all out—down directions. This inward— 
outward sequence is repeated for each row of the 
system, beginning at the top and moving to the 
bottom. Then the inward—outward traverse is begun 
for all upward directicns in the bottom row, and these 
Sequences are repeated a row at a time until the top is 
reached. This entire cycle comprises one iteration. 
Anything that reduces the number of iterations or the 
computational work per iteration shortens the calcula- 
tion. In the above example, it would have been possible 
to reduce the number of operations by using the 
proper boundary condition. For example, if specular 
reflection had been used as boundary condition for the 
left and bottom sides of cell 4, only one-fourth as 
much effort would have been required. 
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Basically the above procedures for solving the 
elementary algebraic recursion are very simple and are 
ideal for implementation on a computer. The programs 
that implement these procedures are not simple, 
however, primarily because of the variety of problems 
each program can handle. We have considered only the 
portion of the computation known as the inner or 
within-group iterative cycle. Here we have allowed 
generality of cross sections, geometry, and boundary 
conditions. Most discrete-ordinates programs apply 
algorithms of varying complexity to accelerate the 
inner iterative process. In two-dimensional geometries, 
this accelerative algorithm may itself be iterative. Both 
this process and the inner iterations themselves may 
form part of an outer iterative cycle. For example, if 
particles are fed from one group to another by fission 
or upscattering, then the flux in the group being solved 
by the inner iteration depends on the flux in groups 
that have yet to be treated. Another iterative cycle, 
usually called outer iteration, is required. If the 
problem is homogeneous, i.e., if it has no sources 
independent of the flux, an eigenvalue can be deter- 
mined during the outer iterative process. In neutron 
problems, Keg is determined during the outer cycle. 
Most discrete-ordinates programs allow searches for 
such system parameters as zone thicknesses or material 
concentrations that guarantee a certain value of Kerr. 
These problems require sequences of outer iterations. 
A useful example of such a computation is the “time 
eigenvalue” problem, where the time-dependent flux is 
assumed to vary as exp (at), and a sequence of outer 
iterations is performed to find that value of a 
corresponding to a prescribed value of Kerr. 

Both these parameter iterative sequences and the 
outer iterative cycle can be accelerated, and many 
algorithms are used to speed up the computation. In 
addition to all this machinery, most programs have 
sophisticated input loading and checking routines, 
provisions for detailed editing and dumping of com- 
puted information, and options for punching, plotting, 
or passing on information to other programs. The bulk 
of computation time is spent executing the basic 
recursion; the bulk of program maintenance time is 
spent tending the peripheral routines. Options are 
relatively easy to build into a new program but are 
much more difficult to thread into existing programs. 

The general theory of the methods used in most 
discrete-ordinates programs is given in Ref. 6; a more 
detailed exposition of the theory for a specific two- 
dimensional program is given in Ref. 7 and for a 
one-dimensional program in Ref. 8. In the next section, 
defects of these methods are discussed. 
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LIMITATIONS OF DISCRETE- 
ORDINATES METHODS 


Geometry 


The discrete-ordinates method is simple, it is 
general, and it is most versatile. The method has 
limitations, however. Foremost of these is its present 
practical restriction to one- and two-dimensional 
orthogonal geometries. The reactor shown in Fig. 4 is a 
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Fig. 4 Real reactor configuration. 


piece of precision engineering tnat was not designed to 
fit the orthogonal grid of a discrete-ordinates calcula- 
tion. Figure 5 shows how reality must be forced into 
artificial boundaries that are the same for each energy. 
Similarly, Fig. 6 shows a planar cross section of a 
reactor and Fig. 7 shows the simplified representation 
when, where possible, advantage is taken of the system 
symmetry to describe only a fraction of the system. 
Although Figs.5 and 6 imply representation of the 
whole reactor, this is only occasionally possible, and 
then only with limited cesciiption of detail. The 
number of unknowns, and hence core storage require- 
ments and computational expense, in a discrete- 
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ordinates calculation is directly proportional to the 
product of the number of groups, mesh cells, and 
directions used. In almost any two-dimensional reactor 
description, the user of a discrete-ordinates program is 
asked to compromise between accuracy, computational 
expense, and realistic description of geometry. This 
need to compromise leads to the popularity of cell 
calculations in which a representative segment of a 
reactor core is calculated in great detail to provide 
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Fig. 5 Orthogonal grid representation of reactor in Fig. 4. 


parameters that may be used in cheaper diffusion- 
theory calculations of the entire system. The need to 
compromise is also the reason no general three- 
dimensional discrete-ordinates program has _ been 
written. Several existing two-dimensional discrete- 
ordinates programs perform more arithmetic opera- 
tions per cell than a three-dimensional (x,y,z) program 
would require, but the number of mesh cells required 
for describing real three-dimensional systems is so large 
that the computational expense would become pro- 
hibitive. 

Orthogonal grids, then, are inefficient. The coming 
decade will see the development of discrete-ordinates 
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programs with improved ability to describe geometric 
detail. Work has already begun? on a program to 
describe planar, regular hexagonal meshes such as that 
in the real system of Fig.6. A proposed general 
triangular mesh for two-dimensional calculations has 
been described by Reed.'® 


Negative Fluxes 


Designers of numerical methods dream of finite- 
difference schemes that are economically cheap, ac- 
curate, and positive. As we have seen, large systems 
(and large is measured in terms of particle mean free 
paths over the energy range of interest) strain the 
capacity of numerical programs with the requirement 
for large numbers of mesh cells. If the arithmetic mean 
relations described earlier are used to solve the finite- 
difference form of the transport equation, large 
numbers of spatial cells may also be required to ensure 
that smoothly varying positive fluxes are obtained in 
situations where the flux should be smoothly varying 
and positive. These arithmetic mean relations are 
astonishingly accurate. Not only do errors decrease 
rapidly as finer mesh intervals are used but also answers 
obtained with crude meshes are often satisfactory. The 
relations also lead to cheaply evaluated recursions, as 
we have seen. They do not guarantee a positive 
solution, however. By examining the recursion relation 
in one-dimensional plane geometry, one can easily see 
that a negative flux could be generated, in the absence 
of sources, whenever 0,Ax/2u,, > 1, that is, whenever 
the effective particle path across the cell is large 
compared to two. Usually, sources are present, or 
negative fluxes occur only in some directions so that 
the angle-integrated flux, while being distorted per- 
haps, does not become negative. Shielding applications, 
where particles stream through large thicknesses of 
source-free, absorbing material, are a notable, impor- 
tant exception. 

For most practical applications, it is feasible to 
decrease Ax (or Ar) in one-dimensional calculations 
until satisfactory fluxes are obtained. This is especially 
true when corrective recipes are used. In many two- 
dimensional problems, it is usually not feasible to 
decrease the size of mesh intervals sufficiently to 
prevent every negative flux, so considerable effort has 
gone into developing corrective recipes.'’ All these 
corrective recipes involve local modification of the 
arithmetic mean relations but can be divided into two 
categories: “‘set-offending-flux-to-zero-and-recompute” 
algorithms and ‘‘switch-to-a-positive- 
difference-scheme” algorithms. The former is more 
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Fig. 6 Actual reactor core configuration. 


accurate than the latter, but each requires continuous 
testing to detect negative fluxes. Thus, in addition to a 
basic recursion, tests must be made to detect a negative 
flux, and, if one is detected, extra arithmetic must be 
performed to calculate a positive flux. 

To avoid this proliferation of algorithms and 
consequent increase of complexity in programs, I have 
attempted to provide strictly positive difference 
schemes for the solution of discrete-ordinates equa- 
tions.'? Although these schemes have been particularly 
useful in parameter surveys and although they interact 
favorably with the iterative process to reduce compu- 
tational time, their basic inaccuracy is such that one 
usually returns to the accurate arithmetic mean 
schemes. In Table 1 the properties of commonly used 
difference approximations are summarized. The assess- 
ment of the properties is based on performance in 
two-dimensional calculations. 

The entries of Table 1 are arranged in order of 
increasingly positive solution. The more positive the 
Fig. 7 Orthogonal grid representation of core cross section in scheme, the better its convergence properties, but the 
Fig. 6. worse its accuracy. This is the kind of tradeoff that 
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Table 1 Properties of Discrete-Ordinates-Difference Approximations 





Interaction with 





Scheme Accuracy Positivity iterative cycle 
Basic arithmetic Excellent No Fair 
mean, or diamond 
scheme! ? 
Same as above, but Good to excellent Yes Fair to good 


with “‘set-to-zero- 
and-correct”’ 
algorithm 
inclided7*® 


Weighted diamond Goou to excellent 
scheme with fixed if weights can be 
weights! 4 chosen close 


enough to diamond 


weights 


Variable weighted 
diamond scheme?! 


Fair to good 


Yes, provided Good to excellent 
weights are chosen 

properly (by 

trial and 


error) 


Yes Excellent 





gives the methods designer nightmares. During the next 
decade the search for the simple, accurate, positive, 
and quickly converging scheme will continue. There is 
some evidence that most of these objectives may be 
obtained if the explicit nature of the algorithm is 
sacrificed’*; i.e., if one is willing to use a relation 
involving angular fluxes from previous iterations. 


Flux Distortions and Ray Effects 


In addition to not being positive, the simple, 
arithmetic mean scheme produces oscillatory distor- 
tions of the flux. These can be considered a subcase of 
the negative-flux syndrome. Note in Fig. 3 that, if a is 
<4, , the cell-centered flux of box 3 will be less than 
the right-edge flux for the same cell. Had we used more 
boxes in the calculation, we could have observed, for 
a <},, a spatially damped oscillation of the flux as we 
progressed inward. Roughly speaking, this oscillation, 
from cell edge to cell edge, is around the exact 
solution. Because the cell-centered flux is on a straight 
line between two edge fluxes, it can be very accurate, 
and all integral results of discrete-ordinates methods 
are based on cell-centered fluxes. Occasionally, for 
instance, when emerging angular fluxes are of interest, 
these oscillations can become important. They are 
controlled by the same steps required to make fluxes 
positive. Normally the oscillations are only an aestheti- 
cally annoying price paid for accuracy. 


A more serious spatial distortion occurs due to a 
defect in the discrete-ordinates formulation itself. 
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When we evaluated Eq. 1 at a set of discrete directions 
to get Eq. 2, we changed the nature of the equation. 
The analytic transport equation describes the streaming 
of particles in all directions; the discrete-ordinates 
equation (even with analytic space derivatives) allows 
particles to stream only in discrete directions. This is 
not important in one-dimensional geometries, and the 
seriousness of the defect was not noticed for many 
years. 

In two-dimensional geometries, for systems includ- 
ing isolated sources in absorbing media, this restricted 
streaming of particles along discrete directions or rays 
can lead to serious distortions in the flux. It is not hard 
to understand why. For the extreme case of an isolated 
source in a pure absorber, it is easy to find (outside the 
source) positions in the system not intersected by an 
allowed straight line from the source. Because the only 
contribution to the flux at such positions is that due to 
streaming, the calculated flux is zero. However, such 
positions can be found adjacent to locations that are 
intersected by allowed directions. Here the flux is 
appreciable, so the flux varies rapidly from zero to 
some significant value. Such a distortion of the true 
flux behavior is called a ray effect.'® It can be 
particularly important if the detailed space variation of 
the flux is important. Two-dimensional shielding cal- 
culations, especially, can be seriously affected by the 
ray effect. Any time the value of a flux at a point is 
required and the flux is attenuated by several orders of 
magnitude in passing from source to detector, the ray 
effect can cause serious errors. 
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There are many remedies for the ray effect.'” 
Using more directions, or more properly oriented 
directions,'*® provides more allowed paths from source 
to detector. Scattering, which provides a local reemis- 
sion source after a collision, mitigates the effect. It is 
often possible to calculate, either analytically or very 
accurately, the flux that has not suffered a collision 
and use this flux to provide an undistorted source of 
first collisions to be treated by discrete-ordinates 
methods. 

Perhaps the most powerful remedy for ray effects 
is the alteration’’ of the equations to make them 
similar to spherical harmonic equations, which are free 
from ray effects. The nature of the alteration is the 
introduction of a special source that emits particles 
selectively into regions not reached by direct streaming 
paths. 

We show examples of ray effects in the section 
titled “Applications.” One reassuring aspect of cal- 
culations performed comparing discrete-ordinates solu- 
tions to exact solutions and other solutions free from 
ray-effect distortions is that solutions that manifest ray 
effects seem to oscillate about the correct solution. 
Consequently global quantities can be estimated ac- 
curately. For example, the total number of particles 
emerging from a two-dimensional shield is likely to be 
accurately predicted, although the precise location of 
emergence is not. 

The next 10 years will see continuing efforts to 
modify discrete-ordinates equations to eliminate ray 
effects. The thrust of these efforts will be to make such 
modifications a minor, easily executed part of the 
overall calculation. 


Convergence Acceleration 


The discrete-ordinates method relies heavily on 
iterative algorithms for solving the transport equation. 
The size of the (NV x N) matrix involved is simply too 
large for direct inversion. In two-dimensional calcula- 
tions, V, the product of the number of space points 
and the number of directions, can be of the order of a 
million. 

In marked contrast to diffusion-theory methods, 
the convergence of discrete-ordinates inner iteration 
algorithms is almost independent of the size of the 
space mesh or the number of angles used. Nevertheless, 
there has been a constant effort to accelerate the 
iterative process. A wide variety of algorithms has been 
applied. The success of these algorithms has been far 
from uniform. No straightforward method has been 
found which works well for all classes of problems. 
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The most promising of all acceleration procedures 
are variations of a particle-conserving renormalization 
suggested by Bell.? Now referred to as “coarse-mesh 
rebalancing,” these algorithms determine local factors 
by which the flux is multiplied to ensure particle 
conservation. Their theory has been examined by 
several authors,’ *?! and Reed?? has recently shown 
that the method is a nonlinear variant of the synthetic 
method analyzed by Gelbard and Hageman.?* By 
examining the easier to analyze synthetic method, 
Reed has been able to develop a coarse-mesh rebalanc- 
ing procedure that works well for even the most 
difficult of one-dimensional problems. Possibly a ver- 
sion of Reed’s technique will also work well in 
two-dimensional calculations. Present applications of 
standard coarse-mesh rebalancing techniques have been 
made to work well in most two-dimensional applica- 
tions,’ but only by adhoc recipes such as periodic 
application and switching from coarse-mesh to whole- 
system rebalancing in the event of trouble. 


No matter how successful such developments may 
be, research on accelerative procedures will continue 
because of the potentially large reward in terms of 
increased problem-solving capability. A relatively un- 
explored area of development that seems fruitful is the 
application of block inversion techniques that rely on 
the direct inversion of small portions of the equations 
of the system and the iteration of the coupling 
between blocks. Such methods have been pioneered in 
discrete-ordinates computations by Wagner.?* 


Data Storage and Retrieval 


One of the big advantages of the discrete-ordinates 
method, the ability to generate simultaneously the 
entire space-, energy-, and angle-dependent distribution 
function, is also a disadvantage of the method because 
all that data must be handled, stored, and processed. 


Most discrete-ordinates programs do not store the 
angular flux but make an extra iterative cycle to store 
it if it is needed for output. Moments of the angular 
flux necessary to generate scattering sources are stored. 
With increasingly detailed calculations, more and more 
such moments must be saved. A four-term Legendre 
scattering expansion requires 10 such moments to be 
stored in two-dimensional calculations. Each of these 
moments must be stored for each space mesh for each 
group. 

All the two-dimensional programs now in use solve 
the transport equation one group at a time. Data 
required for the calculation of other groups are stored 
peripherally. More advanced two-dimensional pro- 
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grams?* solve the equations one line at a time, as if the 
equations were one dimensional. Such treatments have 
been made possibie by efficient use of rapid access 
peripheral storage units. Such developments are in- 
dicative of a general trend that will probably continue 
until the entire angular flux is stored. Moments of the 
angular flux will then be calculated as needed, rather 
than stored. With the angular flux available, new 
difference schemes and new convergence acceleration 
algorithms are possible. 

Part of the data-transfer problem will be alleviated 
as the ability to transfer data from bulk fast memory 
to disk or tape units becomes easier, while calculation 
is in progress. Most discrete-ordinates programs are 
computation bound. That is, they expend most of their 
effort in computation rather than in data transfer. 
Thus it seems feasible that much data transfer could 
take place without a significant increase in the total 
computational time. The trend toward increasing 
spatial and energy resolution in design calculations will 
speed the development of this capability. 


CAPABILITIES OF DISCRETE-ORDINATES 
METHODS AND PROGRAMS 


In previous sections, some capabilities of discrete- 
ordinates methods were indicated. In this section, all 
the available capabilities are reviewed, and easily 
realized additional potentials of the technique are 
indicated. Not all the capabilities are of immediate 
interest in any one application, such as reactor-physics 
design, but most may eventually prove useful. Also 
reviewed are numerical programs that have been 
written to apply the method and to evaluate the most 
generally available modern programs. Programs are 
categorized in several ways, beginning with their 
ancestry, and their effectiveness is assessed both in 
terms of the number of options offered and in terms of 
the ability to cope with the limitations discussed in the 
preceding section. 

Much fine discrete-ordinates work has been done at 
the naval reactor laboratories, Bettis and Knolls, and at 
the weapons laboratories, Livermore and Los Alamos. 
Some of this effort is undocumented, and some of the 
resulting numerical programs are not generally avail- 
able. In this review, no attempt is made to present 
complete coverage of these programs. 


Geometry 


Discrete-ordinates programs are available for both 
time-dependent and static problems in one-dimensional 
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plane, cylindrical, and spherical geometries. Discrete- 
ordinates programs also can be made to calculate 
plane-geometry problems in which the angular flux 
varies as a function of azimuthal angle about the x axis, 
as well as a function of the cosine of the angle between 
the particle direction and the x axis.?® Applications 
for this capability are most common in astrophysics. In 
two dimensions, programs exist for solving (x,y) 
planar, (9,0) planar, (p,z) cylindrical, and (7,8) spherical 
geometry problems. In this latter geometry, the system 
is assumed homogeneous in angle of revolution about a 
pole but can vary as a function of angle from the pole. 
Developmental work is under way for hexagonal planar 
geometry? and general triangular geometry pro- 
grams.'° In the latter work, the triangular grid can be 
applied to planar or cylindrical geometries. The only 
time-dependent, two-dimensional geometry program 
available is for (p,z) cylindrical geometry.” 7 


Boundary Conditions 


Boundary-condition options are quite general. It is 
possible to specify the angular flux at each boundary 
point for each incoming direction. Boundary condi- 
tions involving the flux being calculated (implicit 
boundary conditions that may involve iterative solu- 
tion even in the absence of scattering) are of two 
types: albedo and periodic. The vacuum boundary 
condition is a special case of an albedo boundary 
condition. The most general albedo condition would 
allow the incoming flux in each direction at each space 
point for each energy to be an arbitrary linear 
combination of the outgoing angular fluxes in the same 
space—energy interval. For the most part, such gener- 
ality is neither needed nor allowed. The specular 
reflection (angle of reflection equals angle of in- 
cidence) and white (all outgoing particles returned 
isotropically) boundary conditions are commonly used 
examples of albedo conditions. Periodic boundary 
conditions, useful in the representation of periodic 
lattices in planar geometries or in the representation of 
periodic variables in cylindrical and spherical geome- 
tries, set the incoming angular flux at one boundary 
equal, angle by angle, to the outgoing angular flux at a 
periodically located boundary. 


Source Options 


In addition to being able to specify the incoming 
angular flux at a boundary, it is also possible to 
introduce an angularly dependent source at interior 
surfaces. These are sometimes referred to as “shell” 
sources. Most discrete-ordinates programs will accept 
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an isotropic source of particles for every space interval 
and energy group of the system. Some accept aniso- 
tropic distributed sources in terms of the spherical 
harmonic expansion coefficients of the angular source. 
In principle, a detailed angular source could be 
introduced for every space—energy cell. 


A source due to first collisions can be entered in 
most discrete-ordinates programs by using the regular 
distributed source options, but several programs will 
automatically calculate such a source and use it for 
subsequent calculation of the collided flux. At present, 
the uncollided sources that are calculated are from 
point, line, or plane sources for which analytic evalua- 
tion of the exact flux is simple. Not only do such 
sources improve accuracy and mitigate ray effects but 
they also can speed convergence of the interative 
process. More important, the sources permit the use of 
lower-order angular resolution for a given accuracy. It 
is only a matter of time before more general uncollided 
flux computations are possible. 


Angular Quadrature 


Much work has been devoted to the development 
of quadrature sets for discrete-ordinates calculations, 
that is, to selecting the discrete directions in which the 
transport equation is evaluated and to selecting the 
weight associated with each direction.2* 3? Most 
numerical programs in use accept arbitrary quadrature 
sets as input, including Gauss Legendre (P,,) and 
Double P, for one-dimensional integration, but require 
that such sets meet certain requirements. Perhaps the 
most stringent of these is the triangular arrangement of 
directions on octants of the unit sphere as advocated 
by Carlson.?*-3! This arrangement is economical of 
-points and permits a two-dimensional angular integral 
to be executed as a single sum by sweeping over sphere 
latitudes along which one of the angular variables is 
constant. The arrangement also satisfies basic sym- 
metry requirements, preventing such unsatisfactory 
situations as the calculation of two different results for 
different but equivalent spatial orientations of a 
system. 


Much effort has been devoted to finding an 
optimum or clearly superior direction set from among 
the many possibilities within the triangular, symmetry 
preserving sets. Such a set does not seem to exist. For 
practical computational purposes, it makes little dif- 
ference which of the standardly available sets is used, 
provided the same set is used consistently. In shielding 
or other applications where the angular flux becomes 
collimated in streaming, it may be advantageous to use 
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especially tailored quadrature sets. For example, if the 
emerging angular flux is desired far from a localized 
source, it is known in advance that most of this flux 
will be confined to a narrow band of directions. Hence 
most of the available directions can be clustered in this 
range. Sometimes special sets can be formed within the 
framework of standard point arrangements; otherwise 
minor modifications to programs can be made to 
permit “asymmetric” quadratures.°**** In principle, 
the discrete-ordinates equations can be solved for any 
direction sets, but programming treatment of progres- 
sion through the mesh and setting of boundary 
conditions may impose limitations. 


Convergence and Acceleration Options 


With the exception of a direct-inversion program 
for one-dimensional geometries listed below, all 
discrete-ordinates programs iterate until input con- 
vergence criteria are satisfied. These criteria vary from 
program to program. In general, the more integral the 
quantity being monitored, the faster the convergence. 
For example, kes converges more rapidly than the 
fission-energy distribution, which converges more 
rapidly than the spatial variation of the scalar flux, 
which converges more rapidly than the pointwise 
angular flux. Various programs monitor different 
quantities. Some have options for monitoring point- 
wise convergence or integral quantities as desired. 
Others provide for convergence only in spatial zones of 
interest. A particularly constraining form of con- 
vergence criterion insists that coarse-mesh balance be 
satisfied for the whole system. 

For acceleration of convergence, most discrete- 
ordinates programs apply some form of the particle 
conservation rebalancing procedure discussed in the 
preceding section. Most often this rebalance is per- 
formed for the whole system in each group to 
accelerate inner iteration and for the whole system to 
accelerate upscatter convergence. Other possibilities 
include coarse-space-mesh rebalance for inner itera- 
tions, and coarse-group rebalance for upscatter, fission 
eigenvalue, or fission plus inhomogeneous source 
problems.”°’?> Further refinements to pointwise space 
rebalance*® and groupwise energy?°’*” rebalance are 
available. 

In addition to rebalancing accelerations, use is also 
made of Chebyshev acceleration?’’** and more con- 
ventional extrapolation procedures.*7"*? The 
Chebyshev acceleration works well in one-space dimen- 
sion but failed catastrophicaliy when applied in one 
two-dimensional program.'? 
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Eigenvalue Options 


Because of their common ancestry, most discrete- 
ordinates programs allow calculation of kegs or system 
dimension, concentration, or time eigenvalues to 
achieve desired values of Kee. A “time” or “alpha” 
eigenvalue is the value of the exponent in an assumed 
exponential time dependence of the flux. Most pro- 
grams include provisions for representing leakage in 
transverse directions by bucklings, and some are 
capable of searching for values of the buckling required 
to establish a given Kerr. 


Special Capabilities 


Nonlinear problems can also be treated by discrete- 
ordinates methods. In photon transport, scattering 
rates are enhanced by the photon population in the 
final state. Such effects can be included in the 
discrete-ordinates equations.*® Other nonlinearities, 
such as feedback mechanisms, can also be introduced 
in the equations, As long as nonlinearities can be 
treated as sources and calculated iteratively, there 
seems to be little difficulty in including them. 

In time-dependent one-dimensional geometries, it is 
relatively easy to allow the time variation of the system 
description, for example, to allow boundaries to move, 
compositions to change, and so on. Thus discrete- 
ordinates computations coupled to Lagrangian hydro- 
dynamic calculations are straightforward. In two 
dimensions, moving boundaries do not usually conform 
to orthogonal grids, so the situation is more difficult. 

In general, most discrete-ordinates programs use 
the same space and angle mesh for all energies, which is 
inherently wasteful. It is easily possible to allow the 
angular resolution to be a function of energy. Most 
commonly, a switch to diffusion theory is allowed.*' 
Similarly it is possible (because the inner iterative cycle 
is based only on moments of the flux) to gradually 
increase the number of directions used in inner 
iteration as the iteration progresses to save effort and 
possibly accelerate convergence.*? Also possible, al- 
though less simple, is the changing of the space mesh as 
a function of energy. In the same category is “‘line 
deletion” or the use of local grids. At present, a fine 
space resolution in one region requires mesh lines that 
run throughout the system, again wasteful. An easily 
specified, automatic procedure for allowing variable 
grids would be highly desirable. 

Most discrete-ordinates programs require cross sec- 
tions to be regionwise constant, although a region may 
be as small as a single mesh box. Many programs allow 
cross sections to be multiplied by a pointwise func- 
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tion,’"*' in effect, allowing a continuous density 
variation of cross sections. 

Charged-particle transpor or other transport 
problems with an acceleration term in the transport 
equation can be treated with discrete-ordinates pro- 
grams. For example, it is relatively easy to introduce a 
speed derivative into the finite-difference algorithm. 


Specific Programs 


Many different numerical programs with approxi- 
mately the same capabilities have been written to solve 
the discrete-ordinates equations. To differentiate be- 
tween programs, one will find it helpful to understand 
their genesis. Figure 8 shows the genealogy of the 
codes developed since the method originated at 
Los Alamos. Some of the early programs were written 
in very efficient but special machine languages, such as 
FLOCO.** These languages hid the inner workings of 
the program and kept those who understood the 
physics of particle transport separate from those who 
understood the numerical solution techniques. With 
the gradual emergence of FORTRAN as a common 
language, this barrier disappeared, leading not only to 
the proliferation of new programs but also to the local 
modification of programs for special applications. The 
programs DTF-IV’ and ANISN,*! which have been 
used at installations around the world, exist in dif- 
ferent versions at each installation. In one dimension, 
an up-to-date version of either of these programs 
should satisfy most needs. ANISN was developed in an 
IBM computing-system environment, and DTF-IV 
evolution began on IBM machines, switched to the 
UNIVAC-1108, and continued on CDC-6600/7600 
computers. ANISN offers certain options (e.g., shell 
sources and group selection of diffusion theory) not 
available in DTF-IV and in certain guises provides more 
sophisticated input and output options. DTF-IV has 
perhaps been seasoned longer in a wider variety of 
applications but, except for more sophisticated ver- 
sions like 1DFX,°’ is probably somewhat less eco- 
nomical to use than ANISN. A selection between the 
two might well hinge on how easy the initial acquisi- 
tion would be. 

In Table 2 most of the “‘last-generation” one- 
dimensional programs of Fig. 8 are categorized in terms 
of options. The programs XSDRN** (ANISN) and 
GTF*® (DTF-IV) are programs for calculating space- 
dependent spectra for use in forming multigroup cross 
sections. The speed of one-dimensional calculations has 
progressed to the point that few hundred-group cal- 
culations, even in the thermal range, are feasible. 
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One-dimensional time-independent programs 









































SNAP/SNAG 
LASL-1953 
| 
SNG 
LASL-1958 
| 
DSN 
NIOBE LASL-1959 MIST 
UNC-1960 | + 1 PHILLIPS PETROLEUM 
INTERNUCLEAR-1961 
bi eal DTK GAPLSN MGSN = 
5 i LASL-1961 GA-1961 PRATT & 
GE-1962 1 WHITNEY — ‘966 ee 
f 1 coe CO.-1961 
SNARGID DTF DTF-II 
BE21  ANL-1965 —LASL/UNC-1963 Al-1963 
BAPL-1969 | F 
SNARCID DTF-IV ANISN 
ANL-1970 LASL-1965 ORNL-1965 
| | 
e T T a f ar 1 
DTFXRAY DTF69 GTF IDFX SNID XSDRN ASOP 
NWEF-1968 SANDA-1969 GGA-1968 GGA-1971 GE-1968 ORNL-1969 ORNL-1969 
Two-dimensional time-independent programs 
TDC 
LASL-1959 
t 1 
DDK 2DXY 
LASL-1963 AEROJET-1961 
I 1 
DDF/2DF ODD-K TDSN (also ID) 
LASL-1965/UNC-1964 WANL-1966 NASA-LEWIS-1966 
! : T q 
DOT SNARG 2D TWOTRAN (x,y) IT 
ORNL-1966 ANL-1966 GGA-1967 KAPL-1968 
! 1 
DOT Iti TWOTRAN 
ORNL-1969 LASL-1969 
{ T CS f I tT i ig 1 
DOT-II1 DOT-2DB DOT-IIW TWOTRAN-V TWOTRAN-FC TWOTRAN-PN TWOTRAN-SPHERE 
ORNL-1971 GE-1969 WANL-1970 LASL-1970 LASL-1971 LASL-1971 LASL-1970 
Al ATOMICS INTERNATIONAL 


ANL ARGONNE NATIONAL LABORATORY 
BAPL BETTIS ATOMIC POWER LABORATORY 
GA GENERAL ATOMIC 

GE GENERAL ELECTRIC 

GGA GULF GENERAL ATOMIC 

KAPL KNOLLS ATOMIC POWER LABORATORY 
NWEF NAVAL WEAPONS EVALUATION FACILITY 
ORNL OAK RIDGE NATIONAL LABORATORY 
LASL LOS ALAMOS NATIONAL LABORATORY 
UNC UNITED NUCLEAR CORPORATION 

WANL WESTINGHOUSE ASTRONUCLEAR LABORATORY 


Fig. 8 Ancestry of one- and two-dimensional discrete-ordinates transport programs. Chart is 
reasonably complete but not exhaustive. Year of development is approximate. 
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A similar categorization is given for two- 
dimensional programs in Table 3. The choice of most 
versatile programs is between a version of DOT*” or 
TWOTRAN,’ and again depends to some extent on the 
ease with which the program can be acquired. Both are 
large programs making extensive use of peripheral 
storage, overlays, etc., so that local computer-system 
configurations are much more important than for the 
one-dimensional program. DOT was developed for and 
has been used extensively in shielding applications. By 
contrast, TWOTRAN originated in a thermal reactor 
environment and includes many convergence accelera- 
tion features. The SNARG2D program is greatly 
limited in capability. Its limitations illustrate the 
developmental trend toward general anisotropic scat- 
tering, sophisticated convergence acceleration algo- 
rithms, negative-flux corrections, and generalized 
source options. 

Finally, in Table 4 we tabulate information for 
time-dependent discrete-ordinates programs. 

All the programs noted in Tables 2 to 4 compactly 
and automatically store data for each problem to 
minimize fast-memory usage. This feature permits 
problems with any combination of input parameters to 
be executed, provided total storage does not exceed 
fast memory. As noted in the tables, others do even 
better, automatically utilizing peripheral storage. 


USAGE GUIDE FOR DISCRETE- 
ORDINATES PROGRAMS 


Suppose you have obtained a discrete-ordinates 
program from one of the code centers or direct from 
the author. You have made the changes necessary to 
implement the program on your computer and success- 
fully executed the test problems that accompanied the 
program. What next? A discussion of this question and 
advice about the general usage of such programs 
follow. 


Expectations 


One of the first things to do is to calibrate the 
program in terms of cost, reliability, and accuracy. You 
probably already know how well your other programs 
do; running several problems with your discrete- 
ordinates program will give you a feel for relative 
running time. Because of the convergence insensitivity 
of discrete-ordinates methods to number of space 
intervals or number of directions, it is easy to estimate 
the cost of subsequent computer runs once a problem 
of the applicable type has been run. Execution times 
are almost directly proportional to the total number of 


K. D. Lathrop 123 


space intervals or directions. Provided the addition of 
groups does not alter the upscattering or fission 
coupling, execution time is proportional to the number 
of groups. 

Checking accuracy usually requires that the num- 
ber of both intervals and directions be increased. If the 
arithmetic mean relations are used, doubling the 
number of either should reduce the error by a factor of 
4. It is also possible to estimate absolute errors in this 
fashion. It is also necessary to check convergence 
precisions to see how much solutions change. A too 
tight convergence wastes much machine time. Once 
acceptable accuracy has been established, it should be 
maintained for all problems of a parameter survey. For 
example, if fuel zone widths are being varied, the same 
relative number of space intervals per zone should be 
maintained. 

While becoming familiar with a discrete-ordinates 
program, get to know the physics of particle transport 
as well. The unaccelerated inner iterative process 
consists in accumulating the solution by adding to- 
gether the flux due to no collisions, one collision, two 
collisions, etc. The use of any negative-flux correction 
recipe forces more particles from the affected cell. 
Such forcing may be acceptable when shielding situa- 
tions are considered because it can lead to an over- 
estimate of the escape probability, but may not be so 
satisfactory for predicting breeding ratios. Learn the 
idiosyncrasies of flux streaming in voids in curved 
systems which can either enhance or deplete a void 
central flux, depending on the flux angular distribution 
at the edge of the void. Learn to use the correct 
boundary conditions. Programs are dumb and will 
accept specular reflection as an outer boundary condi- 
tion in cylinders or spheres. The programs will also 
return a solution with infinite gradient in such condi- 
tions, although a zero gradient solution is probably the 
one you seek. (This pathology leads to the develop- 
ment of the “white” or isotropic return boundary 
condition for cell calculations.°°) Learn to choose 
angular direction sets pertinent to a problem. Some 
sets in (x,y) geometry have too many directions 
pointing in the transverse (z) direction and will 
underestimate leakage from such systems. 

Once a general familiarity with the program is 
acquired, most problems can be run acceptably. There 
will always be exceptions. Common problems are 
discussed below. 


Convergence 


The more absorbing and leaky a group, the faster it 
will converge by the ratio of sources proportional to 
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the flux divided by losses proportional to the flux. In 
thick homogeneous zones, this ratio approaches the 
within-group scatter cross section divided by the total 
cross section. The closer the ratio to unity, the more 
difficult a group is to converge. In fact, in “alpha” 
eigenvalue searches [exp (at) time dependence], this 
ratio can exceed unity when a is negative and no 
convergence is possible. 

The more local the convergence criterion, the more 
slowly convergence will be obtained. We have earlier 
noted the oscillations due to the arithmetic mean 
difference equations. Either these alone or in conjunc- 
tion with negative-flux fixups can produce slowly 
damped oscillations as a function of the iterative cycle. 
Hence a pointwise convergence criterion may converge 
slowly, whereas an integral criterion may converge 
rapidly. 

Certain acceleration algorithms can delay or pre- 
vent convergence. Simple overrelaxation can fail cata- 
strophically if the relaxation factor is not kept within 
bounds. The pointwise rebalancing algorithm such as 
used in ANISN*® can fail if not modified as indicated 
by Reed.?? For this reason the algorithm is applied 
only every third iteration in ANISN. The Chebyshev 
acceleration in the earlier TWOTRAN (x,y) program®* 
can prevent convergence. 

Convergence may be prevented by a too-tight 
criteria. On a UNIVAC-1108, a criteria of 10° 
convergence will fail because the result of subtracting 
two nearly equal numbers may differ by this much. An 
outer iterative precision should be larger than an inner 
precision because the error in the outer cycle can be 
the cumulative error of all the inner cycles. 

For the above reasons, most discrete-ordinates 
programs have provisions to limit the number of inner 
iterative cycles. Such provisions not only prevent 
endless iterations but they can also be used to save 
computation time. In any problem expected to run 
through several outer cycles, no more than 10 inner 
cycles should be allowed per group. As the outer cycle 
progresses, convergence of the inner cycles is also 
obtained without undue iteration in early outer cycles. 
When convergence is limited, most programs monitor 
the ultimate residual error, allowing a user to decide 
whether more iteration is necessary or not. 

When one shifts to one program from another, it is 
important that he be familiar with the effectiveness of 
the different criteria. TWOTRAN’ criteria are often 
two orders or more of magnitude more effective than 
those of DDF.°' The DOT program has very useful 
“zone-of-interest” criteria?* allowing selective con- 
vergence. All programs have a hierarchy of convergence 
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precisions, and, if the user has the option, he should 
decrease the required precision in the general order of 
local to global. For instance, inner cycles should be 
converged most tightly, then upscatter or rebalance 
error for the whole system, and then eigenvalue errors. 


Negative Fluxes 


If the program has a fixup or uses a positive 
difference scheme, the only legitimate source of a 
negative flux is the existence of a negative source. 
Negative sources can arise due to the use of a truncated 
spherical harmonic scattering expansion. When this 
expansion is evaluated at the discrete directions of the 
problem at hand, it may be negative. Such expansion 
occurs frequently in first-collision source problems in 
which the uncollided flux is often completely col- 
limated. Consequently its spherical harmonic moments 
do not decrease, and the scattering expansion for the 
collided flux is likely to be negative. Normally only 
angular fluxes are negative in such situations, and the 
angle-integrated flux remains positive, especially the 
total of collided and uncollided fluxes. Usually nega- 
tive fluxes from such causes can be tolerated. 

Another reason for negative fluxes is an error in 
input cross sections. A mispunched exponent in a 
scattering transfer matrix can wreak havoc with the 
size of a scattering source. This leads us to the most 
frequent source of unsatisfactory problem runs. 


Input Errors 


Although most discrete-ordinates programs use 
sophisticated loading routines to count, check, and 
print input data and although most programs check the 
input to guarantee that nonphysical or nonsensical 
options are not specified, input errors remain the 
largest cause of problem failure. Not every type of 
error is anticipated by the checking routines. For 
example, it is assumed, not always correctly, that users 
specify positive convergence precisions. 

Certain classes of input errors can be determined 
from the printed output of the calculation. Advisory 
messages are printed by checking routines whenever 
cross-section sets are inconsistent in the sense that the 
absorption implied by the input total cross section and 
scattering matrix does not correspond to the indepen- 
dently entered absorption cross section. Balance tables 
are prepared during the course of the calculation and 
are a sensitive measure of cross-section inconsistencies. 
If the group sum of the outscattering rate does not 
equal the group sum of the inscattering rate, while at 
the same time group neutron balance is not in error, 
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then some cross section is incorrect (or there is code 
error). The neutron-balance monitor also provides 
clues. It is a global measure of the equality between 
losses and sources. Balance suffers when a group or 
groups are not converged or when, for example, a 
high-order anisotropic scattering source is used with a 
low-order discrete quadrature, which will not correctly 
integrate the spherical harmonic coefficients of the 
scattering source and will thus cause an erroneous 
addition to the balance equation. 

Balance tables and summaries are excellent means 
of monitoring changes to discrete-ordinates programs. 
Neutron or particle balance is extremely sensitive to 
any error in formulation; consequently the balance 
tables provide excellent diagnostic aid in debugging. 


APPLICATIONS 


In this section several examples of the application 
of discrete-ordinates methods are discussed. These have 
been selected from problems that show a variety of 
different, but typical, uses of the method. Because the 
examples are ones with which the author has been 
associated, the programs used are those available to 
him. 


Calibration of the Method: 
A Shielding Bench-Mark Problem 


In addition to executing test problems and com- 
paring results with those from similar programs, it is 
extremely useful to measure the performance of the 
method in real problems where exact solutions are 
known. Naturally it is very difficult to pose realistic 
problems for which exact or very accurate solutions are 
available. The Mathematics and Computations Division 
of the American Nuclear Society has compiled a book 
of such bench-mark solutions®? of which the problem 
below not only models a realistic shielding situation 
but also provides a means for assessing the severity of 
the ray effect. The problem is a two-group (x,y) 
geometry model of a central source emitting in a highly 
absorbing medium. The geometry for the problem is 
shown in Fig.9, and the two-group material cross 
sections including linearly anisotropic scattering are 
given in Table 5. The portion of the system shown is 
approximately 13 by 14 mean free paths in extent, and 
the secondaries ratio, 0;/0;, is a highly absorbing 0.16 
in the first group. A P; 9 evaluation of the exact-kernel 
solution for the infinite medium model of this problem 
is available for comparing solutions. This solution is 
applicable everywhere except near the right boundary, 
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Fig.9 Geometry and boundary conditions for bench-mark 
problem. Solutions are compared along the lines y = 80 and 
y= 140. 


Table 5 Cross Sections and Sources 








Group 1 Group 2 
Isotropic, cm” 
Xa 0.061723 0.096027 
vif 0.0 0.0 
=r 0.10108 0.108529 
Zsogg 0.015923 0.012502 
Zsog— 1g 0.0 0.023434 


Linear Anisotropic, cm™ 
on 0.008976 
Ysig—1-g 


0.003914 
0.009016 


Source Density, neutrons/cm? 
0.0065460 0.017701 





where a vacuum boundary condition is specified. This 
enables the effect of the boundary to be assessed. 

In Fig. 10 we compare the x variation of an S; 
discrete-ordinates solution to this bench-mark problem 
to the x variation of the kernel solution along the line 
y = 80. The approximate solution is naturally lower 
than the infinite medium solution near the right 
boundary; otherwise it is close to the kernel solution 
but oscillates about it. These oscillations are due to the 
ray effect. The further from the source, the larger the 
distortion, as evidenced in Fig. 11, where the approxi- 
mate solution is compared to the exact solution near 
y = 140 or the top edge of the system. Also shown in 
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Fig. 10 Percent deviation of S,, solution from “exact” 
solution at y = 80 cm, Exact solution is for infinite medium so 
that S,, solution, which was computed with a vacuum 
tight-boundary condition, is correctly smaller. Oscillations are 
ray effects. 


Fig. 11 is a spherical harmonic-like solution of order P, 
obtained from a modification of the discrete-ordinates 
method.’! Figure 12 shows the same two solutions 
compared to the kernel solution, except that a reflect- 
ing boundary condition was used with the P; solution. 
In this case the P,; solution deviates less than 3% from 
the kernel solution.* 

In assessing these results, one should note that the 
flux in group one falls by more than four orders of 
magnitude as a function of y from the bottom of the 
system to the top near x =0. The flux at the top 
further decreases by more than a factor of 30 as a 
function of increasing x from x=0 to the right 
boundary. Consequently the problem is a severe test of 
discrete-ordinates methods in a shielding-type situa- 
tion. Except, of course, near the right boundary, the 
S;2 flux solution is nowhere more than 18% in relative 
error, which is remarkably good agreement. Similar 
results were obtained for the second group. 

The problem as a whole provides an excellent 
check-out of a multigroup, two-dimensional, 
anisotropic-scattering transport program. 





*The P, solution was obtained from a modified version of 
TWOTRAN’ using a discrete-ordinates approach.'7 Much 
more arithmetic and iterative effort is required for the P, 
solution. 
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Fig. 11 Percent deviation of S,, and P, solution from 
“exact” solution at y = 140cm, S,, solution shows increased 
ray effect, whereas P, solution eliminates the effect. (Oscilla- 
tions in the P, solution are due to space differencing.) 
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Fig. 12 Same comparison as Fig. 11 except P, solution was 
calculated with a reflecting boundary to provide a direct 
comparison with “exact” solution. 


Reactor Core-Design Optimization 


In the preliminary design of a reactor core, many 
parameters are at the core physicist’s disposal. To 
optimize performance, he is likely to perform a series 
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of cheap, one-dimensional transport calculations. 
Typical of such a design effort are the calculations 
performed and described by Adler®? for the Kinetic 
Intense Neutron Generator (KING) reactor.°* In this 
reactor type the aqueous fuel solution is cooled after 
passing through the active core. For a beryllium- 
reflected, light-water-moderated version of _ this 
system,°* Adler varied the fuel concentration and 
central-flux-trap radius in a series of one-dimensional, 
cylindrical geometry, 16-group Ker critical zone thick- 
ness and concentration search problems. The results 
from these problems permitted the determination of 
the optimum value of the Feinberg®* performance 
index, as illustrated in Fig. 13. The most complex of 
these calculations required no more than a few minutes 
computation (CDC-6600). With the advent of the 
CDC-7600, such parameter surveys are becoming fea- 
sible for two-dimensional calculations. 
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Fig. 13 Variation of performance index with flux-trap radius 
and fuel-solution concentration. 


Reactor Control-Rod Worth and 
Radiation Heating Calculations 


For the Kinglet reactor configuration shown in 
Fig. 14, Malenfant®® has recently reported results of 
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Fig. 14 Kinglet reactor calculation configuration. 


one-dimensional neutronic and gamma-ray heating cal- 
culations and of two-dimensional control-rod worths. 
The composition of the system is shown in Table 6. 
The comparatively routine heating calculations pro- 
duced axial average neutron and gamma heating rates 
for each material for 36 radial intervals. These results 
are summarized by zones in Table 7. 

The more unusual control-rod-worth and central- 
void-worth calculations were performed by Malenfant 
using the TWOTRAN’ program and were compared 
with experimental results obtained by Wimett.°° In 
the experiment the control-rod worth for a 0- to 8-in. 
roa insertion was determined by integrating differential 
worths. During the experiment the system was kept 
near critical by varying the depth of a 1 -in.-diameter 
central void. The calculations indicate that the cumula- 
tive measured rod worth for the full 8-in. insertion would 


REACTOR TECHNOLOGY, Vol. 15, No. 2, Summer 1972 





130 DISCRETE-ORDINATES METHODS FOR SOLUTION OF THE TRANSPORT EQUATION K. D. Lathrop 











Table 6 Kinglet Material Description sured rod worth. Similarly the calculated worth of the 

35.5-cm, 1-in.-diameter central void in the presence 

Composition and absence of the control rod should bracket the 

Region atom/barn-cm g/cm* experimentally measured void worth. The agreement 
between calculation and experiment is excellent. 

Core H, 0.0648 mM, Goes In Fig. 15 the axial shift of the power density due 

O, 0.0351 O, 0.93239 : : 3 : 

U_ 0.000216 U. 0.084349 to control-rod insertion is shown. Figure 16 shows a 

7. abe Zr, 0.0429 Zr, 6.497 contour plot of the two-dimensional flux in the 

Control rod Cd, 0.0162 Cd, 3.02 thermal group. Note the severe perturbation caused by 

Fe, 0.0345 Fe, 3.20 the control rod and the interesting profile of the 

Reflector Be, 0.122 Be, 1.825 reflector flux peak. Contour plots and projective plots 





of fluxes are a standard option of TWOTRAN. 


























Table 7 Summary of Kinglet Radiation Heating l ! l | | ] 
Normalized to 10'® Fissions 3.0/— 4 
Rod in 8 in.—+ 
Gamma-energy Neutron-energy Total energy ud 
AS a ai = 
deposition, deposition, deposition, a 
Region MJ MJ MJ - 
S 2.0 
Core 0.494 0.421 0.915 o 
Zr tube 0.261 0.261 ” 
Be reflector 0.825 0.250 1.075 
ie re Eee aint tes lJ 
Total 1.580 0.671 2.251 a 
- 1.0 
aq 
P| 
w 
a 
be a function of void position, due to a displacement 
of the flux peak from the midplane as the rod was : 
inserted. The experimental results were not corrected 0 
for this effect. Both calculational and experimental 
results are summarized in Table 8. The worth of the 
fully inserted rod was calculated with no void and with Fig. 15 Kinglet axial power-density distributions for control 
the full void. These two results should bracket mea- rod fully inserted and completely out. 


Table 8 Kinglet Calculational and Experimental Results 





Calculated Experimental 





Calculated worth, worth, 
Configuration* keff dollars dollars 
Rod out, no void 1.05397 
Rod in, 1-in.-diameter void 1.04249 
Rod in, no void 1.01866 
Rod in, full 1-in.-diameter void 1.00271 
Rod out, 1-in.-diameter void 1.44 
1-in.-diameter void with no LS 
correction for void—rod 
interaction 
Rod in, 1-in.-diameter void 1.99 
8-in. rod insertion, no void 4.41 
8-in. rod worth with no correction 5.00 
for void—rod interaction 
8-in, rod insertion, 1-in.-diameter void 4.97 





*1-in.-diameter void is understood to extend from top of core to z = 35.5 cm. 
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Fig. 16 Kinglet thermal-group-flux contour plot showing 
effect of control rod and location of reflector-flux peak. Not 
all contours are drawn in the vicinity of the control rod. 


Pulsed-Neutron Calculations 


Suppose a point burst of neutrons is produced on 
the axis of a finite cylinder. What is the space—time 
evolution of the flux? Such questions are answered by 
time-dependent discrete-ordinates calculations, an 
example of which is shown in Fig. 17. Here the time 
dependence of the leakage and absorption rates from a 
discrete-ordinates calculation is compared with Monte 
Carlo results. In this problem a point-burst source was 
located at p=0, z=1.9 inside a cylinder of 1-cm 
radius and 2-cm height in a medium with o,/o, = 0.2. 
Inflections in the curves of the figure are visible at 
times associated with the uncollided wave front reach- 
ing the bottom surface, reaching the right surface, and 
leaving the system. A more complete description of 
this calculation is given in the manual describing the 
TRANZIT program?’ used to perform the discrete- 
ordinates calculation. A similar calculational result 
obtained with this program is shown in Fig. 18, where 
contours of energy deposition in an exponential 
atmosphere are plotted for the time 1 to 2 sec after a 
point burst of neutrons at a height of 26 km at ¢=0. 


Shielding Survey Calculations 


As part of a design effort, it was desired to 
calculate the flux of photons from a ®°Co source and 
other gamma sources located at one end of a cylindri- 
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Fig. 17 Time variation of leakage and absorption rates after 
point burst of particles on cylinder axis. 
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Fig. 18 Energy-deposition contours for time interval 1 to 
2 sec after point air burst of particles. Shape of contours 
reflects the effects of varying air density. 


cal container immersed in water. In particular, the 
objective of the calculations was to shield a small 
region inside the cylinder (housing electronic gear) 
subject to strict weight and size limitations. Figure 19 
shows a picture of one such problem in the sequence as 
it is drawn by the computer output of the TWOTRAN 
program.’ The location of the source and detector 
regions has been added to the figure. An iron barrier is 
interposed between the source and detector, and a 
uranium can surrounds the detector. Water surrounds 
the whole assembly, and a white boundary condition is 
used to ensure an overestimate of the reflected flux 
from the water. This problem was run with 28 radial 
intervals and 88 axial intervals in a three-group 
representation. An Sg quadrature was used. This 
particular problem required about 2 min (CDC-7600) 
of computer time. Although this is very cheap, many 
such calculations were run since various barriers, 
encapsulating configurations, and source spectra were 
tried. To save calculational effort, the following steps 
were taken: 

1. Only one calculation in a particular sequence 
was computed using anisotropic scattering. As it 
happened in this problem, isotropic scattering gave an 
overestimate (because of reflection from the water) of 
the detector flux. 

2. The water end caps were included in only one 
problem in each sequence. 


3.The positive difference scheme version of 
TWOTRAN was used to ensure smooth solutions and 
rapid convergence. As noted previously, this scheme 
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tends to overestimate the flux reaching the detector. 
which is a safe error. 


4.Space meshes and group structures were kept 
crude but consistent. Thus relative comparisons could 
be made without having to obtain absolute accuracy. 


SUMMARY 


The review has covered the theory, capabilities, 
limitations, and applications of discrete-ordinates 
methods and programs for the solution of the transport 
equation. In general terms, the method is more 
expensive but is a more accurate approximation than 
diffusion theory. The method is less geometrically 
versatile than Monte Carlo techniques but is usually 
more economical. For two- or three-dimensional 
geometries, and especially for time-dependent situa- 
tions, a Monte Carlo estimate of localized or integral 
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Fig. 19 Geometry for shielding survey calculations. Picture is 
a redrawn copy of computer output from TWOTRAN with the 
border of the figure indicating boundary conditions: 0, 
vacuum; I, reflective; and 2, white. The number of uniform 
fine-mesh intervals per coarse-mesh interval is shown in the left 
column and bottom row. 
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quantities may be less expensive than a discrete- 
ordinates calculation, but, whenever the detailed solu- 
tion is required in every volume element of the system 
(e.g., power distribution, gamma heating rate, neutron- 
energy deposition, and angular flux), the discrete- 
ordinates method is cheaper. It is also easier to perform 
parameter surveys with discrete-ordinates methods 
than with Monte Carlo programs. The ease with which 
both direct and adjoint solutions are obtained has led 
to the development of many programs to perform 
perturbation calculations with the output of discrete- 
ordinates programs. 

In discussing limitations and usage of discrete- 
ordinates methods, our intention has been to create 
awareness of the pitfalls rather than to inspire mistrust. 
Used with care and intelligence, the method has been, 
and will continue to be, an extremely useful tool. 
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Development of Materials and Systems 
for the Molten-Salt Reactor Concept 


By H. G. MacPherson* 


Abstract: /nterest in molten-salt-fueled reactors stems largely 
from the ease and relatively low cost of purifying and 
reconstituting the fuel. Molten-salt thermal breeder reactors — 
graphite moderated and with salt constituents that absorb few 
neutrons—will have low enough inventories and doubling 
times to greatly ease the future requirements for uranium 
mining and for enrichment facilities. 

Among the reactor designs that have been studied, the 
favorite is a single-fluid breeder fueled with a mixture of ’ LiF, 
BeF,,, ThF,, and **? UF,,. The chemical characteristics of such 
salts and the behavior of fission products in them are quite well 
understood from many years of research and from the 
operation of the Molten Salt Reactor Experiment (MSRE). A 
number of processes for removing the soluble fission products 
(and protactinium) from the salt have been explored in recent 
years. Of these, reductive extraction into bismuth is being 
developed most vigorously. A bubble generator and centrifugal 
gas stripper are being developed to ensure the speedy removal 
of the noble-gas fission products, which are only sparingly 
soluble. 

Unresolved problems receiving current attention include 
the control of tritium arising from neutron reactions with 
lithium, surface fissures in Hastelloy N probably arising from 
fission-product attack, the retention of the radiation resistance 
of Hastelloy N to high temperatures, and the development of a 
steam generator. 


The use of a molten salt as the fuel material has been 
proposed for many different reactor types, much as 
UO, is used as the fuel for a variety of different 
reactors. Conceptual designs have been prepared for 
ship and aircraft propulsion systems as well as for fast 
breeders and thermal reactors for central-station power 
generation. Nearly all of the recent development work 
has been directed toward a graphite-moderated thermal 
breeder. 
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The principal attraction of a molten-salt fuel is the 
ease of fuel preparation and reprocessing. This ease 
leads with certainty toward a low fuel-cycle cost. 
Furthermore, the possibility of rapid processing, in- 
cluding the rapid removal of '*5Xe, makes the 
molten-salt system the only one in which breeding with 
a reasonably short doubling time appears to be 
economically possible in a thermal reactor. 

The fuel salt for these thermal reactors is a mixture 
of 7LiF and BeF, in a mole ratio somewhat greater 
than 2:1, with the addition of as much UF,, PuF3, or 
ThF, as is required to make the reactor critical or to 
provide the desired amount of fertile material. For 
example, the composition of the salt used in the most 
recent conceptual design is 7 Lif -BeF, -ThF,—UF, in 
the mole ratio 71.7:16:12:0.3. It has a melting point 
of about 500°C (930°F) and a vapor pressure of less 
than 0.1 mm Hg at a mean operating temperature of 
620°C (1150°F). This salt has about 3.3 times the 
density of water, about 10 times the viscosity of water 
(11 cP at 620°C), and a thermal conductivity and heat 
capacity per unit volume comparable to water. These 
properties are favorable for heat transport; particularly 
important in this respect is the high heat capacity per 
unit volume which allows large amounts of heat to be 
transported with moderate salt flow rates. 

The high melting temperature is an obvious limita- 
tion, requiring careful design to keep the system hot. 
Concern is also frequently expressed about corrosion 
of the system, containment of the radioactive liquids 
and gases, reliability of the equipment, and main- 
tenance of the highly radioactive system. Operation of 
the 7.5-MW Molten Salt Reactor Experiment (MSRE) 
from early 1965 to December 1969 provided favorable 
first-order answers to these questions. MSRE operation 
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also called attention to two new problems. One 
problem is the necessity for controlling the tritium that 
is produced in large amounts from the lithium in the 
fuel. The other problem is a moderate grain-boundary 
attack on the nickel-based container alloy that appears 
to be associated with the presence of some fission 
product, probably tellurium. 

The present development work is directed toward 
providing the additional technology required for build- 
ing the much larger and longer lasting reactors needed 
for power production. In addition to the tritium and 
grain-boundary problems, particular current attention 
is being devoted to the development of chemical 
processing methods, improvement of the radiation 
resistance of the nickel-based alloy, and the gaining of 
confidence in the handling of the sodium fluoroborate 
salt that has been proposed for use as the intermediate 
coolant to transfer heat from the fuel salt to the steam 
generator. 

An extensive review of molten-salt reactor technol- 
ogy at Oak Ridge National Laboratory (ORNL) 
through mid-1969 is given in a series of articles in the 
February 1970 issue of Nuclear Applications and 
Technology.'* Included in this series is a discussion of 
the rationale for the development of molten-salt 
reactors, showing that they, like the Liquid-Metal- 
Cooled Fast Breeder Reactor (LMFBR), have the 
capability of solving our long-range energy supply 
problem.' Further progress in ORNL’s development 
program and the current reference design of a molten- 
salt breeder are described in Ref.9. Other relatively 
recent articles giving a broad view of molten-salt 
reactors and citing their virtues are those listed in 
Refs. 10 to 14. 

The sections that follow give a somewhat historical 
account of the major reactor configurations and design 
concepts that have been proposed using molten-salt 
fuels, followed by descriptions of the recent progress 
and the current status of important areas of molten-salt 
reactor technology. 


REACTOR CONFIGURATIONS 
AND CONCEPTS 


Since the start of the Molten Salt Reactor Program, 
conceptual designs, revised frequently, have guided the 
development program. The first designs were for 
aircraft propulsion, taking advantage of the high- 
temperature possibilities.'* At the start of the civilian 
power program, homogeneous reactors were consid- 
ered, the only moderation being that provided by the 
fluorine and beryllium in the salt.’® With such 
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reactors, breeding was not quite possible, and only 
converter reactors were considered. By late 1958 
enough experimental information about graphite in 
contact with fluoride fuel salt in a radiation environ- 
ment was available to shift the design effort to graphite 
moderation. It was quickly found that the optimum 
performance usually occurred with 10 to 20% of the 
core volume occupied by molten-salt fuel and the 
remainder by graphite. At that time the most feasible 
reactor still seemed to be a converter, partly because 
complete chemical processing procedures were un- 
known but also because it was not easy to find 
plausible designs for two-region reactors. 

For breeding, the two-region reactor seemed neces- 
sary to prevent leakage from consuming too many 
neutrons. It was realized, though, that, if the reactor 
were made big enough, say a 20-ft-diameter right 
cylinder, the leakage could be kept to about 0.05 
neutron per neutron absorbed in fissile atoms, a value 
low enough to allow breeding if rapid processing 
schemes could be developed.'’ Such a large reactor 
seemed impractical because, with its large inventory, it 
would be suitable only for power levels of 1000 MW or 
greater, a size unheard of at that time. There was a 
tendency to discount the value of reflector regions 
because they would have to be cooled with fuel salt 
and would generate some power and some neutrons. In 
later years all these prejudices were overcome. 

In the early 1960s the major conceptual design 
effort shifted to two-region two-fluid breeder reactors, 
with the 77°U fuel in one fluid and the fertile thorium 
in the other fluid. Reasonable designs were generated 
in which the fuel flowed through the core in graphite 
pipes and the fertile salt not only surrounded the core 
as a blanket but also permeated the core, bathing the 
outside of the fuel tubes and flowing through channels 
in the graphite moderator.'® With uranium and 
thorium in separate salt streams, the processing steps 
could be greatly simplified. A useful breeder could be 
designed using only the fluoride-volatility process, or 
that plus vacuum distillation to recover the carrier salt 
from the rare-earth fission products. An improvement 
on the performance of this two-fluid reactor was 
indicated when it was discovered that protactinium 
could be selectively removed from the blanket salt by 
reductive extraction into a bismuth carrier. 

A continuing concern with the reliability of the 
graphite plumbing in the two-fluid core led to more 
refined calculations of the possibilities of a single-fluid 
breeder. These calculations have evolved into the 
presently favored three-region single-fluid breeder.'? 
The salt that serves as the sing!e fluid contains both 
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thorium and enriched uranium, in a ratio of about 
50:1. A central core zone containing about 13 vol.% 
salt is surrounded by a “blanket” region with 37 vol.% 
salt. With less graphite, the blanket region is under- 
moderated, enhancing the resonance absorptions in 
thorium relative to the fission absorptions in uranium. 
This reduces the leakage, and a 2-in.-thick salt annulus 
at the outside edge of the blanket aids in this function. 
The third region is a graphite reflector about 75 cm 
thick which is cooled with the minimum amount of 
salt, about 1 vol.%. The reflector not only saves 
neutrons but it also reduces the fast flux at the reactor 
vessel. Concurrent improvements in chemical process- 
ing, described in a later section, made the breeding 
performance of this single-fluid system attractive. 

Capital-cost estimates for molten-salt reactors have 
been made a number of times, and they usually fall 
within the range of such costs for light-water reactors. 
The significance of the exact number is usually clouded 
by the inclusion of a rather large contingency, the 
magnitude of which depends on the confidence of the 
estimator. The capital cost estimated for a 1000-MW(e) 
three-region single-fluid breeder of the type described 
above was $203 million, based on January 1970 
construction costs. The number for a light-water 
reactor, on the same basis, was estimated at $201 
million. The direct cost of the molten-salt system is 
about $12,000,000 lower for the turbine plant equip- 
ment because of the higher temperature steam, but this 
advantage is used up in such items as an auxiliary boiler 
($3,000,000), maintenance equipment ($4,500,000), 
and extra allowance for contingencies ($6,000,000). 
The power production cost was estimated at about 5.5 
mills/kWh. 

Along the way to this latest design, a number of 
interesting modifications have been suggested. A simple 
one that is surely worthy of future experimentation is 
a two-fluid two-region reactor with the regions sepa- 
rated by a metal membrane.”° In this case the core salt 
would contain both uranium and thorium to allow 
some breeding in the core, but leakage from the core 
would be usefully absorbed in the blanket salt, which 
would contain thorium but no uranium. Such a reactor 
could have less leakage and a lower inventory than the 
single-fluid reactor, provided the metal membrane can 
be made sufficiently thin. 

A much more inventive modification involving lead 
cooling was suggested by Bettis.?' In this reactor the 
fuel salt is cooled by direct contact with a falling 
shower of lead droplets, and their momentum transfer 
is used to provide circulation of the fuel, thereby 
replacing the fuel pumps. A big advantage would be the 
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reduced inventory of fuel achieved by the elimination 
of the external heat-exchanger circuits. This would be 
particularly advantageous for a fast molten-salt reactor, 
which would use a chloride salt for the single fuel and 
fertile stream. Only a small amount of development 
work has been done on lead-cooled systems, but such a 
system is being considered seriously in the United 
Kingdom for fast reactor use.” ? 

Taube,?? first in Poland and now in Switzerland, 
has proposed fast molten-salt reactors using either 
aluminum chloride or mercury as a boiling coolant, in 
direct contact with the fuel. He has outlined the long 
list of development problems that must be solved to 
make such a system practical but sees great advantages 
in the compactness of the reactor and in its large 
negative temperature coefficient of reactivity. Another 
variant of the family of molten-salt reactors, the 
epithermal MOSEL concept, has been studied at 
Jiilich.?* In connection with this reactor, the pros and 
cons of internal vs. external cooling are discussed by 
Gat.?° 

The next reactor proposed for construction by 
ORNL is a breeder experiment, the Molten Salt 
Breeder Experiment (MSBE).?° The purpose of the 
MSBE would be to serve as a true pilot plant for a large 
breeder reactor. It must be capable of operation at 
power densities and fast-neutron fluxes at least as high 
as would be encountered in a power breeder and would 
incorporate facilities of continuous chemical process- 
ing. The components should be of such a size that only 
modest scaling-up would be required for the power 
breeder. For economy, the size should be as small as 
possible and still meet these requirements. This mini- 
mum size is determined primarily by the requirement 
for high power density and the desire to keep the 
neutron-energy spectrum and fuel concentrations rea- 
sonably close to those for the power breeder. A reactor 
of about 150 MW(t) is planned. 

Although it was desired that the MSBE have a 
breeding ratio approaching 1.00, the actual demonstra- 
tion of breeding does not seem as important as high 
power density from a technical point of view. The 
problem of achieving a breeding ratio above 1.00, still 
maintaining the neutron spectrum soft enough to be 
representative of the ultimate breeder, is that it would 
require a rather large low-leakage core. With the 
high-power-density criterion, this would mean a reactor 
of perhaps 250 MW(t) and extra cost. It should be 
possible to obtain an accurate check on the ultimate 
breeding performance of a large reactor by comparison 
of the calculated and actual performance of a MSBE 
whose conversion ratio is in the range 0.9 to 1.0. So far 
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the reactor-physics calculations seem to be in very 
good shape, as judged in part by the accuracy with 
which the critical concentrations of fuel were calcu- 
lated for the two loadings of the MSRE. The greatest 
uncertainty in the calculation of breeding is in the 
value of eta for 7**U; and the value used has been 
validated by measurements of alpha for 727U in the 
MSRE.?’ 

Recently ORNL has carried out a design study of a 
300-MW(e) “prototype” reactor.?*® The principal cri- 
terion for this reactor was that it should require as 
little development effort as possible and still be 
economically attractive as far as running costs are 
concerned. The purpose of the reactor would be to 
demonstrate the feasibility and delineate the problems 
of construction of a large molten-salt-reactor power 
station. The minimum development requirement led to 
the decision that it be a converter, without sophisti- 
cated chemical processing, rather than a breeder. Also, 
the core size was made large, about 6.4 m in diameter, 
with a correspondingly low power density, so that the 
graphite could be expected to last for the lifetime of 
the reactor. No attempt would be made to seal the 
graphite pores to keep out xenon, and the tritium 
problem would be handled by using a nitrite—nitrate 
salt in a third loop between the primary-coolant salt 
and the steam system. 

Recent reactor-physics and fuel-cycle studies have 
examined the effects of increasing the core size for 
larger reactors to reduce the graphite replacement 
problem.” For a 1000-MW(e) breeder, the principal 
effect of increasing the core diameter to 8.7 m (provid- 
ing a 30-year graphite life) was to increase the fissile 
inventory by 50%, and the fuel-cycle cost by 0.2 
mill/kWh. Other such studies have examined the use of 
plutonium as the startup fuel for a breeder and as 
makeup fuel for a converter. A tentative conclusion is 
that the use of plutonium rather than ?°°U gives a 
comparable fuel-cycle cost.?° The use of plutonium in 
molten-salt reactors is of special interest to the Indian 
Atomic Energy Commission, and they have started 
some work on the chemistry of PuF3 in molten 
fluoride salts. 

The success of the MSRE gave rise to an increase in 
industrial interest in molten-salt reactors. An industrial 
study group headed by Black and Veatch?! reported 
favorably on the promise of the system and proposed 
the construction of a prototype reactor that could be 
completed in 1976. A group headed by Ebasco 
Services, Inc., has made a similar study, and they have 
now been given the role of undertaking the major 
design effort in the United States, under contract to 
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ORNL. They have compiled a rather complete sum- 
mary and evaluation of molten-salt technology*? and 
have concluded that the existing technology is suffi- 
cient to justify the construction of a demonstration 
plant which could be made to work and from which a 
great deal would be learned about the cost of construc- 
tion and operation. 


THE FUEL SALT 


Previous reviews*’??*** have described how the 
presently favored molten-salt constituents and propor- 
tions were chosen. The mixture of 7 LiF and BeF, as 
the major brew to which the fissile and fertile species 
are added is probably unique in its combination of 
properties of low neutron cross section, stability 
against radiation damage and thermal decomposition, 
low vapor pressure, and tolerable melting point. Since 
the fuel is the heart of the reactor, the importance of 
understanding its basic properties and reactions cannot 
be overstressed. 

There are several ways of broadly characterizing 
the chemical state of a fused fluoride-salt mixture. The 
different ways are not completely independent, but 
each emphasizes or correlates certain aspects of the 
chemical behavior and physical properties. One of 
these ways for characterization arises from the relative 
proportions of the major constituents and is related to 
their tendencies to form complex ions. BeF, added to 
LiF tends to form the ion BeF3 , and ThF, tends to 
form the ion ThF3~. In each case the complex ions tie 
up or immobilize fluoride ions coming from LiF. By 
analogy with the silicate system, the BeF3} and the 
ThF3” ions can be thought of as “basic,” and any 
remaining unassociated (or free) F ions can be 
regarded as “‘acidic.”?5 

The excess or deficit in the free fluoride-ion 
content can be placed on a rough scale,?° as follows: 


Excess free fluoride = mole % LiF 
— 2(mole % BeF,) — 3(mole % ThF, ) 


This excess (or deficit) influences many of the proper- 
ties of the salt. For example, it has been shown to have 
a strong influence on the activity coefficients of 
various metallic-ion species dissolved in the salt mix- 
tures. This is particularly apparent in measurements of 
the ease with which bismuth containing lithium can 
reduce and extract rare earths from the fuel salt.>7 
Also, for example, it affects the solubility of BF; in 
the fuel since BF; tends to form the BF, ion, which it 
can do most readily in a salt with free fluoride ions.**® 
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A polymer model developed by Baes*? is related in 


concept. Starting with the observation that solid BeF, 
consists of the tetrahedral ions BeF%° arranged in a 
three-dimensional network, he studied the polymeric 
character of LiF—BeF, melts using activities derived 
from electromotive-force (emf) measurements and 
from liquidus data. In the concentration range of BeF, 
of interest for reactors, the model gave best fit to 
experimental data with most of the beryllium in the 
BeF3 ion, and with only small amounts of Be, F3~ ion. 

With respect to physical properties, salts with a 
high percentage of BeF, are undesirable because their 
polymeric structure results in a very high viscosity. For 
a recent summary of the physical properties of a 
variety of salt compositions, reference should be made 
to a report edited by Cantor.*° 

A second and perhaps more important way of 
characterizing the chemical state of a fuel-salt mixture 
is its effective potential for oxidation or reduction. 
This redox potential is directly related to the valence 
states that are assumed by those metal ions in the salt 
mixture that are capable of existing in different valence 
states. The greatest importance of the redox potential 
is that, to a large extent, it controls the corrosive 
tendency of the fuel mixture toward the Hastelloy N 
container alloy. The major constituents of this alloy 
are nickel, molybdenum, chromium, and iron; of these, 
chromium, constituting about 7% of the alloy, is by far 
the most easily oxidized. Accordingly, any corrosive 
tendency of the fuel salt acts selectively on the 
chromium of the alloy. 

Lithium, beryllium, and thorium in their fluorides 
LiF, BeF,, and ThF, have neither higher nor lower 
valence states. They cannot be oxidized and can be 
reduced only to the metals, and then only by reducing 
agents very much stronger than the constituents of 
Hastelloy N. Mixtures of these materials would not be 
expected to be corrosive, and this was confirmed by 
the almost immeasurably small corrosion of Has- 
telloy N by the LiF—BeF, coolant salt in the MSRE.*! 
Uranium, on the other hand, has several stable va- 
lences; its known fluorides include UF;, UF,, UF;, 
and UF.. Salts containing UF; and UF, have redox 
potentials so high as to be very corrosive and will be 
excluded from the circulating fuel system. In the 
absence of extraneous strong oxidants, such as HF, 
FeF,, and NiF,, the redox potential is determined by 
the relative concentrations of UF, and UF;. 

In this situation the corrosion occurs through the 
reaction 


Cr+ 2UF, = 2UF;3 + CrF, 
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with all fluorides remaining as solutes in the salt. 
Equilibrium is reached with quite small concentrations 
of UF; and CrF,, and the UF3/UF, ratio of an 
isothermal system can be adjusted so that the extent of 
reaction is negligible. However, the reaction is tempera- 
ture dependent, and there are small changes in the 
equilibrium concentrations as the temperature changes. 
In the reactor, therefore, there will be some surface 
depletion of chromium in the hot zones and formation 
of a surface layer rich in chromium in the cold zone. A 
satisfactory equilibrium is established? if the U**/U* 
ratio is above 0.005. In the MSRE, this condition was 
usually maintained (although not always), and the 
corrosion rate was about 0.1 mil/year, with some 
indication that much of this was associated with the 
admission of air during shutdown periods.*? 

The direct relation of the U**/U* ratio to the 
corrosion reaction suggests that it would be the ideal 
measure of the redox potential of the salt. As far as 
corrosion is concerned, the higher this ratio, the less 
corrosive the salt will be, but an upper limit, beyond 
which one should not go, is set by the disproportiona- 
tion reaction by which UF; is converted to UF, and 
uranium metal. This limit is rather high, corresponding 
to 20% or more of the uranium occurring? as UF3. In 
the presence of graphite, disproportionation can pro- 
duce uranium carbide instead of uranium metal, and 
recent experiments*? suggest that this reaction goes 
more readily, placing an upper limit of about 0.06 for 
the UF3;/UF, ratio to prevent formation of uranium 
carbide. There seems to be adequate leeway to main- 
tain the small ratio needed to keep corrosion low. 

In the MSRE, the U**/U** ratio was determined by 
the occasional analysis of a sample of fuel salt 
withdrawn from the reactor. The analytical method 
involved sparging the sample with hydrogen and 
measuring the amount of HF evolved as the various 
species were reduced.** Although the results obtained 
were useful, the method had many difficulties.4° Since 
then a relatively simple voltammetric method has been 
developed that seems suitable for automated in-line 
analysis.4® The method involves measuring the differ- 
ence between the equilibrium potential of the melt, as 
indicated by an inert platinum electrode, and the 
voltammetric equivalent of the standard potential of 
the U**/U** couple. The measurement has been auto- 
mated and is now undergoing test in a molten-salt 
thermal-convection loop.*” It gives promise of being 
applicable to an operating reactor. 

The fission process tends to increase the oxidation 
potential of the salt® since the active metal fission 
products are not produced in sufficient quantity to 
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match up with all the fluorides liberated by the 
destruction of UF,. Therefore it is necessary to adjust 
the redox potential to keep it within the desired range. 
In the MSRE, the redox potential was adjusted by the 
occasional addition of some reducing metal (beryl- 
lium), suspending it in a basket in the pump bowl until 
it dissolved.? Presumably the potential could also have 
been adjusted by introducing hydrogen into the gas- 
stripping system; indeed, this is the method of regu- 
lating the redox potential of the salt as it is prepared 
for use. The designs for breeder reactors include the 
continuous processing of a side stream of fuel so that, 
for these reactors, hydrogen treatment will probably be 
the normal method of adjustment. 

Long and Blankenship*® have carried out measure- 
ments of the concentration ratio of UF3/UF, that is in 
equilibrium with various mixtures of hydrogen and HF. 
From these data an equivalence of HF concentration in 
hydrogen to the U**/U* ratio can be obtained as a 
function of temperature. For example, at 600°C a 
U**/U* ratio of 0.01 in the molten salt is in 
equilibrium with hydrogen at | atm pressure contain- 
ing about 20 ppM HF. 

Baes*? has published a very useful summary of the 
thermodynamic information available in 1969 pertain- 
ing to reactions in molten fluoride salts. He gives 
equilibrium constants as a function of temperature for 
many heterogeneous reactions and indicates how they 
can be applied to the control of impurities in molten- 
salt reactors. Of particular interest to this discussion 
are his curves of the concentration of various corrosion 
and fission products to be expected as a function of 
the redox potential expressed as the UF4/UF; ratio. 

The redox potential also determines the balance 
between different valences of other metals. For exam- 
ple, if protactinium is in the salt, it can exist in either 
the 4+ or 5+ valence state. At the redox potential 
regarded as optimum for the reactor, protactinium 
exists as Pa**, but, by sparging the salt with HF, it can 
be oxidized*® to PaS*. This oxidation reaction is 
favored by increasing the temperature. 

During MSRE operation, there were short periods 
when the fuel was quite oxidizing, with the calculated 
U?*/U** ratio falling below 0.004. During these 
periods the fission-product niobium appeared in solu- 
tion in the fuel salt,°’ indicating the oxidation of 
niobium from Nb°® to Nb?* or Nb**. On one occasion 
the fuel was apparently sufficiently oxidizing for Fe?* 
to appear in solution. The presence of niobium in 
solution in the fuel could perhaps be used as an 
indicator of an excessively high oxidation potential, 
one at which corrosion could become a problem. 


The stability of niobium fluorides in molten salts 
is being investigated by Weaver and Gill.5? They 
find that NbF, is remarkably stable in the mixture 
2LiF—BeF,, much more so than as a solid. The 
niobium gradually disappears from solution at 700°C, 
presumably by disproportionation, producing NbF; 
gas. They conclude that niobium is slightly more noble 
than chromium and considerably more reactive than 
iron. 

The fuel salt has a limited tolerance for oxide ion, 
depending on the proportions of the major con- 
stituents. If the solubility for oxide ion is exceeded, 
metals will precipitate as oxides in a definite sequence. 
If protactinium is oxidized to the 5+ valence, it is the 
first to precipitate.6*'5* This unique position is being 
investigated as a possible processing step for separating 
the protactinium from the fuel.5* The tetravalent 
oxides are the next to precipitate; if there is a 
preponderance of ZrF4 over UF,, it will precipitate 
separately, before uranium or thorium. Advantage was 
taken of this fact in the MSRE, adding ZrF, to protect 
against the possible precipitation of UO,. If the oxide 
content is further increased, the uranium and thorium 
come down together as a mixed oxide, which is high in 
UO, until the uranium is depleted.*® Finally, if 
enough oxide is present, BeO is precipitated. Pu** does 
not precipitate as oxide from fuel salt even in the 
presence” of excess BeO or ThO3. 

The trifluoride species, PuF;, UF3, and the rare 
earths, have a limited solubility in fluoride-salt mix- 
tures and tend to coprecipitate when their solubility is 
exceeded. This naturally caused some concern over the 
use of plutonium as a fuel, so a number of measure- 
ments were carried out. Barton et al.°® measured the 
solubility of CeF3 in a large variety of molten-fluoride 
compositions between 550 and 800°C. In mixtures of 
the two constituents, LiF and BeF,, they found a 
minimum solubility for the composition with about 37 
mole % BeF,. The solubility was greatly enhanced by 
replacing BeF, with ThF,, presumably because the 
Ce** ion, with a charge-to-size ratio z/r of about 3, can 
compete for F ions better with Th** ions (z/r = 4) 
than with Be?* ions (z/r=6). The solubility also 
increased sharply in compositions with a very high 
value of the free fluoride-ion content. 

Measurements at ORNL*® and at the Bhabha 
Atomic Research Center in India®® have shown the 
solubility of PuF3 in the single-fluid breeder salt (72 
mole % LiF, 16 mole % BeF,, and 12 mole % ThF,) to 
be nearly identical to the solubility of CeF; measured 
by Barton. The solubility is about 0.8 mole % at the 
liquidus temperature of about 500°C and increases to 
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about 5.5 mole % at 800°C. These values are adequate 
for the use of plutonium as a fuel and would allow the 
accumulation of rare-earth fission products for as long 
as 10 years without processing if this were thought 
desirable. The measurements of PuF 3 solubility are 
being extended to other compositions by the Bhabha 
Atomic Research Center. 


BEHAVIOR OF FISSION PRODUCTS 


In an ordinary reactor, most of the fission products 
will remain within the fuel can, but the mobility that is 
inherent in a fluid-fuel reactor makes it necessary to 
trace the paths that the individual fission products 
take. For a molten-salt fuel, there are three broad 
classes of fission products: those whose fluorides are 
stable in the salt at its redox potential (the soluble 
fission products), the noble gases, and the noble 
metals. A great deal of practical information on the 
disposition of these different fission-product groups 
was provided by operation of the MSRE. 

The major soluble fission products are the rare 
earths (including yttrium), zirconium, barium and 
strontium, rubidium and cesium, and iodine and 
bromine. Metal coupons exposed in the MSRE core 
showed only very small deposits of these soluble fission 
products, generally corresponding to the amounts that 
might be expected from fission-product recoil or as a 
residue from deposited precursors.°' Slightly more was 
found on the graphite, presumably reaching there via 
motion of aerosols into the graphite pores. Similarly, 
very small amounts of these soluble fission products 
were found in the off-gas system, presumably carried 
there by an aerosol mist and by gaseous precursors. 

Xenon and krypton gases constitute about 27% of 
the stable fission products, but other stable or long- 
lived fission products spend a transient period as noble 
gases. The noble gases are very sparingly soluble in 
molten salt and take the first opportunity to migrate to 
any gas phase in contact with the fuel. The competitive 
sinks for noble gases are the pore spaces in the 
moderator graphite and the circulating bubbles of the 
cover gas. The reactor designer tries to encourage 
migration to the gas bubbles so that the xenon can be 
stripped from the salt to reduce the '** Xe neutron 
poison. The assumption has generally been made that it 
will be possible, by sealing the pores of the graphite 
and by efficient gas stripping, to keep the loss of 
neutrons to '*5 Xe at about one-tenth of the value for 
a solid-fuel reactor. 

In the MSRE the graphite surface was not sealed, 
and only 4% of the fuel flow was sent through a degas- 
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sing spray ring in the pump bowl. The volume frac- 
tion of bubbles in the circulating salt varied from less 
than 0.05 to about 0.7%. The poison fraction for! #5 Xe 
varied from 0.001 to 0.009, depending on the bubble 
fraction and on whether helium or argon was used as 
the stripping gas. For most of the MSRE operation 
(using helium), the '*5 Xe poison fraction was about 
0.003, as compared to an expected value of 0.015 had 
there been no xenon stripping. The variations in xenon 
poison fraction with volume fraction of the circulating 
cover gas and with the nature of the gas have been 
explained adequately by a model developed by 
Engel®? that takes into account the solubility of the 
cover gas as a function of temperature. 

Fission products that spend seconds to hours as 
noble gases and then decay to other elements can 
migrate into the graphite pores as a gas and deposit 
there as a solid. Concentration profiles as a function of 
depth below the graphite surface for such fission 
products are consistent with what would be expected 
from gaseous diffusion during their half-life as a gas.°! 
In similar manner, some of these solid fission products 
were found in the off-gas system after their transitory 
existence as noble gases. 

Another major group of fission products, consisting 
largely of Nb, Mo, Tc, Ru, Rh, and Te, does not form 
fluorides that are stable at the redox potential of the 
fuel salt and is therefore called ‘“‘noble.”’ Because they 
are born within the salt, these elements must act 
somewhat like the noble gases since they are not 
appreciably soluble, are not wet by the salt, and tend 
to migrate to the salt surfaces. The deposits of these 
noble metals on metal and graphite coupons exposed in 
the core of the MSRE were, with the exception of 
niobium, which forms a stable carbide at reactor 
temperatures, heavier on the metal surfaces than on 
graphite surfaces.°? By a postmortem gamma scan of 
the heat exchanger, Houtzeel, Blumberg, and Dyer? ’°* 
found a heavier deposit of noble metals than would 
have been expected from the samples exposed in the 
reactor core. The deposits were further enhanced by as 
much as a factor of 3 in the regions of high turbulence 
caused by the heat-exchanger baffles. 

Noble metals were also found in the pump bowl of 
the MSRE, both asa film on the surface of the salt and 
in the gas space. Evidence for the surface film was 
provided by the fact that sampling capsules always 
emerged from the salt with a deposit of noble metals 
on the outside surface. The gas samples that were 
collected from the space above the liquid level also 
showed the presence of noble metals. Blankenship has 
proposed that the ejection of noble-metal atoms from 
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the surface of the molten salt should be expected as a 
result of the conversion of the surface energy of a 
nonwetting species to kinetic energy.°5 

By observing the relative concentrations of an 
isotope pair (1°? Ru/'°°® Ru) having different half-lives, 
Bohlman and Compere®? were able to measure the age 
of the noble-metal fission-product deposits. Their 
results are consistent with a model that involves a 
“storage pool” from which the noble metals are 
redistributed after a time. The metals leave the salt 
soon after their birth; some are deposited immediately 
on metal and graphite surfaces, but some are delivered 
to the storage pool, presumably in the pump bowl, 
where they are held up for periods of from 45 to 90 
days. After this period, some of the metals from the 
pool are distributed to the graphite and metal surfaces, 
and some go into the off-gas system. 

In a power reactor, it is important to know the 
distribution of these noble-metal fission products so 
that the afterheat generation in various reactor compo- 
nents can be calculated. Briggs®® has estimated this 
distribution using a model based on mass-transfer 
theory. The distribution between the surfaces of 
graphite, metal, and bubbles is then taken as propor- 
tional to the products of the mass-transfer coefficients 
and the surface areas. In this way, it was estimated 
that, in a power reactor, 10% would deposit on the 
graphite, 40% on the metal surfaces, and 50% would go 
into the off-gas system with the bubbles. The model 
gives results reasonably consistent with the data ob- 
tained from the MSRE. 


COOLANT SALTS AND TRITIUM 
CONTROL 


In all serious studies of molten-salt reactors, an 
intermediate coolant has been proposed to carry heat 
from the fuel salt to the steam generators. A steam leak 
directly into the fuel salt would cause rapid precipita- 
tion of oxides, including UO, , and make the fuel salt 
corrosive by the formation of HF. Since such a leak is a 
likely event and the consequences are serious, a 
protective intermediate coolant loop seems necessary. 
Many coolants have been repeatedly considered for this 
role, but only lead and various salt compositions have 
received sustained attention. 

The MSRE used a salt with a molar composition of 
66% 7 LiF—34% BeF, as the coolant, a choice made for 
maximum safety. With such a salt, leakage into the fuel 
salt would be completely harmless and, except for a 
potential health hazard, leakage into the radiator air 
would be uneventful. For power reactors, its disadvan- 
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tages of high cost and relatively high melting point 
(455°C, 850°F) seem to outweigh its advantages. 
Although other salts are occasionally considered, the 
current choice is a eutectic mixture of sodium fluoride 
and sodium fluoroborate with a molar composition of 
8% NaF—92% NaBF,. It is inexpensive, has reasonable 
heat-transport properties, and has a liquidus tempera- 
ture of 725°F. Its vapor pressure is relatively high,°” 
with about ¥; atm pressure of BF3 gas at the expected 
maximum operating temperature. As a result, the cover 
gas must contain an equilibrium partial pressure of BF; 
to avoid a gradual change in composition in the 
direction of NaF. After some slight initial difficulty, 
this has been accomplished on a routine basis in 
circulating loops.°® 

The sodium fluoroborate salt is considerably more 
corrosive than the fuel salt. Corrosion losses are quite 
sensitive to impurities, particularly water vapor, and 
corrosion rates of several mils/year have been observed 
for short times. With reasonable impurity levels, an 
average corrosion rate of 0.7 mil/year has been 
measured. One loop,°®? operated with high-purity salt 
and at rather extreme temperatures (hot leg 687°C and 
cold leg 438°C), has given the best performance with a 
corrosion rate of 0.2 mil/year. During periods of high 
corrosion, caused by high impurity levels, there has 
been a tendency for deposits, which seem to consist 
largely of Na3CrF,, to occur in the cold leg, leading to 
the possibility of plugging. It may be necessary to use 
cold traps to alleviate this situation. Evidently much 
more good experience is needed in the operation of 
loops in which adequate purity of the salt is main- 
tained before one can be assured of obtaining a 
satisfactory life for heat exchangers and steam genera- 
tors. 

When the sodium fluoroborate salt is mixed with 
the fuel salt, it reacts to produce two immiscible 
liquids accompanied by an increase in BF; pressure.’° 
For a small heat-exchanger leak, the excess pressure 
could be handled without difficulty by the fuel off-gas 
system, but larger leaks can cause sudden generation of 
large quantities of gas, increasing the pressure by an 
order of magnitude. A more thorough investigation of 
the mixing kinetics is in order. In similar vein, the 
planned experiments related to steam leaks into the 
coolant salt would seem quite important. 

When public attention began to be focused on the 
environmental effects of the radioactive effluents from 
power reactors, serious calculations of the disposition 
of tritium from molten-salt reactors were begun. A 
1000-MW(e) molten-salt reactor will produce tritium at 
a rate of about 2420 Ci/day, with more than 98% 
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coming from (”,q) reactions with lithium isotopes. At 
the temperature of operation of a molten-salt reactor, 
iron- and nickel-based alloys are nearly transparent to 
hydrogen, and it was feared that the tritium would 
flow through the coolant salt to the steam in the steam 
generators. 

Briggs’’ has studied the problem extensively and 
devised a model to describe the transport of tritium. 
He assumed that tritium would appear in the fuel salt 
either in molecular hydrogen form (T,, HT) or as 
tritium fluoride, with the relative proportions deter- 
mined by the redox potential of the salt. These gases 
would migrate to surfaces: graphite, metal, or gas. At 
the metal surfaces, some of the molecular tritium 
would dissociate into atoms, dissolve in the metal, and 
diffuse to the outer surface. Tritium fluoride would 
release some tritium atoms by reaction with the metal, 
but this is thought to be a comparatively slow reaction. 
The T,, TH, or TF that did not diffuse through the 
walls or absorb on the graphite would find its way into 
the fuel off-gas system. 

With a reasonable choice of parameters, the model 
gave a good fit to the observed distribution of tritium 
in the MSRE. The measurements of the tritium 
distribution in the MSRE were somewhat erratic, but 
between 46 and 56% was found in the fuel off-gas 
system, 15% on the graphite, between 5 and 9% in the 
reactor cell, and a similar amount in the cooling air 
coming from the radiator heat dump. As much as 20% 
of the calculated production is unaccounted for, but 
there is good reason to believe that most of this was 
retained in a hydrocarbon residue in the pump bowl or 
off-gas system. (The hydrocarbons came from a small 
oil leak into the fuel pump.) 

Applied to a large power reactor, this model 
indicated that as much as 60% or 1500 Ci/day could 
end up in the steam system of the turbogenerator 
plant. Handling this amount of tritium in the steam 
system would be awkward, although not impossible if 
one were allowed to approach the stated maximum 
permissible concentration (MPC) values. This does not 
seem desirable, however, and so a tentative goal of the 
development program is to reduce the plant effluent to 
the level proposed as a design objective for light-water 
reactor plants, about 10 Ci/day. A number of possible 
approaches are under investigation, but only two, both 
involving trapping tritium in a coolant salt, will be 
mentioned here. 

In relatively pure fluoroborate salt, OH ions are not 
stable; when small amounts of NaOH or NaBF3;OH are 
added, the OH ions disappear by thermal decomposi- 
tion or by hydrogen-producing reactions with the 
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Hastelloy N container.’* However, there is evidence ° 
that in the presence of higher concentrations of oxide 
and Cr** ions, steady-state concentrations of hydro- 
gen as high as 30 ppM can’be maintained as hydroxyl in 
fluoroborate test loops. Exchange of hydrogen isotopes 
in such hydroxylated fluoroborate melts has been 
observed, although the rates of exchange have not been 
measured. It seems likely that tritium entering the 
coolant salt will exchange with the hydrogen already 
there and thereby be retarded from passage into the 
steam system. Such a holdup mechanism combined 
with a side-stream processing scheme which removed 
the hydroxyl ion would permit the safe storage of the 
tritium. The tritium so obtained would be diluted with 
hydrogen, but, if such dilution were less than a factor 
of 10%, the volume could be handled readily. 

The other method of preventing the tritium from 
getting into the steam system which has received 
considerable attention is the use of a secondary coolant 
salt, a mixture of NaNO3, NaNO,, and KNO;3 known 
as HTS or Hitec, to isolate the steam generator from 
the primary-coolant salt. If the tritium were oxidized 
to water, the Hitec would effectively trap tritium, but 
its use has the obvious disadvantage of adding a third 
loop and increasing the capital cost. 

The Hitec salt as a replacement for the primary- 
coolant salt has not been given serious consideration 
for several reasons. The mixture may not be suffi- 
ciently stable at the higher temperatures it would 
experience in the primary heat exchanger, and exten- 
sive testing would be necessary to ensure its stability to 
the radiation field in the heat exchanger. A leak of 
Hitec into the fuel would cause precipitation of UO, 
and a large leak would probably result in rapid reaction 
with the moderator graphite. 

Current development work relating to the tritium 
problem includes measurement of the solubilities and 
diffusion constants needed in the model describing 
tritium behavior. The hydrogen content of the sodium 
fluoroborate salt is still under investigation, as is the 
means for recovery of the hydrogen and tritium from 
it. The high-temperature stability of the Hitec salt is 
also under investigation, should there be a necessity for 
using it. 


RADIATION BEHAVIOR 
OF HASTELLOY N AND GRAPHITE 


Hastelloy N was developed specifically as the 
container material for use with molten fluorides. It was 
tailored to have good high-temperature strength, good 
air-oxidation resistance, good fabrication properties, 
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and to be resistant to corrosion by the fuel salt. 
However, Hastelloy N has been found to suffer from 
the high-temperature radiation embrittlement common 
to most iron- and nickel-based alloys. This embrittle- 
ment results from internal generation of helium either 
by thermal-neutron absorption in boron impurities or 
by fast-neutron reactions with nickel. At high tempera- 
tures the helium gathers at the grain boundaries and 
greatly reduces the tolerance for strain before fracture. 
Under some conditions the strain before rupture is less 
than 0.5%. 

A development effort to improve the high-tempera- 
ture radiation performance of the alloy has been 
under way for some years.” The general approach has 
been to produce a fine carbide precipitate within the 
grains and at grain boundaries. By this means the 
helium is inhibited from migrating to the grain bound- 
aries, and the grain boundaries are strengthened. The 
first modification consisted in reducing the molybde- 
num from 16 to 12%, reducing the silicon from 0.5 to 
0.1%, and adding 0.5% titanium. This alloy was 
satisfactory at temperatures up to 650°C, but its 
properties deteriorated above that temperature. The 
cause of the deterioration was shown to be the 
conversion of the finely dispersed MC-type carbide to a 
much coarser carbide of the type MC,. (The metal 
“M” in these carbides is a mixture, high in molybde- 
num.) Recent development work has shown that 
improved ductility and strength at higher temperatures 
can be obtained by somewhat larger additions of 
various metals, such as Ti, Nb, Zr, and Hf. The niobium 
seems to provide improvement only in combination 
with the others, and zirconium makes the alloy 
difficult to weld satisfactorily. In laboratory-prepared 
samples, the hafnium greatly improved the ductility, 
but later commercially produced heats did not perform 
so well. At present the best solution seems to be an 
alloy containing about 1.5 to 2.0% titanium, which 
gives promise of having a ductility of 4% or more when 
irradiated at temperatures up to 760°C and of having 
other properties that are satisfactory. There is some 
concern that the phenomenon of conversion of 
MC-type carbides to MC,-type carbides provides a new 
upper limit to the temperature excursions allowable in 
a molten-salt reactor. It seems probable that short-time 
overheating beyond the temperature for which the 
alloy is tailored will allow the formation of coarse 
carbide precipitates, thus rendering the irradiated alloy 
brittle and endangering future reactor operation. If this 
proves true, it will provide an additional engineering 
constraint in design, but one that can probably be met 
without too much difficulty. 
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The appearance of fissures between the surface 
grains was a new problem encountered with the 
Hastelloy N alloy during the operation of the MSRE.”* 
Fissures were first observed in core surveillance speci- 
mens on which tensile tests were run. Although no 
decrease in tensile strength was observed, photomicro- 
graphic examination of the specimens after being 
elongated by the test revealed fissures to a depth of 
about 5 mils, the thickness of the surface grains. 
Extensive examination of many specimens after the 
MSRE was shut down showed that these cracks 
occurred wherever the alloy was exposed to fuel salt, 
but almost never anywhere else. For example, the 
surface of the heat-exchanger tubes exposed to the fuel 
salt showed fissures, but the opposite side of the tubes, 
exposed to the coolant salt, showed no attack. Simi- 
larly the surface of the control-rod thimbles exposed to 
fuel salt was attacked, but not the other side which saw 
a gaseous atmosphere. 

Two possible explanations are being explored. One 
is that the fuel was quite oxidizing during part of the 
MSRE operation, and there might have been excessive 
grain-boundary removal of chromium during this 
period. The other explanation, which now appears 
more likely, is that some fission product, probably 
tellurium, is responsible, probably by forming an 
intermetallic compound. This subject is under intensive 
investigation at present. Surface grain-boundary defects 
have been produced artificially, both by exposure to a 
highly oxidizing salt and by heat treating with tellu- 
rium introduced as a surface impurity. The effects with 
tellurium seem to resemble those seen in the reactor, 
leading to the belief that tellurium is the primary 
responsible agent. Until further information on the rate 
and extent of the attack and its mechanism is known, 
it is difficult to assess whether the matter is serious or 
not. 

In a power reactor the displacement damage to 
graphite produced by fast neutrons will create some 
problems. At the high temperature of the molten-salt 
reactor, the principal harmful effect will be the gross 
dimensional changes. The usual graphites at first shrink 
as a result of a contraction of the hexagonal graphite 
crystals in their “a” direction. After receiving a fluence 
of about 2 x 10?? neutrons/cm? with energy greater 
than 50 kV, the graphites begin to expand, reaching 
their original dimensions at about 3 x 107? neutrons/ 
cm?. After this, they expand at an ever-increasing rate, 
presumably to destruction. 

Most graphites are highly anisotropic, and thus the 
dimensional changes are quite different along different 
axes, but this effect can be kept small by choice of raw 
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materials, primarily by the use of rounded particles 
instead of needle-shaped ones. With this “isotropic” 
graphite,’* the linear shrinkage before the expansion 
starts may be as much as 2.5%. As the temperature of 
irradiation increases, the maximum shrinkage de- 
creases, but the rapid expansion phase occurs at a 
lower fluence. 


The graphite in a reactor will see a nonuniform 
fluence at nonuniform temperatures, so that a careful 
stress analysis is required.© Fortunately the graphite 
also undergoes radiation creep, and this helps to keep 
the stress problem from being serious. Similarly the 
changes in reactor physics resulting from dimensional 
changes have been investigated and found to be 
acceptable. 


How long the graphite moderator can remain in the 
reactor core before it is replaced is a matter of 
judgment. As the expansion from radiation damage 
progresses, the graphite will become more porous, and 
the breeding performance will gradually show a slight 
deterioration due to the absorption of fission products 
and to the displacement of some fuel from the reactor 
core. At some point the expansion could cause binding 
and possible breakage of the graphite. For practical 
purposes, it has been assumed that the graphite life in 
the reactor will be determined by the time it takes for 
the graphite to return to its original volume after the 
shrinkage phase has been completed. With this restric- 
tion the graphite life is about 4 years if the maximum 
power density in the reactor core is 70 kW/liter. The 
replacement of the graphite core every 4 years is 
expected to cost'® about 0.2 mill/kWh. 


Graphites that do not follow the above-described 
pattern of shrinkage and expansion under radiation 
have been made. This gives hope that some improve- 
ment in graphite longevity can be obtained. Develop- 
ment work is currently directed at understanding the 
deviations that have been observed in radiation be- 
havior, of course, with the hope that some improve- 
ment can be made. 


FUEL PROCESSING 


A number of studies have been made of the 
dependence of breeding performance and fuel-cycle 
costs on the cycle times for removing fission products 
from the fuel of a molten-salt reactor.2*75 77 The 
breeding performance, of course, improves with rapid 
fuel processing, but the net fuel-cycle cost shows 
remarkably little variation between the extremes of 
virtually no processing and the rather sophisticated 
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processing required to get the optimum breeding 
characteristics. With little processing the reactor has a 
conversion ratio between 0.8 and 0.95, and, when 
operated as a breeder, the breeding ratio is about 1.07. 
With a specific power of about 1.5 MW(t)/kg, the 
breeder has a fissile production credit of about 0.1 
mill/kWh. As far as fuel-cycle cost is concerned, what 
the converter loses in fuel production credit and fuel 
purchase cost, it makes up for in decreased cost of the 
processing plant and its operation. In both cases the 
fuel-cycle cost has been calculated to lie between 0.6 
and 0.8 mill/kWh, although recent capital-cost esti- 
mates of the processing plant’® will lead to higher 
values. With the wide range of processing speeds 
economically interesting, a number of different com- 
binations of individual processing steps have been 
explored. 

All of the reactor performance calculations assume 
that '° Xe is removed on rapid enough cycle to keep 
the xenon poisoning to 0.005 neutron per neutron 
absorbed in fissile atoms. The problem here is to 
prevent the xenon from entering the porous graphite 
moderator. In the past it has been assumed that this 
could be accomplished by a rapid purge of the fuel 
with helium or other gas, combined with some method 
of sealing the surface pores of the graphite. A 
remarkably good method of sealing the graphite pores 
was developed,* involving heating the graphite to 
about 900°C and alternating between evacuation and 
admission of a hydrocarbon gas such as butadiene. 
Unfortunately the low permeability was not retained 
after irradiation to a fast-neutron fluence of about 
10?? neutrons/cm? . Development of sealing techniques 
is continuing, both with the “pulse-impregnation” 
technique described above and with isotropic pyrolitic 
coatings put on at a somewhat higher temperature. 
Occasional samples have retained their low permeabil- 
ity to fluences of over 107? neutrons/cm?, but it is by 
no means certain that large-scale success will be 
achieved. 

Additional attention is being given to the possibil- 
ity of removing part of the iodine precursor of xenon 
by giving a mild HF treatment to the side stream of 
fuel going to the off-gas system.’? If this proves 
practical, the xenon poisoning can be kept within the 
specified 0.005 limit without any special sealing of the 
graphite pores. If neither pore sealing nor iodine 
removal becomes practical, the xenon poisoning can be 
kept to between 0.01 and 0.02 by the planned 
gas-purging system, and the loss from design per- 
formance will be about 0.01 in breeding or conversion 
ratio and less than 0.02 mill/kWh in fuel-cycle cost. 
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Protactinium is a fairly serious neutron poison for 
the breeder since, for each neutron absorbed, a 
potential ?*°U atom is also destroyed. At a specific 
power of 1.5 kW(t)/g, the present reference-design 
reactor would have a breeding-ratio loss of about 0.05 
if the protactinitum were allowed to decay in the 
reactor fuel. It is therefore desirable to remove the 
protactinium from the fuel soon after it is formed and 
to store it somewhere outside the reactor until it 
decays to 733U. Removing the protactinium on a 
10-day cycle® reduces the breeding-ratio loss to about 
0.01. 

It would be most desirable to remove the protac- 
tinium without also having to remove the uranium, and 
this may be possible through an oxidation process. 
Early work®® showed that the addition of excess 
amounts of BeO, CaO, or ThO, would precipitate 
protactinium from fluoride melts and that the protac- 
tinium could be brought back into solution by HF 
treatment. More recently, it has been shown*>">* that, 
if protactinium is first oxidized to the 5+ valence, it 
can be precipitated as Pa.O,; with a clean separation 
from uranium. Bell and McNeese** have described a 
possible flow sheet for a processing plant using this 
process. 

Another way to remove the protactinium, cur- 
rently the favorite, is by reductive extraction from the 
fuel salt into bismuth in a multistage countercurrent 
extraction column. The reductant is lithium dissolved 
in the bismuth, in such amount that the thorium it 
extracts will not exceed the thorium solubility in 
bismuth. Many experiments have shown that the 
distribution coefficients for uranium, protactinium, 
and thorium are widely different, allowing for clean 
separations. Since uranium is preferentially extracted, 
it must be removed from the salt first, either by prior 
reductive extraction or by fluorination. 

Once the protactinium is extracted into the bis- 
muth, it is converted back into a fluoride by treatment 
with HF, and the protactinium fluoride is transferred 
into a separate recirculating salt stream. This salt 
stream is then fluorinated to separate and recover the 
2330 as it is produced from the protactinium. The 
reductive extraction process has now been demon- 
strated successfully in a flowing system, transferring 
uranium from salt to bismuth in a _ packed 
column.®!°8? 

The reductive extraction process requires a con- 
struction material that is compatible both with fluoride 
salt and with bismuth containing dissolved lithium. 
Corrosion tests have shown that molybdenum and 
graphite should be suitable from this point of view,” 
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but both materials are brittle and present problems in 
fabrication. A considerable development effort has 
been expended on the forming and joining of molybde- 
num. It has been found that short sections of pipe up 
to 12 cm in diameter can be made by back extrusion, 
as can suitable end caps with the desired protuber- 
ances. A series of Fe—Mo—C brazing alloys for joining 
molybdenum has been developed which may be usable 
with bismuth, but their suitability for long-term use 
has not yet been demonstrated. Since welds in molyb- 
denum are brittle at room temperature, all welded 
joints will be either back brazed or strengthened by 
brazed-on sleeves. Tight joints have been made both 
with a gasketed mechanical coupling that can be 
disassembled and by roll bonding.®? The progress that 
has been made in the technology of fabricating 
equipment from molybdenum is impressive. Plans are 
also under way to make increasing use of graphite, 
especially for such equipment as mixer—settler tanks 
that can be fabricated by simple machining operations. 

As noted above, the reductive extraction method 
for isolating the protactinium requires the prior re- 
moval of uranium. The complete removal of the 
uranium from the fuel salt on as short a cycle as 10 
days is an ambitious undertaking, requiring a method 
that offers little chance for loss of uranium before it is 
reconstituted in the fuel. Several methods that appear 
to meet this requirement are available. One is reductive 
extraction into bismuth, and this method was consid- 
ered along with an electrolytic method of putting the 
uranium back into the fuel salt after the protactinium 
had been removed.8* A number of factors, including 
the amount of reductant required and the difficulty of 
developing a continuous electrolytic process, caused 
this method to be abandoned in favor of fluoride 
volatility. 

Fluoride volatility is a proven method for the 
recovery of uranium from molten fluoride salts. This 
process involves sparging the molten salt with a 
mixture of fluorine and helium, converting the UF, to 
gaseous UF,. As a batch process, this method was used 
to recover the *7°U from the first charge of fuel salt 
that was used in the MSRE. With some modifications 
the equipment used for the MSRE could be used for 
the infrequent processing of a converter reactor. For 
the larger-scale uranium removal needed by the breeder 
reactor, a continuous fluorinator is considered to be 
necessary . 

The most troublesome point of such a device is 
corrosion at the liquid—gas interface. Where the metal 
surface is exposed only to gas, a protective fluoride 
film is formed, but this film is dissolved whenever it is 
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contacted by the molten salt. To overcome this 
problem, we plan to operate the fluorinator with a 
layer of frozen salt lining the vessel walls. This requires 
the combination of an internal heat source to keep the 
salt moltei: and external cooling to provide the frozen 
layer on the walls. When the reactor is operating with 
irradiated fuel, the radioactivity will provide the 
internal heat, but, in the development of the equip- 
ment, we plan to use induction or resistance heating of 
the salt. The first such experiment with molten salts 
should be under way soon.®* 

After the xenon poisoning has been limited to a 
low value, more than 85% of the neutron losses to the 
remaining fission products are due to the rare earths. If 
the rare earths are removed continuously from a small 
stream of fuel with equivalent cycle times of 20 to 50 
days, the total fission-product poisoning can be kept to 
less than 0.02 neutron per neutron absorbed in fissile 
atoms, a very satisfactory level for a breeder reactor. 
The preferred method for removal of the rare earths is 
the so-called “metal-transfer process.”°° In this pro- 
cess the rare earths are first extracted from a fuel-salt 
stream into lithiated bismuth in the same manner as 
the protactinium, after all the uranium and protac- 
tinium have been previously removed. An additional 
step is necessary because it is difficult to separate the 
rare earths from the thorium cleanly by this first 
extraction. The second extraction step is from the 
bismuth to a stream of LiCl, which discriminates 
strongly against thorium. The concentration of rare 
earths in the LiCl is kept low by a further extraction 
into a very strongly reducing stream of molten Li—Bi 
alloy. An extra bonus of the metal-transfer process is 
that barium and strontium are also separated. 

The process has been demonstrated in a small 
engineering experiment®’ in which lanthanum and 
neodymium were transferred from a fuel salt to the 
final Li—Bi alloy and concentrated relative to thorium 
by a factor of 10°. 

If the oxidation process has been used to remove 
the protactinium, it will then be necessary to remove 
uranium as an extra step prior to the metal-transfer 
process. This can be a smaller operation than the 
uranium removal prior to the protactinium extraction 
because the cycle time is likely to be more nearly 50 
days than 10 days. The fluoride-volatility process is 
still a valid way to do this, but there is also the 
possibility of removing the uranium by precipitating it 
as the oxide. 

Taking advantage of oxide-equilibria studies by 
Bamberger and Baes,°° which showed that uranium 
might be precipitated somewhat preferentially in the 
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presence of thorium, Bell and McNeese®? have pro- . 
posed a process flow sheet. Experimentation is cur- 
rently being carried out with encouraging results. A salt 
containing both uranium and thorium is being treated 
with water vapor as an oxidant. Early results indicate 
that the initial precipitate, high in UO,, does not 
equilibrate with the lower concentration of uranium 
remaining in the solution. Taking advantage of this fact 
should make possible the precipitation of most of the 
uranium in a single stage without removing much 
thorium. The principle of uranium oxide precipitation 
seems to be well established, and, although its embodi- 
ment into an engineering process is in its infancy, rapid 
progress is expected. 

Although the metal-transfer process appears to give 
the best fuel-salt purification for a high-performance 
breeder, there are other possibilities for rare-earth 
removal that are perhaps worth keeping in mind. If the 
bismuth-contacting systems prove expensive or if un- 
seen difficulties develop, the other methods may be 
applicable, especially at some sacrifice in breeding 
performance. 

When two fluid breeders with separate salt streams 
for fuel and fertile materials were in vogue, distillation 
was the reference method for freeing the fuel salt of 
rare-earth fission products. Many experiments®® on the 
relative volatility of the different constituents indi- 
cated that, by continuous operation of a still at 
1000°C, it would be possible to recover the carrier salt 
as a distillate with its composition 66 mole % LiF and 
34 mole % BeF,. For this composition of distillate to 
appear, the composition of the liquid in the still would 
be about 88 mole % LiF and 12 mole % BeF, at the 
start of distillation, although the concentration of rare 
earths would gradually build up during the course of 
the distillation. If the still bottoms were recycled 
through a 400-ft? holdup tank, it was estimated that 
their disposal could be deferred for the life of a 
1000-MW plant.*®° 

Batchwise distillation was carried out on an engi- 
neering scale, both with nonradioactive salt and with a 
portion of the fuel salt of the MSRE. The results were 
in accord with expectations, except that there may 
have been some carry-over of nonvolatile species by 
mist in the runs on the MSRE fuel.?° In the 
nonradioactive tests the separation factor for rare 
earths from LiF was in excess of 1000, but with the 
MSRE fuel the indicated separation factor became as 
low as 100 (for cerium). 

With the shift of attention to the one-fluid breeder, 
in which the ThF, and UF, are in the same salt, 
interest in distillation diminished because the 1000°C 
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distillation would not recover the thorium, and, for a 
breeder to be effective in conserving fuel, the thorium 
must be recycled. Smith, Ferris, and Thompson?! 
indicate that distillation at 1200°C might make the 
recovery of thorium possible, but the problem of 
construction materials for this higher temperature 
discouraged further consideration at the time. The 
possible construction materials for the 1200°C still are 
molybdenum and graphite; since these are also the 
materials that must be used for the now-popular 
reductive extraction process, and the problems of 
constructing equipment from them are being faced 
anyway, it may be worth considering distillation again. 

In some recent experiments an attempt was made 
to precipitate plutonium from a molten-salt solution 
by oxidation, and it was found to be very stable. There 
was no loss of plutonium from solution®’ even in the 
presence of ThO, and BeO. From the similarity of 
plutonium chemistry to that of the rare earths, it is 
probable that they also are difficult to precipitate as 
oxides. This suggests that, after the protactinium and 
uranium have been precipitated as oxides, the thorium 
may also be separated out in this way, leaving the 
carrier salt containing the rare earths ready for distilla- 
tion. There are obvious difficulties with this suggestion, 
but it is made to indicate that there may be other 
future possibilities for processing. 

Similarly Grimes* has repeatedly suggested that the 
precipitation of mixed trifluorides, perhaps using natu- 
ral CeF, as the saturating fluoride, could be used to 
remove the rare earths. Perry? estimated that the use 
of this method for a relatively high-performance 
breeder would require about 35 tons of cerium per 
year for a 1000-MW reactor. With these and other 
possibilities for chemical processing relatively un- 
explored, there seems to be much room left for 
interesting investigation. 


DESIGN AND DEVELOPMENT OF 
ENGINEERING COMPONENTS 


Most of the recent design and development work 
on components has been directed to items that need to 
be changed significantly from those used on the MSRE. 
A prime example is the off-gas system proposed for a 
power reactor. The MSRE used a simple spray ring in 
the pump bowl to establish contact between the fuel 
salt and the purge gas, which was introduced primarily 
down the annulus around the pump shaft. The plan for 
a power reactor is to have specially designed gas-bubble 
injectors and centrifugal gas separators in a side stream 
of fuel salt. Extensive experiments have been carried 
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out with such devices, using water and water—glycerine 
mixtures. Several different bubble generators have been 
tested and appear workable, and good efficiency has 
been obtained with a gas separator.?* Tests with a 
closed aqueous loop are in progress, and a molten-salt 
test loop is under construction.?4 

The MSRE had no steam generator, so that work in 
this area is new. The preference at ORNL is for a 
supercritical steam generator, primarily because it is 
expected to lead to fewer difficulties in adapting it to 
the high melting point of the fused-salt intermediate 
coolant but also because of its high thermal efficiency. 
The design studies that have been made are based on 
modifications of the TVA Bull Run supercritical power 
station.2*> The Foster Wheeler Company, under con- 
tract to ORNL and the AEC, has started a design study 
leading to a development program. In the meantime 
the compatibility of the Hastelloy N alloy with high- 
temperature steam is being tested in a cooperative 
program? between ORNL and the TVA. These corro- 
sion tests have shown a metal removal rate of less than 
0.00063 cm/year at 538°C (1000°F). Tube burst tests 
to determine if there is a stress effect on the corrosion 
are under way. 

The sodium fluoroborate intermediate coolant salt 
has been circulated for 11,500 hr at 750 gal/min in an 
Inconel loop without excessive corrosion or visible 
attack on the pump impeller.?° The biggest difficulty 
was the nearly complete plugging of the bubbler tube 
used to indicate the liquid level in the pump bowl. A 
new engineering-scale pumped loop of similar size is 
being constructed out of Hastelloy N piping taken 
from the coolant loop of the MSRE. It will have several 
side streams that can be used for a variety of 
purposes,”” but the principal function will be to serve 
as a vehicle for testing the control of the chemistry of 
the coolant salt. High on the priority list will be trial of 
schemes for the control of tritium in the coolant salt. 

In the engineering design of molten-salt power 
plants, a great deal of attention has been given to 
safety and maintainability. The experience in main- 
taining the MSRE has been drawn upon generously. 
The philosophy of maintenance is “remove, replace, 
and repair or discard.”?* To remove and replace 
requires that all service lines be accessible for discon- 
nection and connection by long-handled tools. The 
major equipment items must be laid out so that they 
are visible from the maintenance area, and their 
positions must be accurately fixed so that replacement 
equipment will fit in place. Flanges were used to 
connect the major pieces of equipment in the MSRE, 
but the larger-sized piping for the power reactors will 
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require welding in preference to flanges. A program to 
develop a machine for remotely cutting, facing, posi- 
tioning, and welding pipe has been under way for some 
time, and the machine appears quite promising for 
molten-salt reactor use.? 

Although there is optimism about maintenance 
based on the experience with small pieces of equip- 
ment in the MSRE, it must be recognized that 
replacement of major reactor cell equipment will be 
very difficult and probably time consuming, at least 
until some experience is gained. Consequently even 
greater than usual emphasis must be placed on reliabil- 
ity of equipment. 

Another problem that has received some attention 
is that of replacing the graphite core when its life has 
ended. One plan is to lift the whole core out as a unit 
in a cage,'” along with the roof of the reactor vessel. 
The steps by which this might be done have been 
worked out and appear to be feasible. The magnitude 
of the task is comparable to the proposal in the 
“hot-cell” maintenance concept of the LMFBR to lift 
the reactor head with its shielding as a unit. Other 
designers may prefer to group the graphite sticks in 
bundles and remove the bundles one at a time. 

Particular attention has been given to the design of 
the drain tank and its method of cooling.’? The great 
desire is to provide a safe “‘walk-away”’ situation in case 
of a severe accident, with the fuel automatically 
drained into a tank that is always cooled by natural 
convection and requires no engineered safeguards. 
Several plausible methods have been proposed. 

The safety engineering has naturally been con- 
cerned with accident situations, mostly those in which 
circulation of the fuel salt or the coolant salt has been 
lost. In such cases the afterheat produced by deposited 
noble-metal fission products must be taken into ac- 
count. For example, studies have been made of the 
peak temperature that would be reached by the 
intermediate heat exchanger if it were drained of both 
fuel and coolant salts. Here it would be desired to use a 
tube bundle enclosed in an inner shell so that it could 
be replaced as a unit, but the extra heat-transfer barrier 
presented by the inner shell makes its use prob- 
lematical.!°° 

Some of the problems of safety of a molten-salt 
reactor will be different from those of other reactors. 
No accidents comparable in chain of events to the 
loss-of-coolant accident in light-water reactors or the 
pump-coastdown-without-scram accident in LMFBRs 
have yet been postulated for molten-salt reactors, so 
that there is reason for somewhat less concern with big 
accidents. On the other hand, a large MSBR plant will 
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have its radioactivity normally spread over a larger area 
and in more pieces of equipment, and great care must 
be taken to avoid containment leaks. This problem can 
be faced squarely only when a serious design is 
undertaken for a reactor that will be constructed. 


SUMMARY 


Any appraisal of the current status of molten-salt 
reactor development or any proposal as to what should 
be done next must be a personal one reflecting the 
experience and interests of the writer. However, three 
separate groups of engineers who have studied the 
molten-salt system agree that the present technology is 
adequate to justify starting the design and construction 
of either a power-producing prototype or a reactor 
experiment of significant size and performance, such as 
the proposed MSBE, or both. For the smooth progres- 
sion to a molten-salt reactor power-producing industry, 
these are the most urgently needed steps, for they will 
lead to the development of full-size equipment, to the 
realistic facing of difficult design decisions, and to the 
all-important operating experience. 

That there is adequate information for a decision 
to proceed with plant design and construction projects 
does not mean that the development is completed; all 
such proposals allow a long time, such as 7 or 8 years, 
for the expected completion of such a plant. During 
that time the larger pieces of equipment would be 
designed and tested, and a number of opportunities for 
altering the details of design would be available. New 
knowledge from a concurrent research and develop- 
ment program could be utilized. 

I will try to summarize the status of the develop- 
ment program in the approximate order of the pre- 
ceding technical sections. 

The understanding of the basic chemistry of the 
fuel salt seems to be adequate at present for any 
normal, and for some abnormal, operating conditions 
of a reactor. Further work on the effects of an 
increased oxidation potential (to very low ratios of 
UF3/UF,) on the compatibility limits of the fuel 
should be valuable since some of the fuel-processing 
steps may leave the fuel salt quite oxidizing. The 
compatibility with both Hastelloy N and graphite 
should be explored under these oxidizing conditions, at 
both normal and high temperatures. 

The voltammetric method of measuring the redox 
potential of the fuel is a very important development. 
Some of the chemistry research should be devoted 
toward improved or alternate methods of determining 
the chemical condition of a fuel salt under the highly 
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radioactive conditions that exist in a reactor plant. 
Along the same line, more work seems to be justified 
on methods for adjusting the redox potential of a 
reactor, including the kinetics of the reactions. 

The behavior of the gaseous and the soluble fission 
products seems to be well enough understood. For the 
noble metals the situation is not as good, and much 
more experimentation is required. Such experimenta- 
tion is difficult because it is not easy to simulate the 
physical situation that a newborn fission-product atom 
sees. Subjects needing attention are agglomeration, 
adhesion to surfaces, transport in purge gas—in fact, 
any knowledge that bears on the mode of transport 
and the place and permanence of deposition. Also, the 
evidence in the MSRE of different conditions of 
surface tension of the fuel salt, which may have been 
associated with fission products, needs attention. Reac- 
tions of the various fission products with the graphite 
should be studied to learn the conditions of tempera- 
ture and redox potential that will cause carbides to be 
formed. 

Although sodium fluoroborate will probably turn 
out to be satisfactory as the coolant salt, much more 
needs to be learned about it before one can predict a 
30-year life for the heat exchanger or the steam 
generator. Basic chemistry studies related to its corro- 
sive tendencies are needed, leading to knowledge of 
purity requirements and methods of control. The 
events that occur when it is mixed with fuel salt also 
need further elucidation. The sodium fluoroborate salt 
was chosen because of its low cost, low melting point, 
and because one could recover readily from some 
contamination of the fuel salt. Recently the fact that 
its small hydrogen content may provide a means for 
control of tritium has counted in its favor. Other salts 
have been proposed, and I believe it is worthwhile to 
keep up some activity exploring the other possibilities. 
It may be possible to find a different set of compro- 
mises that are as acceptable as those associated with 
the sodium fluoroborate salt. 

Research on the habits and control of tritium is 
already intense and should be continued until its 
behavior is understood and a method of control is 
assured. Experience with the MSRE left some ques- 
tions about absorption of tritium in graphite and about 
a reaction of tritium with hydrocarbons that should be 
explored. In my opinion the tritium can be controlled 
satisfactorily; the question is just how and at what 
cost. 

Improvement in the resistance of Hastelloy N to 
helium embrittlement is well along, and this work 
should be completed without too much difficulty. 
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More needs to be learned about the effects of 
inadvertent overheating of irradiated Hastelloy N. The 
investigation of its suitability for steam-generator use 
needs to be continued, and, if it should turn out not to 
be suitable, some modification or substitute must be 
found. 


The investigation of the surface fissuring of 
Hastelloy N in the MSRE fuel salt is in too early a stage 
to estimate the outcome. It is encouraging that 
out-of-pile experimentation methods are showing the 
effect, so that much more can be learned quickly about 
the reactions involved, their dependence on fission- 
product concentration, their kinetics, and the ultimate 
extent of the attack. Some consideration had already 
been given to using a cladding metal for the primary 
system, and work along this line should probably be 
encouraged as a possible way to cope with the surface 
fissuring. 

Graphites suitable for molten-salt reactors can be 
manufactured with present technology, but their 
limited life in the fast-neutron environment adds 
inconvenience and cost. The fact that different graph- 
ites show different radiation behavior provides hope 
that some improvement in radiation resistance can be 
achieved. The subject should be explored until the 
causes of the observed differences are better under- 
stood because the subject is of long-term importance. 
However, one should not expect more than modest 
improvements over the performance of presently 
known types of graphite. 


It has not been many years since the time when 
molten-salt reactor engineers could only hope that 
suitable chemical processing methods would be devel- 
oped, basing that hope on the expectation that the 
liquid state of the fuel would make chemical processing 
easier and on faith in the ingenuity of the chemists. 
Now there seems to be a relative abundance of possible 
methods, at least two for each needed step. At present, 
all that can be advised sensibly is to continue the 
development of the reductive extraction and the 
oxidation methods for removal of protactinium, the 
metal-transfer process for the periodic removal of rare 
earths, and the oxidation and fluoride-volatility meth- 
ods for extraction of uranium. Other processes, such as 
distillation, should not be completely forgotten. It is 
important to keep several options open because the 
relative cost of the different chemical processing steps 
will increase in importance as commercial reactor 
operation nears. 


The responsibility for keeping a low xenon level in 
the reactor lies with the engineers, and current experi- 
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ments on the introduction and stripping of bubbles 
should do the major job. The work on sealing the 
graphite pores should not be abandoned until the 
utility of depositing thinner layers of pyrolysis prod- 
ucts is fully explored. Lighter deposits may retain their 
integrity under radiation better. In this connection it 
should be remembered that the gaseous permeability 
measurement used as a measure of the degree of 
surface sealing may yield a pessimistic estimate of the 
effect of the salt—graphite interface in the prevention 
of xenon entry into the graphite. 

A great deal more experience is needed in the 
handling of the coolant salt, including control of its 
chemistry during operation. Any corrosion products 
must be kept from plugging small passages. One of the 
experiments with the coolant salt should involve a 
small steam generator so that any compatibility prob- 
lems wiil be encountered early. The full development 
of a steam generator is of great importance; experience 
in this area is perhaps the most glaring omission to 
date. 

The MSRE operated well and provided much useful 
information on the reliability of its equipment. Its size 
was small, for its main purpose was to demonstrate 
that a molten-salt reactor could be run successfully. It 
accomplished its purpose admirably, and, in the normal 
course of events, the next step would be to proceed to 
the construction of a larger reactor that could serve as 
the proving ground for larger and higher power-density 
equipment. The construction of a larger reactor costs 
much more money and raises the question, “Is it worth 
it?” 

As far as can ve determined, without constructing 
and operating a power reactor, the answer is “Yes, if it 
is worth developing an alternate to the LMFBR.” Some 
possible reasons for wanting an alternate breeder reactor 
are: (1) The molten-salt system is devoted to the use of 
the thorium reserves and may be better for this 
purpose than the fast reactor; (2) although the LMFBR 
will surely be a technical success, with all the resources 
being devoted to it throughout the world, there may 
still be some reservations about its universal use, such 
as siting restrictions, cost of power, or just plain 
prejudice, and all of these are in evidence to some 
extent; (3) with its complete on-site operation of the 
fuel cycle, practical even for a single unit, the 
molten-salt reactor offers a different mode of industry 
operation, avoiding the necessity for fuel-element 
fabrication and keeping complete control of the fuel 
cycle in the hands of the operator; (4) competition 
between reactor suppliers would have a stronger basis if 
the products they sold were significantly different. 
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There has been gradually increasing industrial 
interest in the United States and Europe’®’ and 
increased interest by the Indian government. Although 
the exact mechanism is not apparent to me, this 
interest could, and I believe should, be focused on a 
reactor construction project, either an MSBE or a 
prototype power reactor. 
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IAEA SYMPOSIUM ON NUCLEAR POWER-PLANT 
CONTROL AND INSTRUMENTATION 


Prague, Czechoslovak Socialist Republic, Jan. 22—26, 1973 


This symposium is designed to provide specialists from Member States an opportunity for a 
comprehensive exchange of information on the current international status of the field. It will emphasize 
problem areas and recent experience in the design, operation, and maintenance of control and 
instrumentation systems, including relevant aspects of safety and economy. 

Topics to be included are: 


@ Nuclear and process instrumentation: design, reliability, and economic considerations 

e Application of experience gained from operating nuclear power plants 

@ On-line computer applications and related operational, reliability, and safety considerations 

© Engineering standards for control and safety systems: considerations of quality control 

@ Installation testing and in-service inspections 

® Potential failure detection: measurement techniques, signal analysis and interpretation, and 
equipment design 


The nomination of a participant will be accepted only if it is presented by the government of a 
Member State of IAEA or an international organization invited to participate. 
For details, persons in the United States should write 


Mr. John H. Kane 

Office of Information Services 
U. S. Atomic Energy Commission 
Washington, D. C. 20545 


and persons outside the United States should write 


Miss C. H. de Mol van Otterloo 
Administrative Secretary 

Division of External Relations 
International Atomic Energy Agency 
Karntner Ring 11 

P. O. Box 590 

A-1011 Vienna, Austria 





CORRECTION 


Wayne K. Lehto, one of the authors of the article on pages 345-353 of Vol. 14, No. 4, of this journal, has 
submitted a correction for the last sentence of the first paragraph in the right-hand column on page 345. 
That sentence says “‘.. . this technique . . . has not as yet been applied to thermal reactors at power... .” 
Lehto has learned that such a reactivity-measurement technique was in fact used at Peach Bottom, 
operating at 25 MW(t). The reference for this added information is: H. Long, R. Schlicht, K. Van Howe, 
and A. Weiman, Use of On-Line Reactivity Computer To Calibrate Control Rods at Power, USAEC 
Report GA-8104, Gulf General Atomic, July 12, 1967. 
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Effects of Metal Swelling and Creep on 
Fast Reactor Design and Performance 


By Paul R. Huebotter* 


Abstract: Dimensional instability of stainless steel, due to 
radiation-induced swelling and radiation-enhanced creep, pro- 
foundly affects the design and performance of fast reactors, 
This article reviews the design literature in an attempt to 
clavify the impact of swelling and creep and to scope the 
options available for accommodating these recently discovered 
fast-neutron irradiation effects, In this critical review the 
problems of swelling and creep are stated for selected reactor 
components, and contemporary-design thinking on each com- 
ponent is presented, The subject is then discussed on a 
reactor-by-reactor basis with either a summary description of 
the design strategy to accommodate swelling and creep or, in 
the case of reactors designed prior to discovery of these 
phenomena, the impact of swelling and creep that has been 
experienced or is foreseen, Finally the probable effect of 
swelling and creep on fast reactor economics is conjectured, 
This effect is translated into an incentive for research and 
development focused toward a metallurgical solution to the 
problem, 


Austenitic (AISI 300 series) stainless steels have long 
enjoyed favored status as structural materials in the 
high-flux high-tempe*-ture regions of fast reactor 
power plants or test facilities.’ Many fast reactors 
(including Experimental Breeder Reactors! and II, 
Enrico Fermi 1, Dounreay Fast Reactor, Rapsodie, and 
BR-5) were designed before irradiation effects on 
stainless steel had been experienced to even 1% of the 
fast-neutron fluence now thought to represent the 
threshold of potentially economic fuel-cycle costs. In 
the liquid-metal-cooled fast breeder reactor (LMFBR), 
this commercial threshold is generally regarded as 
~2 x 107% neutrons/em? (£>0.1 MeV), although 
twice that value is desired for target plants in the 
21000-MW(e) size range. The early reactors, of course, 


*Materials Science Division, Argonne National Laboratory, 
Argonne, Ill. 60439. 
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were built as irradiation test facilities and/or to 
demonstrate breeding, operability of sodium systems, 
technical feasibility of certain fuel-cycle concepts, etc. 
In the design of these early fast reactors, stainless steel 
was assumed to be dimensionally stable other than 
from thermal expansion. The reduction in ductility due 
to irradiation was acknowledged; however, _ this 
ductility loss was generally recognized by designers to 
be partially compensated by strength enhancement (in 
solution-treated material) from neutron bombardment. 
The major objectives in the strategy of core structural 
design were: (1) to maintain stresses sufficiently low to 
not challenge the uncertain limits of ductility and 
(2) to maintain a spatial fuel configuration that is 
constant or changes predictably and safely during 
operating transients. 

The first of the early reactors to pass the fluence 
milestone of 5x 107?  neutrons/em? was the 
Dounreay Fast Reactor (DFR), and, as it did so, a new 
and unexpected irradiation effect became apparent in 
the stainless-steel cladding on experimental mixed- 
oxide fuel elements. Diametral increases that greatly 
exceeded the presumed ductility limit were seen. On 
detailed examination, including immersion-density 
measurements and transmission electron microscopy 
(TEM), swelling of 7% AV/V at 7.8 x 107? neutrons/ 
cm? (all energies) and 500°C was revealed and re- 
ported.” Within months the phenomenon was con- 
firmed in stainless steel irradiated in Experimental 
Breeder Reactor II (EBR-II).?"* 

Since these early papers, additional in-reactor data, 
metallurgical theories, extrapolations, and  out-of- 
reactor simulation in ion-bombardment facilities have 
become abundant. These activities, and the trend of 
progress toward a metallurgical solution to the swelling 


EFFECTS OF SWELLING AND CREEP IN FAST REACTORS 


problem, were recently reviewed by Kangilaski.* 
Although today’s designer has a more complete under- 
standing of fast-neutron irradiation effects than that 
enjoyed by his predecessors, considerable extrapolation 
of existing data is required t_ predict swelling in future 
reactors. Illustrating this point, Fig.1 shows® the 
status ot swelling data vs. the objectives of the U.S. 
LMFBR program (in terms of fast fluence) and the 
EBR-II operating time needed to acquire fluence. 
Actually the required extrapolations are greater than 
suggested in Fig. 1 because all of the high-fluence data 
from EBR-II safety- and control-rod thimbles apply to 
solution-treated 304 stainless steel at temperatures 
below 850°F. Fast Flux Test Facility (FFTF) and 
demonstration plant designers require data on 316 and 
other stainless steels, annealed as well as cold worked, 
to temperatures around 1400°F to cover cladding 
hot-spot conditions with uncertainty factors. 

Thus, since swelling was discovered in 1967, the 
major challenge in the field of fast reactor core design 
has been to predict the nature of dimensional changes 
in reactor components and to accommodate these 
changes by design provisions compatible with the other 
properties of reactor materials, the general charac- 
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program objectives. 
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teristics of fast reactors (including neutronic sensitivity 
to fuel movement), and the objectives and constraints 
of specific reactors. The literature on design implica- 
tions of, and design solutions to, the swelling 
phenomenon was sparse during the first 2 years. This 
silence was due in part to the traumatic and unsettling 
effect of this discovery on “preswelling” approaches to 
core design, including the approaches revealed in the 
two rounds of AEC-sponsored 1000-MW(e) LMFBR 
design studies.71°® Another reason for the dearth of 
literature in the 1968—1969 era was the tendency for 
design judgments to be withheld as having proprietary 
value, in contrast with metallurgical data, theories, and 
extrapolations which largely flowed from government- 
sponsored base-technology programs. In any case, the 
past 2 years have produced a more open disclosure of 
the position of core designers on the swelling issue. 
The phenomenon of radiation-enhanced creep, 
poorly characterized and generally ignored by designers 
before ihe discovery of swelling, has emerged as a very 
important property in certain strategies of designing 
for swelling. It is now the usual practice of core 
designers to discuss both subjects together, as will be 
done in this article. The metallurgical literature on 
radiation-enhanced creep has been recently reviewed 
by Gilbert."7 Core designers of the FFTF have mostly 
used a creep equation furnished by Gilbert and 
Blackburn.’® As in the case of swelling, considerable 
extrapolation of existing data (at least a factor of 10) 
was required to produce a design equation for target- 
fluence conditions in the FFTF and the U. S. demon- 
stration plants. Unfortunately the models for this 
extrapolation are not completely resolved, a fact that 
results in wide differences in high-fluence creep esti- 
mates among researchers in this field. This point is 
illustrated by Fig.2, which shows the 1970 FFTF 
creep expression in terms of creep rate as a function of 
temperature for a single value of stress and two values 
of fluence. The fluence dependence arises from a 
transient term in the design equation which saturates at 
relatively low fluence and contributes little to the total 
creep strain at high fluence and hence nothing to the 
creep rate at that time. Figure 2 also shows some 
Argonne National Laboratory (ANL) estimates of 
radiation-enhanced creep rate under equivalent con- 
ditions.* These estimates, based on a “climb- 


*The comparison is made at an arbitrary stress of 30 ksi. 
At 10 ksi the ANL creep rate would be reduced by a factor of 
9 and the FFTF estimates by a factor of 3. At high 
temperatures (>900°F), this brings the two models into closer 
agreement, but, in the important 600 to 900°F temperature 
range, the models are even more divergent at low stress. 
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controlled” model of radiation-enhanced creep pro- 
posed by Harkness, Tesk, and Li,’°® have changed in 
minor detail as a result of recent work, but the 
comparison illustrates the present dilemma of the 
designer with respect to this important property. 
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Fig. 2 Estimated radiation-enhanced creep rates from two 
different models. 


The present design climate, therefore, is charac- 
terized by uncertainty in the amount of swelling and 
creep that will occur in the next generation of fast 
reactors. Designers with a tendency to use design 
equations literally are reminded by Hafele et al.?° that 
“the existence of a mathematical expression 
should ... mislead nobody about the large errors and 
uncertainties involved.’ Most core designers, including 
those attached to the FFTF project, employ upper and 
lower confidence limits on the swelling and creep 
design equations. Ideally, although this climate of 
uncertainty persists, key reactor components should be 
designed for flexibility to operate satisfactorily over a 
wide range of combined swelling—creep behavior. 
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EFFECTS ON KEY REACTOR 
COMPONENTS 


Although swelling and creep of core structural 
material can affect the design strategy of the entire 
plant, it is convenient to analyze separately the effects 
of these phenomena on individual reactor components. 
In doing this analysis, it is helpful to divide the reactor 
components into two categories as follows: (1) items 
whose service lives are short compared to that of the 
plant, so that they are packaged for routine replace- 
ment by the fuel-handling machine; and (2) items not 
easily changed during the life of a particular plant, so 
that they are designed for 20 to 40 years of service. 

Generally items in the first category have no 
inherent design permanence, and inadequate provisions 
for swelling and creep can be rectified at a cost of a 
few years’ operation with somewhat reduced service 
from these expendable units. Although this cost is not 
trivial, it is considerably less than that of foreshortened 
plant lifetime; thus, where uncertainties in materials 
properties exist, the designer can employ a “boot- 
strapping” strategy rather than one of ultraconservative 
design. Components in this category include: fuel, 
blanket, and control elements; fuel, blanket, and 
control assembly internals; and ducts,* although some- 
what less flexibility may exist here since duct design is 
closely tied to core-restraint design, thus implicating 
hardware in category 2. 

Components in the second category require special 
caution and conservatism in design because, if the 
effects of swelling and creep are badly estimated, the 
operation of the reactor may be interrupted for a long 
maintenance outage or be aborted completely. The 
core-restraint system, load-support structure, and 
instrument tree are examples of components in this 
category. The effects of swelling and creep on the 
design and performance of selected reactor com- 
ponents are reviewed below. 


Fuel Elements 


Prior to the swelling discovery, American fuel- 
element designers were divided in their relative 
preference for strength or ductility in the inevitable 
trade-off between these properties when selecting a 
cladding alloy. The British, on the other hand, have 
always preferred strength, which accounts for the 
appearance of cold-worked stainless steel in the devel- 





*Also known as wrapper tubes, channels, housings, 
shrouds, hex cans, and boxes. Throughout this article, the 
nomenclature used in the FFTF project has been adopted. 
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opment program of the U.K. Atomic Energy 
Authority (UKAEA) prior to the discovery of swelling. 
Now, Bishop generalizes that “cladding of low strength 
and greatest ductility appears to show greatest swelling, 
and vice versa,””' and so the UK designers’ preference 
for strength is even greater. Alloys FV 548, a stabilized 
stainless steel of composition similar to AISI 318, and 
Nimonic PE16, a high-nickel precipitation-type alloy, 
are being developed as future cladding materials. The 
low-swelling tendencies of the latter, if ion- 
bombardment results are confirmed by high-fluence 
neutron irradiation, could make it especially valuable 
as a future duct material. 

The empirically established swelling resistance of 
cold-worked 316 stainless steel has produced a general 
consensus favoring its use for fuel-element cladding in 
near-term facilities in the United States, although other 
stainless steels, especially carbon-stabilized ones, are 
being evaluated.” ? 

The Edison Electric Institute (EEI), speaking for 
potential U.S. customers for fast reactor power plants, 
sounded an early alarm about the implications of the 
swelling discovery.2* The EEI panel focused their 
concern on the issue of fuel-element lifetime, sug- 
gesting the need to discount successful irradiations of 
mixed-oxide fuel elements in EBR-II and DFR by a 
factor of 3 to 4 in burnup, equivalent to the 
neutron-flux differential between these small test 
reactors and the target commercial plants. Figure 3 





rryttl 











a 
= Assumptions 

[— Peak linear power = 16 kW/ft 
|— Fuel-element OD = 0.25 in. 

|— Cladding thickness = 0.015 in. 
Smeared density = 9.35 g/cm3 
— (85% T.D.) 

Cladding temp. at place of max. 
|__ Swelling = 1000°F 
AD/D =4 AVIV 


ae Seay 


! 


_ 
o 


ie i 





Li til 








| 
| 





MAXIMUM AD/D, % 

Ss 

So 
= 
= 
a 
o 
=) 
" 
bd 
o 
Pad 
r=) 
a 


Litiiul 


| 
| 











TT) ee oe 
103 104 105 108 
PEAK FUEL BURNUP, MWd/metric ton 





Fig. 3. Cladding dilation from radiation-induced swelling in 
three fast reactors. 
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illustrates this point showing swelling-induced cladding 
dilation as a function of burnup for a typical mixed- 
oxide fuel element in the EBR-II, the FFTF, and a 
commercial LMFBR. The estimates of both absolute 
and relative swelling depend on the swelling correlation 
in use,* but clearly more cladding swelling per unit 
burnup is to be expected in reactors with greater 
high-energy flux. 

For the past 3 years, two important issues relevant 
to the EEI concern have been debated. These are: 


1. Is the additional cladding @‘lation good or bad in 
terms of fuel-element lifetime? 

2. How much of the diametral increase measured in 
test elements is due to swelling and radiation-enhanced 
creep (thought by many to be innocuous), and how 
much is due to fuel—cladding mechanical interaction 
and thermal creep from gas pressure (the accepted 
failure causes)? 


Selected excerpts from that debate are given 
below: 


Bump?‘ : “One consequence of cladding swelling is 
the provision of more room for gas in the fuel, thus 
lowering the pressure of that gas and, in turn, the 
pressure exerted on the cladding by fuel swelling.” 

Boltax et al.?°: “Figure 4a shows the [cladding] 
stress accumulation when stainless-steel swelling is 
comparatively small ...e.g., in EBR-Il and DFR.... 
Figure 4b shows the behavior of the hoop stress when 
stainless-steel swelling is large, as in the case of high 
fluence/burnup pins [in a commercial LMFBR].” 

Bagley et al.?°: “It has been found ... that pin 
strain normally correlates well with clad swelling so 
that void-development tendencies may validly be in- 
ferred from fuel pin [diametral increase] .” 

Jackson et al.27 The measured diameter 
increase of these two fuel pins [PNL 2-3 and 2-6] is 
almost entirely a result of cladding swelling, and the 
small difference between measured and predicted 
diameter is error.” 

Neimark et al.?® : “Mechanisms other than cladding 
swelling alone are contributing to the diametral in- 
creases [of NUMEC Group C elements with burnup of 
~10%] .” 

Craig et al.2°: “Smeared fuel density appears to 
have a major effect on the diametral increases of these 
[NUMEC D-5 and General Electric F2 series with 
burnups up to ~12%] fuel pins.” 


*The estimates in Fig. 3 conservatively assume a fluence 
exponent of 1.5 and no saturation of swelling beyond the 
range of current experience. 
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Boltax et al.>°: “The results [of computer studies] 
reveal the major effect of fuel pellet density on total 
plastic strain of the cladding. The benefits obtained 
from lowering the fuel pellet density are reduced 
significantly at cladding temperatures where thermal 
creep becomes important (T > $50 C),” 
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Fig. 4 Postulated cladding-stress history in EBR-II] and a 
commercial LMFBR. 


In a recent review paper on mixed-oxide fuel 
technology for the LMFBR, Leggett?! concludes that: 
“Cladding swelling presents few significant problems to 
fuel pin performance [emphasis in original] but must 
be completely defined for assembly design and com- 
patibility with the rest of the core components.” In 
another recent assessment of the critical factors in 
mixed-oxide fuel-element design, Roake?? observes 
that: “We have come full circle since 1966; i.e., pin 
design life was governed (in chronological sequence) by 
(1) fission gas pressure and fuel swelling, (2) stainless- 
steel swelling, (3) radiation-enhanced creep and swell- 
ing gradients, and now, but at a higher burnup, (4) fuel 
swelling.” 
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It is interesting that, other than the recognized 
need to consider swelling resistance as a major criterion 
for cladding selection, the debate over the effects of 
swelling on fuel-element performance has produced 
negligible fuel-element design changes. Table 1 gives 
some of the principal design parameters for fuel 
elements in oxide-fueled sodium-cooled fast reactors of 
comparable size. Note that there is no obvious pattern 
of design difference between fuel elements designed 
before the swelling—creep problem came to light 
(BN-350, PFR, and Phenix) and those designed after 
the problem was recognized (SNR, FFTF, and the 
U.S. demonstration plants designed by Atomics Inter- 
national, General Electric, and Westinghouse).* 

The current consensus of both development and 
design people is that swelling should not be regarded 
primarily as a “cladding problem” or a “fuel-element 
problem.” It is also recognized, however, that a burnup 
limit on fuel elements could result if designers do not 
properly cope with swelling in the design of other core 
components. Fuel-element designers now work on the 
optimistic assumption that core-swelling effects will 
not limit fuel burnup, and fuel elements are designed as 
the components pivotal to core lifetime; hence, essen- 
tially the same as before swelling was discovered. If this 
assumption proves wrong, and allowable fuel burnup is 
dictated by core-swelling effects, the economic loss can 
be partially recouped through the use of higher 
smeared density and the relaxation of fabrication 
specifications that purchase burnup capability beyond 
the usable range.° 


Fuel Assemblies 


Table 2 presents selected fuel-assembly design 
parameters for the same group of oxide-fueled sodium- 
cooled fast reactors. As in the case of fuel elements, no 
revolution in design practice appears to have resulted 
from the discoveries of swelling and creep. 

The wire-wrap spacer system has been chosen in 
both the preswelling and postswelling design eras, as 
has the principal alternative—compliant grids. Neither 
is free of problems arising from swelling and creep. 
Both systems must be carefully designed for accom- 
modation of differential swelling between the fuel- 





*One of the few apparent correlations in Table 1 is that 
most reactors in the first group have bottom-end gas plenums 
and most in the second group have top-end plenums. Since the 
bowing of fuel assemblies is exaggerated by top-end plenums, 
this difference is clearly not a conscious reaction to the 
swelling problem. Instead, it appears to be a U. S. vs. European 
difference of opinion on the safety implications of gas stored 
upstream from the core. 
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Table 1 Fuel-Element Design Parameters in Oxide-Fueled LMFBRs Under Design and Construction* 





Reactors designed after swelling was discovered 














Reactors designed before —— 
swelling was discovered : nae 
Al demon- GEdemon- demon- 
Parameter BN-350 PFR Phenix SNR FFTF stration stration stration 
Fuel form Annealed Annealed Solid Solid Solid Solid 
pellets pellets pellets pellets pellets pellets 
Smeared density, % 73 80 78 80 85.5 80 
Target burnup (av.), 55 45 tb 407/844 ie 
MWd/kg 
Target burnup (peak), 46 69 50+/100¢ 80 80 110 124% 103 
MWd/kg 
Linear power (peak), 13.4 BS;1 14.1 13.8 13 14.5 11.7+/14.5¢ 13.8 
kW/ft 
Location of fission-gas Top Bottom Bottom Bottom Top Top Topt/vent 
plenum at tops 
Fission-gas plenum <3 28 (red. 16.5 25.6 42 26 30 
length, in. diameter) 
Fuel length, in. 41.8 36 Bo A 37.4 36 44 30 35 
Cladding material “Russian 316S.S., 316 S.S. Open 316 S.S., 316S.8., 3218S. 316SS., 
316” 20% C.W. 20% C.W. C.W. 20% C.W. 
Fuel-element diameter, 0.240 0.23 0.260 0.236 0.23 0.27 0.25 0.25 
in. 
Cladding thickness, 0.014 0.015 0.018 0.015 0.015 0.016 0.015 0.016 
in. 
*Missing values are not available in the open literature. + First core. Target core. 


Table 2 Fuel-Assembly Design Parameters in Oxide-Fueled LMFBRs Under Design and Construction* 





Reactors designed after swelling was discovered 











Reactors designed before waateg- 
swelling was discovered : — 
Al demon- GEdemon- demon- 
Parameter BN-350 PFR Phenix SNR FFTF stration stration stration 
Type of hold-down Hydraulic Hydraulic Hydraulic Hydraulic Hydraulic Hydraulic Hydraulic Mechanical 
Sodium velocity (max.), 26 20 (free 24 20 (free 25 312 20.6 
ft/sec stream) stream) (free 
stream) 

Type of duct material “Russian S.S., C.W. 316S.S. Open 316.85. 316SS., 3218S. 316SS., 
316” 20% C.W. C.W. 20% C.W. 

Fuel-element pitch/ PIS 1.26 1.16 1.29 1.24 1.17 1.28 1.24 

diameter (dimensionless) 

Type of fuel-element support Wire wrap Honeycomb Wire Grid Wire Wire Wire Honeycomb 
(alt. grid wrap wrap wrap wrap grid 
pins) (ref.) 

No. of fuel elements/ 169 325 ZiT 169 217 271 127 217 

assembly 

Fuel-assembly width, 3.780 5.54 4.87 4.34 4.575 5.570 3.849 

A/F, in. 

Fuel-assembly in-core 3.859 ar 5.01 4.715 5.887 4.009 §.922 

pitch, in. 

Duct thickness, in, 0.079 0.115 0.138 0.110 0.120 0.125 0.086 

(min.) 





*Missing values are not available in the open literature. 
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element bundle and the duct, and the wire-wrap system 
must also accommodate some differential swelling 
between the wire and cladding. Greater swelling of the 
cladding than an initially taut wire has produced 
pronounced helical distortion of high-burnup fuel 
elements?® but no known failures in the U.S. fuel- 
development program. These swelling differentials can 
arise from temperature differentials (typical values are 
shown in Fig. 5, taken from Ref.33) or from dif- 
ferences in composition or temper possibly so subtle as 
to be presently unrecognized as significant variables in 
swelling behavior. A problem in grid design is how to 
attach it to the duct. Creep will tend to round the 
ducts by as much as 0.070in. per flat, based on 
contemporary creep estimates* and fuel-assembly 
design parameters, and this rounding could distort a 
grid attached directly to the duct. Hanger straps and 
other concepts that permit duct—grid differential 
deflection are under development to solve this recently 
recognized problem. In Core A of Enrico Fermi 1, grids 
are fixed to vertical tie rods, thus creating a “birdcage” 
structure which slips into the duct and is anchored to it 
at the lower end.** This avoids the problem of direct 
attachment to the duct but tends to promote more 
peripheral flow than is desired. Grids were chosen for 
Fermi 1 because of the use of a square pitch in which 
the wire wrap is unstable. The EBR-II, designed in the 
same era, uses wire wraps and a triangular pitch.?° 


Speaking for the British designers, whose Prototype 
Fast Reactor (PFR) uses grids attached to the duct for 
fuel-element support, Bishop?! offers the opinion that 
it is less difficult to avoid the undesirable extremes of 
mechanical interference and vibration in the gridded 
concept than in the concepts that employ non- 
compliant fins or wires on the cladding. Westinghouse, 
with considerable experience with gridded spacer 
systems in pressurized-water reactors and in the U—10 
wt.% Mo core of Enrico Fermi 1, shares the British 
preference for grids in LMFBRs. They contend that, if 
the duct, cladding, and wire should follow swelling 
laws that differ from one another in the adverse 
direction, metal swelling could impose a burnup limit 
on the wire-wrap concept as low as 30,000 MWd/ 
metric ton in a high-flux LMFBR.+?7 Wire-wrap 





*The EBR-II driver-fuel duct rounding has been observed 
at a total fluence of only 1.9 x 107? neutrons/cm?, and it 
agrees with calculated rounding within the range and accuracy 
of measured values.?4 

*This analysis did not include the effects of radiation- 
enhanced creep which tend to mitigate the consequences of 
adverse differential swelling. 
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Fig. 5 Typical axial profile of average metal temperatures in 
an LMFBR fuel assembly. 


proponents contend that this conclusion is valid only if 
the wire-wrap design is unnecessarily handicapped by 
little or no internal clearance at beginning of life. 
General Electric (GE) core designers have selected the 
wire wrap as their reference demonstration-plant fuel- 
element support concept, although a grid spacer 
concept which is “potentially more attractive than the 
wire wrap with respect to differential swelling, fuel 
performance and thermal hydraulics” is under develop- 
ment.?® The reference GE wire-wrap system provides 
radial compliance to accommodate differential swelling 
by: (1) initial clearances built into the assembly, 
(2) staggered-start wires, and (3) omitted wires on 
selected elements. Atomics International (AI) has 
unequivocally chosen the wire-wrap system as a result 
of their fabrication and test program, which satisfied 
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them that the small-amplitude bowing of the mul- 
tiplicity of fuel elements within the duct will exert 
sufficient mechanical interaction to avoid vibration, 
fretting, and wear even when the planar clearance is 
enough to absorb differential swelling.*° Atomics 
International prefers the wire wrap because of lower 
cost, twice the interchannel mixing, and only two- 
thirds as much pressure drop compared with gridded 
designs.*° 

As part of a general assessment of the design, 
research, and development implications of metal 
swelling in fast reactors, Huebotter et al.® concluded 
that: (1) sufficient design accommodation could safely 
be made for the differential axial and radial swelling 
likely to occur in wire-wrap support systems in 
near-term LMFBRs; (2) the fuel-element support sys- 
tem is unlikely to emerge as the burnup-limiting 
component of near-term cores; (3) the choice of a 
fuel-element support system is not an irrevocable 
30-year decision but is one that can be reconsidered for 
follow-on cores when metal swelling and other 
pertinent behaviorisms are better characterized; and 
(4) since the adequacy of the wire-wrap system at 
fluence targets for commercial fast reactors is still a 
matter of industrial disagreement, opportunities should 
be sought to test representative grid designs. On the 
basis of such reasoning, the FFTF project has chosen 
the wire-wrap concept for initial use, but the U.S. 
LMFBR program is supporting grid-development 
projects at GE and Westinghouse, which will both 
include testing in EBR-II as adjuncts to fuel-element 
irradiation programs. 

In principle, differential dilation between the duct 
and the contained fuel elements can produce cladding 
overheating by restriction of coolant passages. In 
practice, however, this effect is of minor concern. It 
has been estimated, for example, that the hot-spot 
temperature of the PFR cladding is increased only 
11°F due to a uniform 1% increase in fuel-element 
diameter over the entire 36-in. core length with no 
compensating duct dilation.?’ In actual reactor opera- 
tion the effect would be even less for either whole-core 
or partial-core refueling strategies. With whole-core 
refueling the differential dilation develops in all as- 
semblies at approximately the same rate. Therefore a 
general increase in hydraulic resistance would be 
compensated by increased pump speed to maintain 
constant core flow and mixed-mean outlet tempera- 
ture. With partial-core refueling, increased hydraulic 
resistance in the older fuel assemblies would tend to be 
offset by their decline in heat generation due to 
fissile-atom burnup. Leduc etal.*' point out that 
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bowing of the duct relative to the fuel elements 
increases the overall hot-channel factor for the cladding 
by changing random (statistical) geometric aberrations 
into systematic (multiplicative) ones. Other than in 
fuel-element bundles with a very low pitch/diameter 
ratio, however, swelling is not regarded as a significant 
factor in normal heat removal from the LMFBR core. 

A design-sensitivity study by Burchill and 
Menzel*? has shown that, if the successful accommo- 
dation of swelling and creep requires some form of 
positive core restraint (see below), fuel-assembly 
dimensions might well be influenced. Hoop stresses, for 
example, vary as the duct diameter cubed and inversely 
as the pad thickness squared. In large assemblies, duct 
stresses can be the limiting factor in design. 


Core Restraint 


Observers on both sides of the Atlantic? !*4? regard 
the core-restraint system as presenting the greatest 
design challenge with respect to successful accommoda- 
tion of swelling and creep in core stainless steel. In 
fact, the term “core restraint” has found common 
usage only since the discovery of the swelling 
phenomenon, implying the recognition of new 
problems and a higher degree of design emphasis. 
Unlike fuel elements and fuel assemblies, the core- 
restraint system implicates reactor internals not easily 
changed during the life of the plant. Therefore the 
strategy of periodic upgrading of design as more and 
better materials data are acquired does not pertain to 
the core-restraint system as it does to expendable core 
components. 

The functional requirements of the core-restraint 
system can be summarized as follows: 

1. Reactivity Control, Produce a safe, calculable, 
and reproducible structural response of the core to 
steady-state and transient operation. 

2. Alignment, Maintain desired juxtaposition 
between the in-core and above-core components of the 
nuclear control system and between reactor-assembly 
effluent streams and downstream instrument sensors. 

3. Refueling. Provide adequate clearance and 
achieve adequate duct straightness for reactor assem- 
blies to be inserted and removed within an acceptable 
range of vertical forces. 


Fast reactor core-restraint systems, including those 
designed before swelling and radiation-enhanced creep 
were discovered, fall into three categories. For each of 
these categories, the essential design features and 
general response to swelling and creep are discussed 
below. 
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Positive Restraint. Designsin this category employ 
external mechanisms, sucn as the hold-down assembly 
in Enrico Fermil or the peripheral clamps in the 
FFTF, to gather at least the high-worth reactor 
assemblies into a tight array before criticality is 
approached. With this tight cold array, in which the 
centerline location of each assembly is known at 
specific elevations, one of which is in or near the core, 
structural reactivity effects are calculable with an 
accuracy limited only by uncertainties in spatial 
temperature distribution. Alignments are assured, 
usually by making the downstream end of the assembly 
one of the elevations of fixed radial position, and 
refueling clearances are produced by release or with- 
drawal of the restraint mechanism. 

In positive-restraint systems the characteristic inter- 
action between thermal expansion effects, swelling, 
and creep is shown** in Fig.6. When only thermal 
strains are considered, the bowing profile of all 
assemblies except that at the core center line is 
qualitatively represented by Fig. 6(a). Differential 
swelling alone would produce such shapes as that 
shown in Fig. 6(b), and the additive effect of thermal 
and swelling strains is represented by Fig. 6(c). In 
suppressing free bowing in the direction that would be 
produced by differential thermal elongation and dif- 
ferential swelling elongation, the core-restraint system 
imposes high bending stresses. The relaxation of these 
stresses, by radiation-enhanced creep, alters the bowing 
profile markedly [Fig. 6(d)]. The shapes qualitatively 
depicted in Fig. 6 are a function of time, location in 
the reactor, and the design equations for swelling and 
creep. 
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Fig.6 Assembly-deflection patterns characteristic of positive 
core restraint. 
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When the core is released for fuel handling, the | 
assemblies assume shapes of a basically different 
character, as shown** in Fig. 7. Again, it is convenient 
to analyze separately the three factors responsible for 
distortion. With the reactor at zero power, the thermal 
strains are removed, and the assemblies would be 
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Fig. 7 Typical assembly-deflection patterns after removal of 
positive core restraint. 


straight [Fig. 7(@@)] if no other distortion mechanisms 
were operative. Differential swelling is not reversible, 
however, and so the upper ends would tend to splay 
outward [Figs. 7(b) and 7(c)] when freed of restraint. 
With sufficient stress relaxation at power, however, the 
profile at zero power can actually be reversed 
[Fig. 7(d)|. The position of the top end of each 
reactor assembly, which must be known within fairly 
close limits for blind fuel handling,* is a strong 
function of swelling and radiation-enhanced creep. 


Passive Restraint, Designs in this category, which 
include EBR-II, Rapsodie, and Phenix, have the follow- 
ing generic features: 

1. Core assemblies are cantilevered from their inlet 
nozzles and are vertically aligned by two-point support 
thereon. 

2.The top ends of core assemblies are radially 
unrestrained other than by mechanical interference 


*In existing fast reactors the allowable uncertainty in 
handling-head position in the 7,6 plane is 0.25 to 0.5 in. 
Improved fuel-handling designs are being developed to increase 
the allowable uncertainty to 1 in. or more. 
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with, and the longitudinal stiffness of, outboard rows 
of blanket and/or reflector assemblies. Forces tighten- 
ing the core at a spacer-pad plane are generated when 
“free-flowering” edge-core assemblies are opposed at 
the top ends by these straight outboard rows. 

3. Refueling clearance is provided exclusively by 
the cold clearances between spacer pads, which are 
typically at one elevation in or near the core. 


Mechanical interactions are produced between 
assemblies in the passive-restraint concept, and these 
are capable of inducing sufficient stress to affect the 
final shape of assemblies through radiation-enhanced 
creep. Figure 8 illustrates*® the calculated time- 
dependent mechanical interactions affecting a high- 
fluence test assembly in EBR-II. Compared to the 
positive-restraint designs, however, one is less likely to 
encounter reactor assemblies with creep-dominated 
residual shapes because the thermal bowing is unre- 
strained over at least part of the power range. An ideal 
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design strategy has been suggested*® to be one in 
which all or most of the thermal bowing is allowed to 
occur without restraint while resisting the additional 
bowing caused by differential swelling. This permits 
creep to mitigate the permanent bowing of core 
assemblies without risking the concave (relative to the 
core center line) bowing produced by relaxation of 
stresses incurred in suppressing thermal bowing. These 
concave curvatures can produce an undesirably loose 
core at low power. 


No Restraint. In designs of this category, assem- 
blies are cantilevered from their nozzles* or from a 
support point between the nozzles and the corey and 
are unrestrained at the top end over the entire power 
range. The structural reactivity effect is dependably 





*An example is the Gulf General Atomic gas-cooled fast 
reactor demonstration plant.*7? 
+An example is the British Prototype Fast Reactor.*® 
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Fig. 8 Time-dependent mechanical interactions between assemblies in EBR-II. 
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negative since thermal bowing moves fuel away from 
the core center line. 

In this approach, creep is not a design factor 
because bending stresses are not produced. Thus the 
assemblies assume the residual swelling-dominated 
curvatures [Fig. 7(c)] resulting from flux and tempera- 
ture gradients. Periodic rotation is generally envisaged 
as the operational strategy to preclude an unacceptable 
degree of bow from the standpoint of alignment and 
fuel handling. Whether periodic rotation will be an 
acceptable operational practice in commercial fast 
reactors will depend on the frequency with which it 
must be done. 

Both thermal bowing and bowing from differential 
swelling are a strong function of duct length down- 
stream from the core. Figure 9 shows® the typical 
contribution to overall bowing from the portions of 
the duct containing the fuel, the upper axial blanket, 
and a downstream gas plenum. For minimization of 
bowing, particularly important in this category of 
core-restraint design concepts, the upstream gas 
plenum is advantageous, 


Other Reactor Components 


The extent to which other reactor components are 
implicated in the general problem of radiation-induced 
swelling and radiation-enhanced creep of stainless steel 
is difficult to assess in general terms because of wide 
variations in the design of reactor internals other than 
those discussed above. Some specific observations, such 
as those given below, can be made, however. 

1. Elongation of assemblies, and differences in 
elongation, will complicate the design of fuel-handling 
equipment somewhat. Figure 10 shows® end-of-life 
elongations of reactor assemblies in the AI 1000-MW(e) 
design’? based on the 1969 FFTF swelling predictions 
for 20% cold-worked 316 stainless steel.2% More 
flexibility in the fuel-handling equipment is needed 
than when the assembly lengths are constant and 
uniform. With only a 14-in.-long core, EBR-II duct 
elongations have been limited to % + 4 in., but these 
have caused handling difficulties because they were 
unanticipated in the design.4° 

2.Some reactors employ a sweep arm or similar 
device to assure that assemblies are properly seated on 
the grid plate. Length increases will tend to mask the 
evidence of improper seating, particularly in a core 
containing both old (long) and new (short) assemblies. 

3. Aside from its incompatibility with significant 
length increases, the EBR-II fuel-handling system has 
been given much of the credit for successful operation 
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Fig. 9 Typical contributions of fuel and above-fuel length to 
free bowing. 


despite unanticipated swelling effects.5° The key 
feature is the hold-down mechanism that spreads the 
six assemblies adjacent to the one being inserted or 
withdrawn, thus creating the necessary local clearance. 

4. Swelling of the grid plate, instrument plate, and 
other core-proximity “permanent” hardware will 
probably be avoided by the fluence limits imposed to 
retain the desired ductility (generally < 107? neutrons/ 
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Fig. 10 Component elongations due to swelling in a 1000-MW(e) LMFBR. 


cm?), Extra shielding of components beyond the ends 
of the assemblies is conveniently purchased with small 
increases in axial blanket lengths. Figure 11 shows® the 
effect of blanket length on grid-plate and instrument- 
plate flux in the AI 1000-MW(e) design.'? Grid plates, 
whose designs are generally governed by deflection 
criteria rather than stress criteria, cannot be permitted 
to undergo significant differential swelling. Fortunately 
theoreticians>! believe that low flux will produce less 
swelling than high flux to the same fluence. This works 
to the advantage of the core-proximity hardware of 
long life expectancy. 


EFFECTS ON SOME PRESENT AND 
FUTURE FAST REACTORS 


The effects of swelling and creep have been felt in 
operating fast reactors, which in no case were designed 
to accommodate these effects but can still yield 
valuable engineering insight into how to design future 
reactors. This group is represented by EBR-II and 
Enrico Fermi 1. 


In the next grouping are the Russian BN-350, the 
British PFR, and the French Phenix, all scheduled for 
criticality between late 1971 and early 1973. These 
facilities were in an advanced stage of design when the 
swelling—creep properties gained recognition. In 
general, these phenomena were not permitted to 
perturb the designs significantly. 

Finally, indications of how fast reactor designers 
have responded to the problems of swelling and creep 
in an overall sense are available in the literature 
describing reactors still in design or in the early stages 
of construction. This group of reactors includes the 
German-Benelux SNR, the FFTF, and the U.S. 
demonstration plants being designed by Atomics inter- 
national, Combustion Engineering, General Electric, 
Gulf General Atomic, and Westinghouse. 

As we have already seen, major differences in 
design philosophy are not found in the design of fuel 
elements and assemblies. Therefore this reactor-by- 
reactor assessment of the impact of swelling and creep 
on design and performance will concentrate on the 
core-restraint system and the other reactor internals 
sensitive to the swelling—creep phenomena. 
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Fig. 11 Effect of axial blanket thickness on flux at above- and 
below-core components. 


EBR-II (Refs. 36, 50) 


The EBR-II core plan employs hexagonal ducts of 
2.290 in. across outside flats, with a nominal cold 
clearance between assemblies of 0.030 in. One %-in.- 
diameter spacer pad (produced by indenting the inside 
surface of the duct) is provided at the center of each 
flat about 1 in. above the core midplane. The nominal 
height of each pad is 0.014 in., and so the nominal gap 
between pads is 0.002 in. cold and zero at operating 
temperature. 

‘The lower adapter (nozzle) of each assembly fits 
into locating holes in two horizontal support (grid) 
plates. The upper plate supports the weight of the 
assembly and provides a pivot point by means of a 
spherical seat. The amount of free pivoting is limited 
mainly by the 0.004-in. radial clearance in the lower 
grid-plate hole. Core radial restraint of the passive type 
is provided by the eight to nine rows of radial reflector 
and/or blanket assemblies similarly cantilevered from 
their respective inlet nozzles and also equipped with 
spacer pads. A retaining ring 0.030 in. beyond the final 
row, near the top of the assemblies, is also present. 
This ring was installed as an aid to peripheral-row fuel 
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handling and is believed to play no essential role in 
core restraint. The use of hydraulic hold-down, 
pioneered in EBR-II, circumvents the need for top-end 
fixity, so the assemblies are free to deflect under the 
influences of differential swelling, creep, and thermal 
gradients. The latter (the only effect recognized during 
the design phase of the EBR-II project) produces 
top-end interferences just beyond the core edge, 
creating inwardly directed forces at the pad plane 
according to calculations employing the BOW-V 
code.5* These forces maintain a tight core during 
operation and also produce a small positive structural- 
reactivity coefficient owing to the compressibility of 
mid-flat pads on the relatively thin (0.040-in.) ducts. 
This positive component is greatly overridden by fuel 
and sodium expansion and other negative components 
of the power coefficient.5* Most of the above- 
described design features and operating characteristics 
of EBR-II are illustrated® schematically in Fig. 12, a 
preswelling-era conception of EBR-II core-component 
bowing behavior.5* 

The reactor has performed well in the face of 
radiation-induced swelling and radiation-enhanced 
creep, neither phenomenon having been recognized 
until 10 years after design, and the plant factor has 
increased steadily during 7 years of operation. Core 
components not limited by fuel-burnup restrictions, 
ie., control- and safety-rod “thimbles” (guide tubes), 
have functioned to world-record exposures up to 
1.7 x 10? neutrons/cm? (all energies) before swelling 
effects forced their removal. Although the thimbles 
operate at relatively low temperature (700 to 850°F), 
this solution-treated* 304 stainless steel has 
exhibited** up to 11.7% AV/V, believed to be the 
highest yet seen from neutron exposure. 

Curiously, elongation and dilation, rather than 
bowing, eventually forced the removal of these high- 
fluence components. Elongations of % +16 in. were 
experienced,*® well beyond anything anticipated in 
the design of the fuel-handling equipment. Dilation up 
to 0.090 in. has been seen,*? more than enough to fill 
the nominal 0.030-in. gaps separating thimbles from 
neighboring assemblies which, for the most part, 
contain low-burnup driver fuel and therefore ex- 
perience little duct dilation to further aggravate the 
local congestion. 

Examination of high-fluence ducts and thimbles 
from EBR-II has produced two “surprises,” both 





*These hexagonal thimbles, like the ducts on driver and 
fuel-test assemblies, contain some cold work at corners and 
spacer pads but are essentially in an annealed condition. 
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Fig. 12 Thermal deflections of EBR-II assemblies. 


relevant to the search for a simple and reliable 
core-restraint design to accommodate the calculated 
effects of swelling. These are (1) the amazing degree of 
longitudinal straightness of these components and 
(2) the near constancy of duct diameter measured 
across spacer pads (although these pads are located on 
the highest swelling portion of the duct) while con- 
siderable duct dilation is found % in. away. 


The first phenomenon is illustrated*® by Fig. 13. 
In the three high-fluence thimbles removed from the 
reactor in July 1970, note that measured bowing was 
less than predicted bowing by factors of 6 to 12, 
depending on the swelling expression used for the 
predictions. Thimbles discharged more recently show 
the same tendency of straightness beyond any reason- 
able expectation. On the basis of the only thimble 
completely characterized by destructive measurements 
of density as of this writing (5A3), radiation-enhanced 
creep is not responsible for this straightness, although 
the straightening effects of creep have been calculated 
on a high-fluence fuel-test assembly.** The second 
surprise, that associated with spacer-pad behavior, is 
illustrated*® in Fig. 14. This is being evaluated as a 


creep phenomenon, and one that gives rise to hope that 
reactors with in-core spacer pads are not doomed to 
low-burnup capability. 

Both of the above-described phenomena in EBR-II, 
as well as fuel-element longitudinal bow and duct 
rounding, are being studied at ANL to test the 
reliability of contemporary swelling and creep ex- 
pressions, albeit on solution-treated 304 stainless 
steel.24°45»5°63 Also, of course, the surprising 
capability of the passive EBR-II core-restraint system 
to maintain straight components with absolute swelling 
up to 11.7% by volume must be understood before it 
can be reliably extrapolated to the next generation of 
fast reactors, as some experimenters have suggested 
doing.6-46+64 


Enrico Fermi 1 (Refs. 35, 65) 


Enrico Fermi 1, designed in the same era as EBR-II, 
has a more positive means of producing core radial 
restraint, not unlike the clamped-core concepts that 
have emerged since the discovery of swelling, except 
that radial restraint is produced by fixing the top ends 
of core assemblies rather than by direct side loading. 
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Fig. 13 Measured and calculated bowing of EBR-II safety and 
control thimbles. 


In this reactor the square cold-worked 347 
stainless-steel ducts of 2.646-in. width and 0.096-in. 
thickness are vertically oriented in a double support 
plate that receives tubular nozzles in close-fitting holes. 
A plane of spacer pads exists 4.5in. above the 
midplane of the 30.5-in.-long core. These pads, made 
of Colmonoy, are welded to the corners of the ducts 
and machined to close tolerances to produce the 
2.693-in. square pitch of core assemblies. 


The handling head of each core and inner radial 
blanket assembly is captured by a stiff hold-down finger 
projecting downward from a horizontal hold-down 
plate. This finger prevents lifting from hydraulic forces 
and also precludes radial displacement of the 
assembly’s upper end from thermal bowing or dif- 
ferential swelling. The center-to-center spacing of 
hold-down fingers is slightly less than the center-to- 
center spacing of support-plate holes. Accordingly, 
when the hold-down plate is lowered (as it always is 
before reactor startup) and the handling heads are 
captured by their respective hold-down fingers, an 
inward gathering of assemblies occurs, and all clearance 
at the spacer-pad plane is eliminated. 


This design, known as the “tight-core concept” is 
illustrated®® schematically in Fig. 15. The design has 
produced very dependable reactivity effects, including 
a negative structural component of the power coeffi- 
cient, which was the primary design objective. Un- 
fortunately the operating experience of this reactor has 
been limited to low fluence, and so the capability of 
this design to accommodate swelling effects has not 
been tested. The most pertinent analogous experience 
resulted from the October 1966 fuel-melting incident, 
during which brief period a few core assemblies became 
badly distorted as a result of high thermal gradients.°® 


The potential of this reactor to accommodate 
swelling and creep effects has been estimated.® No 
significant swelling would be expected during the 
remainder of the reactor operations with low-burnup 
metallic Core A fuel, although creep effects may be 
seen, probably as residual longitudinal bowing concave 
with respect to the reactor center line. When the 
proposed UO, core is installed and operated at its 
intended 400-MW(t) power level, the peak total flux 
(6.5 x 10'5 neutrons/cm?) will be such as to produce 
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Fig. 15 Core-restraint design in Enrico Fermi 1. 


a burnup/fluence ratio similar to that of future 
LMFBRs, and swelling could become a major factor. 
Duct elongation due to swelling can be accommo- 
dated by compression of the nozzle spring, a feature 
originally intended to provide for differential thermal 
expansion. Clearance for dilation accommodation, 
assuming the use of the full 0.047-in. interassembly 
gap, is slightly greater than EBR-II.* Because of the 
in-core spacer pads, however, the usefulness of the full 
interassembly gap is not assured. Substantial creep 
distortion around the pads is required. Although this 
phenomenon is not understood sufficiently to assure 
that dilation accommodation will not be restricted by 
in-core pads, the indications are favorable. The stress 
required for the creep rate to match the swelling rate in 
Fermi! has been estimated® to be less than that 
required in EBR-II.+ At this time it appears that the 





*The interassembly gaps in Fermi 1 and EBR-II represent 
linear clearances of 1.77 and 1.31%, respectively. 

+This result is highly dependent on the creep and swelling 
expressions used. The calculations referred to were based on 
expressions that yield lower swelling and higher creep estimates 
for the cold-worked Fermi ducts relative to the solution- 
treated EBR-II ducts under the same conditions. 


limiting factor in the ability of the Fermi design to 
accommodate swelling and creep is the necessity of the 
handling heads to mate with the hold-down fingers 
when the hold-down plate is lowered following refuel- 
ing and prior to startup. This operational requirement 
constrains the allowable bowing due to swelling and 
creep. The nominal allowable mismatch is ¥ in. An 
inspection facility has been built, which, among other 
things, can provide direct measurements of assembly 
straightness. The facility is a portable, heated, shielded 
hot cell that can be stationed at the mouth of the 
reactor exit port to receive assemblies one at a time. 
Originally designed and built for the requalification of 
irradiated fuel assemblies°? suspected of having 
sustained minor damage in the October 1966 fuel- 
melting incident, this facility can be used to monitor 
the effects of swelling and creep, although a more 
efficient procedure has been developed to monitor 
bowing alone. This procedure involves the use of an 
under-sodium periscope built into a diving-bell struc- 
ture that is lowered down the exit port onto the 
handling head of an assembly that has been placed in a 
special gauging pot in the transfer rotor. Bowing is 
measured optically, ie., by photoanalysis, to a 
quoted®® precision of +0.015 in. in magnitude and +5° 
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in angular direction. This procedure, already used to 
measure bowing in a complete core loading, offers the 
capability for further study of integral swelling—creep 
effects as reactor operation proceeds. 


BN-350 (Refs. 69, 70) 


The BN-350, being built at the Caspian Sea in the 
USSR, is designed to produce 150 MW(e) and an 
equivalent 200 MW(e) for seawater desalination. The 
plant was scheduled for criticality?’ at the end of 
1971 and should be the first of the generation of 
high-flux intermediate-size [250:' to 500 MW¢(e)] 
reactors to operate at power. 

The core design of BN-350 (like that of the test 
reactor BOR-60) was frozen without allowance for 
stainless-steel swelling, although the design schedule for 
BN-600 (the next in the Russian series of LMFBR 
power plants) allows ample time for such provisions to 
be made within the constraints imposed by a fixed-size 
vessel.”7° BN-350 uses the passive type of core restraint 
provided by the annulus of radial blanket and shield 
assemblies cantilevered from long inlet nozzles sup- 
ported at two points in a conventional double grid 
plate offering hydraulic hold-down. Spacer pads, which 
occupy 0.059in. of the nominal 0.079-in. inter- 
assembly gap, are apparently located on the fueled 
length of the duct. These similarities with EBR-II 
might be expected to promote similar response to 
swelling and creep, except that the dissimilarities in 
spacer-pad design, shown7° in Fig. 16, could produce 
two important differences.® 

First, the 0.020-in. clearances in the spacer-pad 
plane (vs. 0.002in. in EBR-II) could produce a 
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Fig. 16 BN-350 core-assembly pattern. 
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significant positive structural-reactivity coefficient over 
part of the power range during early core lifetime. This 
results from thermal bowing and bowing due to 
differential swelling, producing top-end contact 
between adjacent rows of assemblies, hence inward 
deflection of the fueled span of the assemblies to the 
extent allowed by pad-plane clearances. After metal 
swelling absorbs this clearance (at a peak burnup of 
~30,000 MWd/metric ton based on U.S. swelling 
estimates®), the positive feedback potential is ex- 
changed for another threat, this one to fuel-handling 
feasibility. Note that in this design, which uses “one- 
sided” spacer pads (see Fig.16), duct dilation of 
20.020 in. below the spacer-pad plane will trap 
assemblies in the core and prevent their removal other 
than by forceful extrusion or sequential unloading, 
beginning with the outermost assemblies. 

An upper limit on burnup capability, based on 
U. S. swelling estimates,° is ~70,000 MWd/metric ton. 
This assumes utilization of the full 0.079-in. gap 
accommodating the dilation produced by ~6% 
volumetric swelling. Burnup capability of the early 
cores is calculated to be only ~46,000 MWd/metric 
ton, owing to fission-gas pressure.7° (Note, in Table 1, 
that the BN-350 fuel element has a negligible fission- 
gas plenum.) Assuming that swelling limitations permit, 
later cores may be extended to higher burnup with the 
BN-600 fuel-element design featuring integral axial 
blankets and a long fission-gas plenum below the core. 


Prototype Fast Reactor*®*72°7% 


In the patented’* PFR core-restraint concept, each 
fuel assembly, in groups of six, leans on a central 
structural member externally similar to the fuel- 
assembly duct. In most cases the “leaning posts” also 
serve as guide tubes for control or safety rods. This 
design approach, therefore, produces a characteristic 
rosette pattern of fuel and absorber assemblies, as 
shown”? in Fig. 17. 

Contact between the assembly and its assigned 
leaning post is brought about by an out-of-plumb 
alignment of two nozzle-support points. A bending bar, 
concentrically located in the fuel-assembly nozzle 
(“spike” in PFR nomenclature), deflects elastically 
under the force needed to insert the nozzle into its 
receiving socket,* producing a contact force of at least 
50 lb between duct and leaning-post spacer pads 
located well below the core. Among the relatively few 





*It is intended that this force be no more than the 
assembly’s deadweight.” * 
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Fig. 17 Core plan of the PFR. 


PFR design alterations resulting from the discovery of 
swelling, the spacer-pad location was chosen as a 
compromise between minimizing swelling at the pads 
(by keeping them well below the core) and minimizing 
flow-induced vibration (by keeping the support points 
as high as possible). To allow this location to be 
changed in the future, possibly as a result of better 
swelling data, the designers have made the leaning posts 
more readily removable. The finite life of the leaning 
post itself, due to swelling effects, also prompted this 
design alteration.”° 

Swelling-induced elongation of fuel assemblies is 
apparently not restricted, and dilation accommodation 
is relatively generous* for a reactor designed before 
swelling was recognized. Swelling-induced bowing of 
core assemblies is the greatest threat to the success of 





*The 0.16-in. gap in the 5.7-in. lattice spacing allows 
nominal dilation accommodation for a peak swelling of 8.7% 
AV/V, sufficient for the estimated swelling of 5% in PFR.?! 
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this design concept as a viable engineering solution to 
the swelling problem. Without restraining stresses, 
radiation-enhanced creep cannot play a correcting role 
to the effects of swelling. Periodic rotation of fuel 
assemblies, and possibly of leaning posts as well, is the 
operational “‘fall-back” position if bowing becomes 
excessive. The required frequency and operational 
practicality of duct rotation can be meaningfully 
assessed only through reactor operation, which is 
scheduled to begin (initial criticality)’ at the end of 
1972. 

The PFR design strategy places a high premium on 
the development of a duct material with low-swelling 
tendencies. The current British favorite is Nimonic 
PE16 (Ref. 21). If a low-swelling duct material is 
found, the absence of restraint (eliminating radiation- 
enhanced creep as an important determinant of core 
mechanical behavior) would be a design asset to go 
with the many other attributes of this approach.® 
Phenix 7° 7° 

This French LMFBR prototype uses the pot-type 
layout (like the PFR), is rated at 250 MW(e), and is 
scheduled for criticality in early 1973 (Ref. 71). The 
core design is an extrapolation of Rapsodie, the highly 
successful mixed-oxide-fueled test reactor operated at 
20 and 24 MW(t) from 1967 to 1970 and more 
recently at 40 MW(t) with the replacement 
“Fortissimo” core. 

With respect to core restraint, Rapsodie and Phenix 
are similar to EBR-II. That is: 

1. Reactor assemblies are cantilevered from a two- 
point nozzle support employing a double grid plate 
with hydraulic hold-down. 

2. Unrestrained top ends of reactor assemblies 
“flower outward,” in response to thermal gradients 
produced by the characteristic radial flux profile, until 
top-end interference is met. 

3. Lateral support is provided by the longitudinal 
stiffness of many rows of “‘straight”’ assemblies radially 
beyond the core. 


The Phenix core consists of 103 hexagonal fuel 
assemblies 4.87 in. across flats and six control assem- 
blies in guide tubes of the same dimensions. The 
high-pressure (110-psi) plenum also serves 90 radial 
blanket assemblies, 216 steel-reflector assemblies, and 
41 storage positions. Radially beyond the cylindrical 
grid box, comprising the high-pressure plenum, the 
core load support, and the two-point lateral nozzle 
support, is a “false grid.” This false grid provides 
vertical orientation, although no forced convection, to 
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937 lateral shield assemblies, which theoretically are 
also part of the lateral restraint system. These contain 
graphite, boron graphite, or steel rods and have a 
circular cross section either 4.88 or 6.22in. in 
diameter. 

The interassembly gap is 0.218 in., offering a 
2.63% linear dilation accommodation, equivalent to a 
maximum of 7.9% volumetric swelling. The spacer 
pads, as in EBR-II and Rapsodie, are on the mid-flats 
and formed out of the 316 stainless-steel duct. Origi- 
nally, the pads were to be placed at or near the core 
midplane. One of the few concessions to the swelling 
phenomenon* has been to relocate these pads,*? 
probably to above the core where they will result in a 
small positive bowing coefficient well overridden by 
the Doppler coefficient. 

The Phenix fuel elements are designed for a peak 
burnup of 100,000 MWd/metric ton; the minimum 
expectation is half of that.77 Irrespective of fuel 
performance, the operation of Phenix will provide the 
earliest indication of whether the EBR-II-type core- 
restraint approach, so successful in EBR-II and 
Rapsodie, is practical for commercial LMFBR power 
plants. 


SNR®°°8! 


This 300-MW(e) LMFBR station is being designed 
by an industrial consortium that includes firms in West 
Germany, Belgium, and the Netherlands. Construction 
of this loop-type prototype is scheduled to begin in 
1972 on its site in West Germany. In contrast with the 
other European prototypes, the timing of this project 
has permitted deliberate design steps to offset the 
effects of swelling and creep. At this writing, however, 
those steps have not been publicly revealed in detail. 

The 150 hexagonal fuel assemblies are relatively 
small (~4.3 in. across flats) and contain 169 elements 
supported by a system of intermittent grids. The 
assemblies are freestanding and cantilevered from a 
double grid plate incorporating hydraulic hold-down. 
An overhead instrument plate, hung from the top 
shield, provides an emergency stop against lifting of 
assemblies but does not interfere with swelling-induced 
elongation. Dilation is accommodated by interassembly 
gaps. A value of 8% has been mentioned?°’®? as the 
maximum expected volumetric swelling of core stain- 
less steel. 





*Other important modifications have been made to the 
fuel-handling system to tolerate distorted (bowed) assemblies 
and to facilitate periodic rotation if this should prove 
necessary .*1°78 
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Spacer pads are attached to all core and blanket 
assemblies, apparently outside the high-swelling region. 
Outside the radial blanket, which averages five rows in 
thickness, and at elevations matching the spacer pads, 
are ferritic steel rings that provide core radial restraint. 
Owing to differential thermal expansion between fer- 
ritic and austenitic steel, this passive system will 
produce fuel-handling clearances under isothermal reac- 
tor conditions and peripheral restraint against thermal 
and swelling-induced bowing at operating tempera- 
tures. As in the other European designs (PFR and 
Phenix), fuel-assembly bowing is reduced through the 
use of upstream plenums for fission-gas storage within 
the fuel elements. 


FFTF 


Most of the open literature on the effects of 
swelling and creep on reactor components, and the 
possibilities for mitigation of these effects through 
core-restraint design, has been produced by the FFTF 
project. Initially a literature review on the subject of 
fast reactor core restraint was conducted.** When 
early swelling correlations predicted that irreversible 
differential elongation of edge-core assemblies would 
produce several inches of free top-end deflection at the 
burnup and temperature targets of the Fast Test 
Reactor (FTR),** previous core-restraint design effort 
was seen as having limited relevance. A peripherally 
clamped core was chosen as the reference-design 
solution, and an extensive engineering-mechanics effort 
was inaugurated. 


At first, only the effects of swelling were ana- 
lyzed.85°8® Dilation and bowing of core assemblies, as 
free bodies and in an idealized r,z plane of the core, 
were computed based on a procession of swelling 
predictions. These predictions have trended downward 
from 1969 to the present, at least for 20% cold-worked 
316 stainless steel, which had been chosen for cladding 
and ducts. 

The recognition that radiation-enhanced creep 
would be an important determinant of FTR mechani- 
cal behavior stimulated the development of analytical 
tools not formerly required for fast reactor design. 
Responding to this need, Sutherland et al. produced 
the finite-element codes CRASIB®”’ and AXICRP.®*® 
CRASIB is specifically oriented to statically indeter- 
minate beams of arbitrary cross section. The code 
calculates time-dependent loads, stresses, and shapes of 
ducts using any functional representation for swelling 
and creep, although interassembly mechanical inter- 
actions of the type evaluated by the BOW-V code? 
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cannot be simultaneously solved. AXICRP is more 
general than CRASIB in that it can evaluate creep 
effects in plane-strain, plane-stress, and axisymmetric- 
bodies-of-revolution problems. Of the open-literature 
codes, cross-sectional distortion of hexagonal ducts can 
best be analyzed with AXICRP. In addition to these, 
several private and proprietary finite-element codes are 
known to exist in the United States and Europe for 
estimating the time-dependent effects of swelling and 
creep on core components. 


Thus evolved the capability to analyze the complex 
interacting effects of differential thermal expansion, 
differential swelling, and radiation-enhanced creep in a 
reactor with positive core restraint. During this evolu- 
tion the FTR fuel-assembly design had become rela- 
tively firm. Figure 18 shows the main features of 
core-assembly ducts. Not shown is the interassembly 
gap of 0.140in., allowing a 3% linear dilation of 
adjacent ducts before contact is made other than at the 
two above-core planes of load pads. Analyses of the 
response of these ducts to swelling, creep, and thermal 
gradients have been published** and are still in 
progress. 


Perhaps because of the many contributions of the 
FFTF project toward a coordinated engineering attack 
on the problems of dimensional instabilities of stainless 
steel in fast reactor cores, the peripherally actuated 
core clamp approach is currently the favored design 
strategy in the United States (see below). Concerned 
with the reliability of these complex active clamps, 
however, the AEC has sponsored a study to evaluate 
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the feasibility of a simpler passive restraint system 
patterned as much as possible after the unexpectedly 
successful EBR-II design approach. 


U.S. Demonstration Plant Designs 


Most of the design descriptions and analyses 
pertaining to U.S. demonstration plant designs have 
not yet been widely distributed, and the designs are 
evolving rapidly as of this writing. There is, however, 
sufficient open literature on each plant to provide a 
general understanding of how swelling and creep have 
been factored into the respective designs. The follow- 
ing synopses are gleaned from the open literature 
available as of December 1971. 


Atomics International.*®’®°:°° At 500 MW(e) the 
Al plant is the largest of the “first-generation” LMFBR 
demonstration plants being designed under the current 
U.S. breeder program. The AI plant also has the tallest 
core (44 in.) and the largest duct (5.57 in. across flats), 
which are important variables in the design of the 
active-clamp system of core radial restraint, which is 
the focal point of the design strategy to combat the 
effects of dimensional instability of stainless steel. Use 
of the active clamp with above-core and top-end load 
pads makes this design resemble the FTR to a greater 
extent than the others in this group. The active clamp 
was selected to achieve the dual benefits of: (1) large 
clearances to facilitate fuei handling and (2) opera- 
tional tightness at the load-pad planes to minimize 
reactivity anomalies, ensure a negative power coeffi- 
cient, and allow reproducible core performance. 
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Fig. 18 Driver-fuel-assembly duct in the FTR. 
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The 0.3-in. interassembly gap was sized to preclude 
assembly-to-assembly contact (other than at the pad 
planes) by dilation and bowing. With AI’s current 
swelling and creep estimates (nominals), this gap is 
sufficient to prevent duct—duct mechanical interaction 
for the entire 2-year (four-cycle) residence time, during 
which the peak total fluence will reach 4x 107? 
neutrons/cm?. Nominal swelling under these con- 
ditions is 4.2% AV/V. At the 30 upper limit of swelling 
(6.6% AV/V), in conjunction with nominal to mini- 
mum creep, assembly rotation is required to avoid 
contact of ducts near the core periphery, where bowing 
tendencies are greatest. The maximum swelling- 
accommodation potential of the design is given as 15% 
AV/V, based on the dilation allowance for the peak- 
fluence (central) assembly. 


A half-scale full-core mockup including the six 
clamping segments has been built and operated in 
room-temperature air to study the mechanical response 
of the core to peripheral clamping forces with straight 
assemblies and simulated bowed assemblies. On the 
basis of calculation, the required clamping force is 
33,000 Ib per segment, considerably higher than the 
FFTF, due primarily to the increased duct diameter. 
The clamp is actuated by a system of hydraulic 
cylinders using sodium as the working fluid. 


Combustion Engineering.?''?? The only current 
U. S. LMFBR design advertised as a “second-generation 
demonstration plant,” this 500-MW(e) design uses 
high-performance mixed-carbide fuel that, with a peak 
linear power of 32 kW/ft (19 kW/ft av.), permits a 
more compact core than the AI design with the same 
electrical rating.* 


For design purposes, peak metal swelling was taken 
to be 10% AV/V at a peak discharge burnup of 
100,000 MWd/metric ton, corresponding to a peak 
fluence of 2.3 x 107? neutrons/cm? (£ >0.1 MeV). 
Dilation is accommodated by a ~0.3-in. gap between 
hexagonal assemblies measuring 5.932 in. across out- 
side flats. This contrasts with 0.050-in. gaps for 
fuel-handling allowance used in previous Combustion 
Engineering designs.°’'* For the more challenging 
design task of accommodating bowing within the 
constraints of acceptable reactivity effects, four alter- 
natives were evaluated, that is (1) active peripheral core 
clamps and an articulated assembly, (2) decreased duct 
width to eliminate the articulation, (3) “unrestrained” 
core that permits free mechanical interaction of ducts 





*Height, 36 vs. 44 in.; diameter, 74 vs. 76 in. 
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and some lateral movement, and (4) development of a 
suitable low-swelling duct material. 

Design options were affected by the following 
ground rules: (1) the bowing coefficient must be 
negative, (2) no duct—duct mechanical interaction 
except at load pads, (3) no load pads in regions of 
significant swelling, and (4) assemblies must be re- 
strained at a minimum of one level above the grid 
plate. Under the constraints of these requirements and 
certain economic criteria, the first alternative was 
chosen as the basic concept for swelling—creep ac- 
commodation. 

The six mechanically actuated compliant peripheral 
clamps apply a constant force on load pads located on 
the ducts at elevations 6 and 48 in. above the top of 
the core. Bottom support of assemblies uses a conven- 
tional double grid plate and hydraulic hold-down. 
Because of the large duct diameter and the resultant 
high stresses that the core-restraint system would 
impose on a conventional duct, it was deemed neces- 
sary to articulate the fuel assemblies with a double 
joint between the upper grid plate and the bottom of 
the fuel elements. The articulation is essentially a 
hollow open-ended dumbbell supported in two sockets. 
Radial blanket assemblies do not have the articulation 
joints, but instead have a floating-collar type of upper 
load pad to maintain stresses at allowable values by 
permitting 0.25in. of unrestrained top-end radial 
movement. This type of load pad was avoided in the 
core because it tends to produce a positive bowing 
coefficient which the lower-end articulation scheme 
does not. 


General Electric.?*'°*? The 350-MW(e) GE demon- 
stration plant has the smallest duct size (3.849 in. 
across hex flats) of any plant in the present generation 
of U.S. LMFBR designs. With the resulting large 
number of fuel assemblies (276 in two enrichment 
zones in the target core), the unit worth is low— 
prototypic of larger assemblies in proportionately 
larger plants. Although safety is the prime motivation, 
clamping loads and duct stresses are greatly reduced by 
selection of this duct size, other factors being equal. 
The interassembly gap in the GE design is 0.160 in. 

The design strategy for core restraint involves 
clamping, but it differs significantly from the active 
peripheral clamps used in the FFTF, AI, and Combus- 
tion Engineering reactor concepts and the peripheral 
clamping system designed by GE for SEFOR.?* Three 
interacting systems are involved: 

1. Banding ring segments positioned to form a hard 
outer boundary. 
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2. Spring clamps within this outer boundary and 
surrounding the radial blanket. 

3. Internal core clamps that are expanded to take 
up the duct-spacing clearance. 


The distinguishing feature is the third system. The 
internal clamps are distributed in a uniform pattern 
with each clamp surrounded by six fuel assemblies. In 
the core the clamp mechanisms are mounted in 
absorber-assembly guide tubes, thus distributing the 
absorber positions in a symmetry resembling that of 
the PFR. 


The reference core-restraint design has, as one of its 
aims, increased flexibility for accommodating the wide 
range of possible materials behavior. With compliance 
to absorb swelling, thermal expansion, and dimensional 
tolerance provided at the blanket periphery, the 
mechanical response of the core to temperature in- 
crease, for example, should be as a single structure with 
calculable and reproducible neutronic effects. The 
actual clamping forces, however, are generated locally 
and need not be transmitted across many rows of 
assemblies to achieve the intended planes of tightness 
at the load-pad elevations. The design should, there- 
fore, be tolerant of higher friction coefficients at the 
load pads. 


Like SEFOR, the GE demonstration plant uses the 
open-reactor, hot-cell refueling approach. Among other 
reasons for selecting this concept, the feasibility of 
remapping fuel-assembly handling heads, should they 
become misplaced due to swelling and creep effects, 
has been mentioned.?® 


Gulf General Atomic.*7’°* The only gas-cooled 
fast reactor (GCFR) concept to be reviewed here, the 
300-MW(e) GGA (Gulf General Atomic) demonstration 
plant is also the smallest in the current generation of 
U.S. breeder reactor designs. Owing in part to the 
differences between sodium and helium as reactor 
coolants, this design shows some major differences in 
core support and hence in the gross strategy of 
accommodating swelling and creep. Along with these 
differences, there are many similarities with respect to 
core materials, temperatures, and fuel-element designs. 
The similarities are often deliberate to capitalize on the 
LMFBR research and development program in areas of 
fuels and materials development and reactor physics. 
An FFTF correlation for 20% cold-worked 316 stain- 
less steel?* has been used to estimate the amount of 
swelling in the GGA plant. The maximum swelling, at 
GCFR demonstration plant conditions, is estimated at 
10% AV/V for a fluence of 1.8 x 10?* neutrons/cm? 
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(E>0.1 MeV), corresponding to a maximum fuel 
burnup of 100,000 MWd/metric ton. 

The reactor loading contains 118 fuel and 93 radial 
blanket assemblies featuring very large (6.5 in. across 
flats) hexagonal ducts. They are supported from, and 
rigidly attached to, a top-mounted single grid plate. 
The assemblies project downward from this grid plate, 
and fuel handling is performed from beneath the core 
(control rods are driven from above). In this “free- 
hanging” array of core assemblies, as in the freestand- 
ing array in the PFR, duct bowing due to differential 
thermal expansion and differential swelling is un- 
suppressed. Therefore radiation-enhanced creep does 
not have the importance it has in restrained-core 
concepts. Instead, duct dilation and bowing are accom- 
modated by a 0.25-in. interassembly gap that extends 
along the entire 10-ft length of duct, and swelling- 
induced bowing will be held to acceptable values by 
periodic rotation of assemblies. It is anticipated that 
annual rotation, at times coinciding with partial-core 
refueling, will reduce duct—duct interferences to an 
acceptably small degree. Both 180 and 120° rotation 
increments have been evaluated, with no significant 
difference in the calculated degree of interference. 


Westinghouse.°°®°* This 330-MW(e) (gross) plant 
was designed to accommodate 20% maximum stain- 
less-steel swelling (by volume). All structural material 
in the core, including cladding, spacer grids, and ducts, 
is 20% cold-worked 316 stainless steel; and the 
maximum burnup, based on a 6-month refueling cycle, 
is 103,000 MWd/ metric ton. 

The reference core-support and core-restraint con- 
cepts were chosen after a trade-off study that evaluated 
various methods of accommodating swelling and creep 
phenomena. Both hydraulic and mechanical hold-down 
designs were studied. Mechanical over-the-core hold- 
down was chosen, a selection that runs counter to the 
present trend of U.S. and European LMFBR design 
practice. Simplicity and improved predictability are 
cited as the reasons for this choice. With the top ends 
of assemblies securely fixed in a massive hold-down 
structure, and the lower nozzles captured in receptacles 
that offer two-point positioning, radial restraint can be 
somewhat less elaborate than designs utilizing active 
clamps. In the Westinghouse concept, radial restraint is 
assisted by a single row (72) of longitudinally flexing 
“restraint-spring” assemblies surrounding the core 
system of 367 hexagonal assemblies.* Beyond that, a 


*The design consists of 153 fuel, 16 control (based on a 
6-month refueling cycle), 132 radial blanket, and 66 radial 
reflector assemblies. 
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nonyielding “former plate” provides a fixed outer 
diameter and limits total radial free space. The radial 
restraint strategy also requires rotation of some core 
and blanket assemblies every 12 months, based on 
existing data. 

As in Fermi 1, axial thermal expansion and 
swelling-induced elongation of assemblies are accom- 
modated by a nozzle spring, in this case made of 
Inconel X-750. Dilation is absorbed by an interas- 
sembly gap that is unobstructed over the length of duct 
that will experience significant swelling. There are two 
levels of load pads above the core and one below. Pads 
are located at the same elevation on all assemblies from 
core center to the peripheral restraint system. The 
protrusion of the pads beyond the duct varies with the 
type of assembly and the axial elevation, but dimen- 
sions across pads are always slightly less than the cold 
pitch established by the spacing of holes in the support 
plate and hold-down plate. If this embodiment pro- 
duces any positive bowing effects, they are over- 
shadowed by negative constituents of the power 
coefficient. The overall power coefficient is monotoni- 
cally negative over the whole power range, having an 
average value of —0.29 cent/MW(t) and a value of 
—0.13 cent/MW(t) at full power. 

Although the reference design requires semiannual 
refueling, with about 30 fuel and 13 radial blanket 
assemblies replaced each time, the ultimate objective is 
annual refueling, although this will require more 
favorable materials irradiation data than are now in 
hand, or mechanical releases in the passive restraint 
assemblies. Uneven swelling of core components is said 
to be the greatest mechanical disadvantage of long 
refueling intervals. 


EFFECTS ON FAST REACTOR 
ECONOMICS 


The general view of fast reactor designers is that 
the need to accommodate swelling and creep will 
increase the cost of energy production, but not by so 
much as to undermine the role of the breeder reactor 
as a major future energy source for the world. 

Capital costs are subject to increase from (1) larger 
reactor size resulting from space needed to accom- 
modate swelling in the core and (2) complex reactor 
internals not heretofore utilized, e.g., core clamps. 

The only component of the energy-cost increase 
calculable on first principles is the cost of the enlarged 
interassembly gap reflected through decreased breeding 
ratio and specific power. This cost has been calculated 
as ~0.08 mill/kWh for 15% swelling accom- 
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modation®’®® and proportionately less for lesser 


amounts of swelling. The upper boundary on fuel-cycle 
cost increase has been suggested as 0.35 mill/kWh, 
based on the worst-case assumption that fuel burnup 
would need to be reduced by half.® 

Several design approaches have been suggested that 
represent more radical departures from preswelling 
LMFBR design practice. All involve penalties in fuel- 
cycle costs which generally fall between the limits 
given above.° Some of these concepts, and the nature 
of the economic penalty, are given below. 


1. “Cold-Wrapper Tube’.®''°® The duct is main- 
tained in a temperature range below the swelling 
threshold, a task requiring a reactor-inlet sodium 
temperature of only ~600°F. The penalty arises from 
the less-than-optimum core-outlet temperature of prob- 
ably ~950°F. 

2. Refractory-Metal Core Structure. 
fractory metals, such as niobium or molybdenum, 
which apparently resist swelling at temperatures in the 
range of interest for fast reactors,'°* are substituted 
completely or strategically for stainless steel in the 
core. Penalties arise through adverse neutronic proper- 
ties of these materials and higher fabrication costs. 

3. Periodic Annealing.®’'°® Analogous to the 
method of coping with graphite swelling in some 
gas-cooled reactors, swollen core components are peri- 
odically shrunk to original size by a brief annealing 
treatment at >1400°F, which is effective in eliminating 
voids.'°” The penalty arises through reactor unavail- 
ability and/or increased on-site fuel inventory. 

4. Periodic Rotation.®*!°*® This concept is a cur- 
rently favored strategem to limit the extent of swell- 
ing-induced bowing. Loss-of-plant availability is trans- 
latable to an energy-cost penalty. 

5.?733U_Th Cycle.©'1°°*!!° Higher fission cross 
sections and the feasibility of breeding at intermediate 
neutron energies (E<0.1 MeV) greatly reduce the 
swelling/burnup ratio in the 723 U-—Th cycle relative to 
the ??°Pu—U cycle. The penalties are a substantial loss 
of breeding ratio and a heavy sacrifice of existing 
technology on which to design commercial breeders for 
the 1980s. 


6,101-104 Re- 


Returning to the (U,Pu)O,-fueled LMFBR with a 
stainless-steel core structure, it appears that the swell- 
ing phenomenon has thwarted the development of core 
features that could have produced a step increase in 
core-fuel fraction, with resultant increases in breeding 
ratio and specific power. Figure 19 shows® the effect 
of breeding ratio and specific power on fuel-cycle cost 
(other parameters constant) and a thumb-shaped zone 
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Fig. 19 Effects of specific power and breeding ratio on fuel- 
cycle costs and potential impact of swelling. 


defining these values in all LMFBR studies reported to 
date. Notice that the leading edge of the thumb, 
representing most favorable economics, is defined by 
two GE designs''' that featured nonfreestanding 
ducts and essentially no interassembly sodium gaps, 
and hence the highest possible fuel fraction. Since 
these features now seem impractical in view of the 
need to accommodate swelling, the portion of the 
thumb whose leading edge is defined by the AI’? and 
Babcock & Wilcox (B&W)'? 1000-MW(e) concepts 
now seems to be a more practical design regime. The 
B&W design had considerable subsequent study related 
to the impact of swelling and the evaluation of various 
modifications." 1? 

Taking the “physics” value of increased interas- 
sembly sodium gaps as the minimum cost of accom- 
modating swelling and a 50% burnup reduction as the 
maximum cost, a cost—benefit analysis for develop- 
ment of low-swelling alloys has been made.° The 
analysis indicated that the cost savings of having to 
accommodate only 5% AV/V in U.S. breeder reactors 
to be built over the period 1970—2020 (Ref. 113), 
rather than 15% AV/V, have a 1970 present worth of 
from $860,000,000 to $5,600,000,000. 


CONCLUSIONS 


Fast reactor designers have now had 4 years in 
which to evolve design strategies to cope with the 
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various problems posed by radiation-induced swelling 
and radiation-enhanced creep in core stainless steel. At 
this time, no worldwide consensus can be presented, 
since (1) the open literature on design strategies is 
almost entirely American and (2) no reactor specifi- 
cally designed for swelling and creep effects is even 
close to operation, and only one (the FFTF) is under 
construction. Certain trends are apparent in the open 
design literature, however. Those which seem note- 
worthy to this reviewer are: 


1. Despite the sophistication with which swelling 
and creep effects are analyzed and predicted in fuel 
elements, these phenomena have affected fuel-element 
designs very little. There is no clear trend to higher 
smeared densities, as might be justified by enhanced 
cladding dilation; nor have fuel-element designs be- 
come noticeably less conservative, as might result from 
a general resignation that swelling and creep effects, 
rather than intrinsic burnup capability of high-quality 
fuel elements, would limit core lifetime. 


2. Internal design features of fuel assemblies have 
also been relatively unaffected. The helical wire wrap 
still dominates the field of candidates for fuel-element 
support, but more built-in clearances or other schemes 
to produce radial compliance are used than in the 
““preswelling era.” 


3. Fuel-assembly dilation is accommodated by in- 
creased gaps between assemblies, with the gaps con- 
tinuous over the duct length where significant swelling 
is expected. Designers have abandoned the in-core pads 
featured for safety reasons in earlier and smaller fast 
reactors.* 


4. Swelling-induced elongation of core components 
is considered a minor problem, especially in the 
majority of designs where the use of hydraulic hold- 
down leaves the top ends of reactor assemblies free to 
grow. 

5. Restraint of swelling-induced bowing of reactor 
assemblies is favored over no restraint. Active periph- 
eral clamps, which generate comfortable refueling 
clearance, are the most popular current concept for 
providing this radial restraint. Reliable mechanisms and 
reproducible complete or near-complete compaction 
are system requirements under development. 





*A redeeming characteristic of large fast reactors, offset- 
ting their higher swelling/burnup ratio (Fig. 3), is their de- 
creased neutronic sensitivity to radial movement of fuel. This 
neutronic sensitivity, expressed as 6k/k per mil change in core 
radius, is inversely proportional to the core radius as a first 
approximation.’ '4 
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6.In the restrained-core concepts, uncertainties in 
swelling and creep combine in equal measure, accord- 
ing to the principles of superposition, to uncertainties 
in duct deformation during service. 


7. Although a metallurgical solution to the prob- 
lems of swelling and creep is widely desired by reactor 
designers, it is generally hoped and expected that this 
solution will be found within the Fe—Cr—Ni family of 
alloys, rather than in the refractory metals. 


8. Temperature and burnup goals for commercial 
LMFBRs have not been lowered due to swelling and 
creep phenomena, and the economic impact of these 
properties is considered tolerable even without a 
metallurgical breakthrough. 


As for the nearly completed high-flux reactors 
designed with little or no deliberate provision for 
swelling, Hafele appears to represent European opinion 
in saying that “6—8% AV/V can be lived with.”®? This 
is a plausible prediction in view of the successful 
operation of EBR-II with over 11% swelling of certain 
“permanent” core internals, although the reasons for 
this unexpected success are not well understood and 
could be unique to EBR-II. Prevailing metallurgical 
opinion is that cold-worked and/or stabilized stainless 
steels will not reach 6 to 8% AV/V until somewhat 
beyond 107? neutrons/cm?, an adequate achievement 
for the European demonstration plants. 


For the future, optimism or pessimism on the 
swelling—creep issue seems to pivot around materials 
data rather than design solutions. Some investigators 
are more optimistic about cold-worked stainless steel 
than given above, predicting as little as 5% AV/V for 
burnups of the order of 100,000 MWd/metric ton 
(Ref. 71), which is far below the estimates promul- 
gated 2 to 3 years ago. Significantly, this “interim 
metallurgical solution” has been found within the 
family of preferred commercial alloys for fast reactor 
construction. As a systematic search for improved 
alloys is conducted, assisted by larger and more 
prototypic irradiation facilities and an expanded theo- 
retical base, it appears probable that fast reactors 
designed near the turn of the century will be little 
affected by swelling. 


To put the swelling problem in perspective, 
Bishop?! observes that stainless steel heated from 
room temperature to ~1000°F undergoes a 3% volu- 
metric expansion. Therefore, if this recent optimism 
proves accurate, the swelling problem has been reduced 
to magnitudes of dimensional change with which 
designers of complex equipment are familiar. 
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